Plasma-Assisted Ignition in the Flow-Through Combustion Chamber
Esakov I.I. 1, Ravaev A.A. 1, Grachev L.P. 1, Volobuev I.A. 2
1 Moscow Radiotechnical institute of Russian Academy of Science 2 Saint-Petersburg National Research University of Information Technologies, Mechanics and Optics (ITMO University), Saint-Petersburg, Russian Federation
Abstract. The work is aimed to the problem of the organization of combustion in a high-speed flow of a fuel mixture using subcritical streamer microwave discharges. The problem is relevant for engines using combustion at a constant volume, in high-speed subsonic and supersonic flow of the fuel mixture, as well as in detonation engines with forced ignition. The goal is to investigate the possibility of igniting of a series of attached subcritical microwave discharges located at some distance from each other in a high-speed stream near metal or dielectric screen, and also to investigate the ability of these discharges to ignite a fuel mixture. This goal is achieved by the fact that numerical calculations of microwave discharges near the screen are performed at a flow velocity of 5 m / s to 500 m / s. Then, the calculations are repeated for the fuel mixture with a fuel excess ratio n from 0.55 to 1.5. The results of calculations of a mixture of propane with air are compared with experiments that were carried out on an apparatus with a wavelength of X = 12.3 cm with a microwave beam power of P = 1 kW and a microwave pulse duration of t = 0.2 sec. The significant result is that the deeply critical discharge connected to the resonator on the dielectric screen is not blown away by the high-speed flow and ignites the mixture. Importance - a microwave discharges ignite fuel over the entire speed range. The calculation results are confirmed by experiment.
Keywords: combustion chamber, detonation engine, constant volume combustion, plasma - associated combustion, supersonic flow combustion, microwave discharge, streamer discharge, nonequilibrium plasma.
DOI: 10.5281/zenodo.3562187
Combustia plasmatics intr-o camera de ardere prin flux Yesakov I.I. 1, Ravaev A.A.1, Grachev L.P. 1 Volobuev I.A.2
1Institutul de Inginerie Radio din Moscova al Academiei de Stiinte din Risia, Moscova, Federatia Rusa,
2Universitatea Nationala de Cercetari a Tehnologiilor Informational, Mecanica si Optica din Sankt Petersburg
(Universitatea ITMO), Sankt Petersburg, Federatia Rusa Rezumat. Lucrarea este dedicata problemei aprinderii intr-un flux de mare viteza a unui amestec de combustibil folosind descarcari strimer de microunde. Problema este relevanta pentru motoarele care utilizeaza combustia la un volum constant, in fluxul subsonic si supersonic de mare viteza al amestecului de combustibil, precum si in motoarele de detonare cu aprindere fortata. Scopul este de a cerceta posibilitatea de aprindere a descarcarilor de microunde subcritice atasate intr-un flux de mare viteza in apropierea unui ecran metalic sau dielectric, precum si de a investiga capacitatea acestor descarcari de a aprinde un amestec de combustibil. Acest obiectiv este atins prin faptul, ca se realizeaza calculele numerice ale descarcarilor de microunde in apropierea ecranului pentru vitezele de curgere intre 5 m / s si 500 m / s. Se calculeaza concentratia de electroni, temperatura gazului si temperatura electronilor. Se determina caracterul descarcarii privind faptul ca exista sau nu exista descarcarea. La pasul urmator se repeta calculele pentru amestecul de combustibil cu un raport de exces de combustibil n de la 0,55 pana la 1,5. Un amestec de gaz lampat cu aer este inlocuit cu un amestec model. Se cerceteaza de asemenea un amestec de propan cu aer. Rezultatele calculelor unui amestec de propan cu aer sunt comparate cu experimentele efectuate pe o instalatie cu o lungime de unda de X = 12,3 cm cu o putere fascicului de microunde de P = 1 kW si o durata a impulsului cu microunde t = 0,2 sec. Cel mai semnificativ rezultat stiintific consta in faptul, ca o descarcare subcritica profunda conectata la un initiator localizat pe un ecran dielectric nu este suflata de un flux de mare viteza si aprinde cu stabilitate amestecul. Semnificatia consta in faptul, ca o descarcare aprinsa pe un rezonator in apropierea unui ecran dielectric sau metalic nu este suflata de un curent cu o viteza cuprinsa intre 50 m / s si 500 m / s. Descarcarea la aceste viteze aprinde cu siguranta atat amestecului de combustibil bogat, cat si pe cel sarac. in apropierea limitei de concentratie „slaba" a combustiei se aprinde si arde numai combustibilul care a trecut direct prin descarcare. Rezultatele calculului sunt confirmate prin experiment.
Cuvinte-cheie: camera de ardere, motor de detonare, ardere in volum constanta, combustie plasmatica, ardere intr-un flux supersonic, descarcare de microunde, descarcare in flux, plasma fara echilibru.
© EcaKOB H.H., PaBaeB A.A., rpaneB H.n., Bo^o6yeB H.A. 2019
Плазменное горение в проточной камере сгорания Есаков И.И.1, Раваев А.А.1, Грачев Л.П.1, Волобуев И.А.3
Московский радиотехнический институт российской академии наук, Москва, Российская Федерация, 2 Санкт-Петербургский национальный исследовательский университет информационных технологий,
механики и оптики (Университет ИТМО), Санкт-Петербург, Российская Федерация Аннотация. Работа посвящена проблеме поджигания в скоростном потоке топливной смеси при помощи подкритических стримерных СВЧ разрядов. Проблема актуальна для двигателей, использующих горение при постоянном объеме, в скоростном дозвуковом и сверхзвуковом потоке топливной смеси, а также в детонационных двигателях с принудительным поджиганием. Цель - исследовать возможность зажигания присоединенных подкритических СВЧ разрядов в скоростном потоке вблизи металлического или диэлектрического экрана, а также исследовать способность этих разрядов поджигать топливную смесь. Поставленная цель достигается тем, что выполняются численные расчеты СВЧ разрядов вблизи экрана при скорости потока от 5 м/с до 500 м/с. В расчетах вычисляется концентрация электронов, температура газа и электронная температура. Определяется сдувается разряд или нет. Затем расчеты повторяются для топливной смеси с коэффициентом избытка горючего п от 0,55 до 1,5. Смесь керосина с воздухом заменяется модельной смесью, также исследуется смесь пропана с воздухом. Результаты расчетов смеси пропана с воздухом сравниваются с экспериментами, которые проводились на установке с длиной волны X = 12,3 см при мощности СВЧ пучка Р = 1 кВт и длительности СВЧ импульсов т = 0,2 сек. Наиболее существенный научный результат состоит в том, что глубокоподкритический разряд, присоединенный к инициатору, расположенному на диэлектрическом экране, не сдувается скоростным потоком и уверенно поджигает смесь. Значимость - показано, что разряд, зажженный на резонаторе вблизи диэлектрического или металлического экрана не сдувается потоком со скоростью в пределах от 50 м/с до 500 м/с. Разряд при этих скоростях уверенно поджигает как богатую, так и бедную топливную смесь. Вблизи «бедного» концентрационного предела горения поджигается и горит только топливо, непосредственно прошедшее через разряд. Результаты расчетов подтверждены экспериментом.
Ключевые слова: камера сгорания, детонационный двигатель, горение при постоянном объеме, плазменное горение, горение в сверхзвуковом потоке, СВЧ-разряд, стримерный разряд, неравновесная плазма.
ВВЕДЕНИЕ
Актуальна задача повышения эффективности сгорания топлива в относительно скоростных потоках. В настоящей работе рассматривается возможность создания и применения неравновесного СВЧ разряда для поджига и интенсификации горения топлив в высокоскоростных потоках [1, 2]. Создание плазмы разрядами различного вида, а также различные подходы к организации горения в присутствии холодной неравноценной плазмы рассмотрены в обзорных работах [3], [4], [5]. В типичных случаях в камере сгорания авиационного воздушно-реактивного двигателя скорость потока воздуха на входе составляет порядка 100-200 м/с. Возникает задача стабилизации горения топливно-воздушной смеси, повышения эффективности смешения топлива с воздухом, которая решается путем введения в поток фронтовых устройств с плохо-обтекаемыми телами, за которыми образуются области рециркуляции топливной смеси. Такие устройства хорошо отработаны, но они создают внутреннее газодинамическое и гидравлическое сопротивление, что снижает коэффициент полезного действия (КПД) двигателя.
Изучаются иные способы стабилизации горения скоростного потока, например, помещением в поток линейки СВЧ разрядов. В работе [6] показано, что подкритический стри-мерный разряд в отличие от диффузного разряда равного размера эффективно поглощает энергию СВЧ излучения, быстро разогревается и поджигает топливную смесь при достаточно больших скоростях потока.
При этом энергоэффективность такого разряда существенно выше иных способов [7] создания плазмоидов и стабилизации с их помощью горения.
КПД газотурбинного двигателя или энергетической установки увеличивается по мере снижения доли топлива в составе сжигаемой топливно-воздушной смеси [8], поэтому современной тенденцией является переход к сжиганию всё более бедных топливных смесей вплоть до достижения режима так называемого «холодного» горения, когда выделяется не более 20% энергии, запасенной в топливе. В предыдущих работах авторов было показано, что в среде холодной неравноценной плазмы возможно горение с концентрацией топлива ниже «бедного» предела горения, см. также работу [9].
По мере повышения скорости полета летательного аппарата увеличивается и оптимальная скорость потока сжигаемой топливной смеси. Сжигание топлива в скоростных потоках сопровождается значительными трудностями. Начиная со скорости полета летательного аппарата, соответствующей числу Маха М=6 (число Маха - отношение скорости к местной скорости звука), сжигание должно происходить в сверхзвуковом потоке топливной смеси. Такую смесь необходимо предварительно сжать в изоэнтропических волнах сжатия [10] или оптимальных ударно-волновых структурах [11].
Рассматриваются также различные концепции двигателей, работающих по термодинамическому циклу горения при постоянном объеме или детонационного горения [12]. При этом может поджигаться сверхзвуковой поток, набегающий на стационарную ударно-волновую структуру, либо по покоящейся топливной смеси распространяться детонационная волна [13]. Такие задачи отличается друг от друга только способом вычисления температуры торможения за фронтом горения. Потенциально детонационные двигатели на 20-25% превосходят по КПД обычные газотурбинные двигатели, работающие по циклу Брайтона горения при постоянном давления, поэтому к ним последние 30 лет приковано большое внимание. Особенно популярна концепция ротационного детонационного двигателя (rotational detonation engine - RDE). Однако еще никому не удалось создать RDE, превосходящий по КПД двигатель с термодинамическим циклом Брайтона. Причина, как показано в работе [14], заключается в том, что в RDE образуется не сплошная прямая детонационная волна, а более сложная система,
состоящая из тройных конфигураций ударных волн. Оптимальные с точки зрения повышения полного давления тройные конфигурации [15,16] в КОЕ реализовать не удается. Наоборот, после возникновения детонации ударно-волновая структура в процессе своего движения по кольцевому зазору камеры сгорания КОЕ подстраивается таким образом, что количество выделяемой энергии детонирующей топливной смеси стремится к минимуму, а не максимуму.
Представляется привлекательным и в этом случае использовать подвижный фронт горения в неподвижном потоке топливной смеси или наоборот неподвижный фронт в скоростном потоке, создаваемый, например, разветвленным стримерным разрядом. Такой метод называется слабой детонацией и может быть реализован при помощи подкритического стримерного разряда [17] (рис.1). Предыдущие работы коллектива показали, что это возможно, но затраты энергии на создание стри-мерного разряда с развитой структурой при этом слишком велики, что делает такую схему нерациональной для применения в камерах сгорания энергетических установок.
Цель работы - исследовать возможность зажигания присоединенных подкритических СВЧ разрядов в скоростном потоке вблизи металлического или диэлектрического экрана, а также способность этих разрядов поджигать топливную смесь.
Ниже излагаются результаты исследования в соответствии со сформулированной целью, наиболее существенный научный результат состоит в том, что глубокоподкрити-ческий разряд, присоединенный к инициатору,
а) б) в) г)
Рис.1. Стримерный СВЧ разряд в свободном пространстве (а), на поверхности из текстолита (б) и
стекла (в), на капроновой сетке (г).1
I. МОДЕЛИ И МЕТОДЫ
расположенному на диэлектрическом экране, не сдувается скоростным потоком и уверенно поджигает смесь.
Полная нестационарная математическая модель формирования стримерного разряда в поле электромагнитного импульса СВЧ излучения, а также поджигания этим разрядом
1 Appendix 1
67
топливнои смеси и её последующего горения весьма сложна. В рамках настоящего исследования сделано несколько упрощении.
Для расчета химической кинетики горения смеси пропан - воздух используется так называемая квази-глобальная модель горения [18], включающая 12 компонентов (СНв, О, Н, ОН, О2, N2, Н2, СО, НО, НО2, Н2О2, СО2), одну глобальную реакцию
СН + Ю ^ 3СО + 4Я2
(1)
и элементарные реакции, взятые из детального химического механизма [19]. Скорость реакции (1) определяется по формуле
ю = А ■ ехр(-Е / ЯТ)[СН Г [О ]Р. (2)
В формуле (2) эмпирические константы равны: А = 1.5-1012, Е = 30 ккал/моль, а = 0.1, Р = 1.65 ; Я - газовая постоянная, Т - температура в градусах Кельвина. В настоящее время жидк керосин является основой широко использующихся авиационных и ракетных реактивных топлив. Основной трудностью является моделирование первой стадии процессов - переход керосина и продуктов его сгорания в газовую фазу. Часто достаточно иметь возможность получения газового состава, эквивалентного составу сгоревшего керосина. Эта возможность может быть реализована при моделировании сгорания жидкого авиационного керосина сгоранием смеси газообразных топлив, при котором образуется газ такого же состава и с такими же параметрами, как и при сгорании авиационного керосина в газовой фазе.
Горение паров газообразного керосина в воздухе в рамках настоящей работы заменяется в расчетах горением модельной
смеси метана СН4 и ацетилена СИ2, как это рекомендовано в работе [20]. Для горения также используется квазиглобальная модель горения.
Поскольку при горении топлива температура может существенно превышать 1000К в рамках настоящей работы используется формула, в которой учитывается зависимость показателя адиабаты от температуры [21]
у(Т) = 1 +
У, -1
1 + (У, -1)
Т
в/Т
(1 - ев/Т )2
. (3)
В формуле (3) у - показатель адиабаты при нормальных условиях, а в - энергетическая постоянная, равная 3056,4 К. При больших температурах выражение (3) имеет свой предел. Например, при у = 1.4 зависимость
у(Т) стремится к 1.286.
Как показали предыдущие исследования [22], глубоко подкритический разряд имеет такую же природу, что и стримерный подкри-тический разряд с развитой пространственной структурой, но стримеры не выходят за пределы плазменного образования. Это позволяет использовать для расчета модель формирования стримеров и поглощения ими энергии электромагнитного СВЧ излучения [23].
Показано, что при внешней схожести с диффузным разрядом, присоединенный глубоко подкритический стримерный разряд эффективно нагревается СВЧ излучением, а диффузный разряд СВЧ излучение не поглощает и остается холодным [24]. В условиях резонанса в проводнике возникает большой индуцированный ток, а в окрестности инициатора происходит ионизация газа, образуются плазменные каналы, сильно нагревающие окружающий газ. Если перегретый канал не может вырваться за пределы диффузного облака, то разряд остается присоединенным. Если напряженность достаточно большая, то из перегретого канала развивается пространственная древовидная структура, показанная на рис.1. В настоящей работе рассматривается случай, когда стример остается внутри плазменного образования.
Резонатор располагается на некотором расстоянии от проводящего экрана (рис.2-а), а в случае диэлектрического экрана - прямо на его поверхности со стороны СВЧ излучения (рис.2-б) или с обратной стороны (рис.2-в).
Резонатор - металлическая проволока или полоска металла. Чем меньше радиус закон-цовки а и чем ближе резонатор к экрану, тем больше коэффициент усиления поля КЕ .
В случае диэлектрического экрана резонатор можно разместить непосредственно на поверхности (см. рис.2-б, в). Максимум интенсивности поля в этом случае будет в пространстве, образованном поверхностью экрана и радиусной законцовкой резонатора.
Хотя слой диэлектрика ослабляет поле и коэффициент усиления КЕ в случае, показанном на рис. 2-в ниже, чем в случае рис.2-б, такой вариант для камеры сгорания актуален, поскольку позволяет отделить волновод и источник СВЧ излучения от зоны горения.
Ранее в работах авторов были получены численные результаты расчета К и их аппроксимация для резонатора в свободном
пространстве. Выполнен расчет скорости распространения стримера, его диаметра.
Показано расчетным путем и доказано экспериментально, что стримерный разряд с развитой пространственной структурой как в свободном пространстве (рис.3), так и на поверхности экрана (рис.4) не сдувается потоком со скоростью до 1 км/с, а лишь слегка размывается. При постановке задачи расчета поджигания топливной смеси присоединенным глубокоподкритическим разрядом, когда стримеры не выходят за пределы плазменного облака, возникает вопрос - а не сдует ли такой разряд скоростным потоком.
а) б) в)
Рис.2. Схема размещения инициатора СВЧ разряда в проточной камере сгорания над металлическим экраном (а), на диэлектрическом экране (б), под диэлектрическим экраном (в).2
А Ж
У
'1Ж ■ ш
уч
а) б)
Рис. 3. Развитый стримерный разряд, зажженный на инициаторе в неподвижном воздухе (а) и в сверхзвуковом потоке воздуха с М=2 при давлении р=133 кПа, при напряженности исходного электрического поля Ео=1,7 кВ/см, Х=8.9 см.3
а) б)
Рис.4. Стримерный разряд на поверхности диэлектрика в неподвижном воздухе (а) и в потоке с М=2, Ео=3,7 кВ/см, Х=8.9 см.4 Однако СВЧ поле - периодическое и задача сводится к уравнению Гельмгольца в форме интегрального уравнения Поклингтона [25].
Наличие экрана усложняет задачу расчета. В общем случае, для моделирования резонансных процессов нужно использовать уравнения газовой динамики и уравнения Максвелла для термодинамически неравновесной плазмы в трехмерной постановке.
2,3,4ЛррепШх 1
4л г
Б(г) = Б0 (г) + i — Га'Б(г')F(сa')...
п
О ( г, г', а') ёг',
О ( г, г', а ) =
F (с) =
ехр (/'кЯ )
к3 Я3 ка
(1 - т)
2 - 3—
\
Я2
+к2 а2
■ ^
2 1 + i
4л ■с
1 + i
4л ■с
ю
Л
(4)
ю
Я = ,1 а2 +(г - г ')2.
В уравнении (4) Б (г) - комплексная амплитуда электрического поля на оси плазменного канала ( г = 0 ), а - электрическая проводимость, ю - круговая частота, к = 2л / Я - волновое число, F - функция, учитывающая скин - эффект, Зх - функция Бесселя.
Как указано в работе [26] поглощение энергии и разогрев окружающий среды производится только стримером, а не всем диффузным плазмоидом, следовательно, для тонкого стримера уравнения могут быть заменены на обыкновенные дифференциальные уравнения. И, при отсутствии сносящегося потока, задача решения уравнений плазмоди-намики сводится к одномерной.
Ранее авторами в такой постановке на основе методики, изложенной в работе [27], были получены решения для радиуса стримера, скорости его распространения, плотности поглощаемой энергии излучения. В трехмерном случае и при наличии сносящегося потока приходится решать полную систему уравнений.
Как известно, плазма разряда в общем случае термодинамически не является равновесной средой. Энергия электромагнитного поля воспринимается прежде всего электронами и затем через упругие и неупругие столкновения передается тяжелой составляющей.
Поэтому необходимо рассмотреть двух-температурную модель плазмы, различающую температуру электронов и температуру тяжелого компонента.
Электронная температура определяет многие константы реакций. Она определяет константы в уравнении баланса ионизации, частоты ионизации, присоединения, рекомбинации и так далее. Электронная температура
определяет время релаксации в энергообмене между электронными и тяжелыми компонентами, потери на рекомбинацию и тормозное излучение, активное и реактивное сопротивление плазмы.
Диссоциация молекулярного газа зависит, в основном, от температуры тяжелых компонентов и в разработанной модели учитывается через уравнение ионизации Саха.
Теплоемкость и теплопроводность тяжелых компонентов определяются в предположении их термодинамического равновесия, характеризующегося температурой газа. В этом же предположении определяются теплопроводность, коэффициент диффузии и теплоемкость электронных компонентов.
После прохождения реакции горения число молекул практически не изменяется, что позволяет использовать уравнение состояния воздуха для описания процессов в воз-душно-пропановой смеси.
Полная система уравнений выглядит следующим образом
Уравнение баланса ионизации:
/ ч ёТ
с„ (ш,р,Т,Т )—^ = -Т ^■У +
е> г) л 8
(^ р,Те ,Тг )уте)+
т - г и
. (5)
_г_г__|--
(^ Р,Те Т ) Ти (р,Т8 )
В (5) ш - коэффициент ионизации; у1. - сумма частот ионизации, присоединения электронов и рекомбинации; п - плотность электронного газа (концентрация электронов); Т - температура электронов; Т^ - температура газа; Д -коэффициент диффузии электронов.
Уравнение баланса мощности (энергии) для электронного газа:
В (6) сте - удельная теплоемкость электронного газа при постоянном объеме, и - скорость потока; т - время химической релаксации электронного газа; Б - амплитуда электрического поля.
сТ &
г - Т
^ — = -Т V-и
Уе А е
—--:--+
(ш,п,Т ,Т )
: V ' е' 8)
+ V■( бе (Ш, П,Те )VTe ) +
(ш П,Те , Тг )■ |БГ
т
п
Уравнение баланса мощности тепловой Q энергии для газа:
ёТ _
ёг
= -Т V-V +
(у, р,Те ,Т?)
V•( Qg (у Р,Те Т )VTe ) + Т- Т + Н (^Р,Т.,Тв ) ТН (Р,Тв )
(7)
Т* (У Р,
В (7) ^ - удельная теплоемкость газа при постоянном объеме, р - плотность газа, Н -внутренняя энергия.
Уравнение движения для газа:
ёи ёг
- ^(рг + упТ).
(8)
Уравнение неразрывности для газа:
ёр
—!- = -рV• и .
ёг
(9)
В воздухе диссоциация кислорода начинается при Т = 2000 К, а азота - при Т = 4000 К. Образование свободных электронов в процессе ионизации имеет место при Т > 10000 -12000 К.
В рассматриваемом случае продолжительность разряда т полагалась примерно равной или чуть большей времени индукции модельной топливной смеси, что позволяло рассматривать процесс «разряд - индукция - горение» как квазинепрерывный.
Поскольку к началу воспламенения высокотемпературные процессы в разряде уже заканчиваются, а горение пропан-воздушной смеси происходит при температурах порядка 1200К, то все перечисленные выше эффекты в расчетах можно не учитывать.
Таким образом, задача постановки начальных условий для каждого периода горения состоит в расчете разряда, температуры электронов Т и газа Т в сносящий потоке, а
также вычисления подведенного за время т количества тепла д . За время т стример поглощает излучение [25]
Уравнение релаксации внутренней химической энергии для газа:
ёН ёг
Н
В (10) т - время релаксации химической энергии.
Уравнение Максвелла для электромагнитного поля:
В (11) - Н - амплитуда магнитного поля, с -скорость света. Решая систему (5-11) численно, можно рассчитать параметры разряда, температуру электронов Т и газа Т , причем
с учетом сноса плазмы спутным потоком газа.
Е 2
Qs = с • Е^ • Sв-т 4 •л
(12)
(10) где с - скорость света, Е0 - эффективная напряженность поля, ^ - эффективное сечение поглощения.
Тогда задача сводится к нахождению среднего за время импульса в сносящий потоке. В общем случае, для этого необходимо решить систему (5-11) на промежутке времени [0, т].
Я =
11.1.5
30 • р
(13)
. 4лс .ю
V х Н =--7 — • Е,
V х Е = -— Н с
Рассмотрим модель подвода тепла и воспламенения топливной смеси. При высокой температуре возбуждаются колебательные степени свободы молекул, могут происходить диссоциация и ионизация.
д = 2.2 •10-
3
X2
р • Е
30 • р
(14)
8
с
р
В формулах (13-14) длина волны 1 - в см, давление р в торр, Е - напряженность поля в В/см, с - теплоемкость при постоянном давлении. В покоящейся среде плазмоид имеет сферическую форму с радиусом, равным радиусу Дебая Яе, откуда следуют, удобные для постановки начальных условий формулы:
Канал подсоединялся к вакуумной камере.
Модель мгновенного подвода тепла [27] в данном случае неприменима, т.к. процесс во
времени сопоставим с длительностью индукции топливной смеси, а параметры температуры (электронной и температуры газа) и количество подводимого в единицу времени тепла q изменяются во времени.
Расчеты сравнивались с результатами экспериментов, которые проводились на установке с 1=12,3 см, схема которой приведена на рис.5.
(9x9) с
Пропан (10x70) мм
ргоря
Рис.5. Схема постановки эксперимента с аэродинамической трубой вакуумного типа.5
II. РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ
Расчеты выполнялись в следующей последовательности. Сначала для каждого из типов экрана находилась резонансная длина вибратора при его размещении вблизи поверхности экрана, на его поверхности, под ней и в свободном пространстве. Вычислялись коэффициенты усиления поля КЕ . Затем для резонансных длин вибратора производился расчет развития разряда и разогрева плазмы в потоке газа. По характеру роста температуры и поглощения разрядом энергии определялась модель подвода энергии к топливной смеси. После этого выполнялся расчет воспламенения топливной смеси при различном соотношении горючего и топлива.
Выполненные расчеты резонансной длины вибратора выявили следующее. В свободном пространстве вибратор пробивал воздух при давлении примерно 370-380 торр.
При помещении вибратора над металлическим экраном в пучность поля на расстоянии Н = 1/4 от экрана пробой происходил при атмосферном давлении. С приближением вибратора к металлическому экрану на расстояние меньше чем /к<1/4 коэффициент усиления К заметно увеличивается по сравнению
с вибратором в свободном пространстве. При этом резонансная длина вибратора растет.
При размещении вибратора непосредственно на поверхности диэлектрического экрана из плексигласа толщиной 4.25 мм поле в зазоре между концом резонатора и экраном в два с лишним раза выше, чем на концах вибратора, размещенного у экрана или вдали от него (рис.6).
Рис.6. Зависимость коэффициента усиления от длины вибратора L, Х=12.3 см, на концах вибратора в свободном пространстве (1), на концах вибратора (2) и в зазоре между вибратором и экраном (3) при размещении вибратора на поверхности экрана из плексигласа.6
Постепенное удаление вибратора от экрана приводит к уменьшению поля в зазоре между ним и экраном. Любопытно, что
' ЛррепШх 1
72
уменьшение расстояния между диэлектрическим экраном и вибратором приводит к уменьшению резонансной длины, а не к её увеличению, как в случае с металлическим экраном.
Так резонансная длина в свободном пространстве примерно 54.2 мм, а на экране из плексигласа только 45.2 мм, из керамики -34.5 мм. Для керамики, измеренные на рабочей частоте / диэлектрическая проницаемость е=10.3 и тангенс угла диэлектрических потерь ¿>=0.015, для плексигласа: е = 2.5 и tgS = 0.01.
Размещение вибратора под экраном приводит к ослаблению поля и увеличению ^ . Максимальное значение поля в зазоре линейно увеличивается с уменьшением частоты от 257 В/м (2.8 ГГц) до 490 В/м (2.1 ГГц), см. рис.6. При частоте использовавшегося в экспериментах генератора 2.45 ГГц Е > 350 В/м.
Пробой наступает при атмосферном давлении, усиление поля слабо зависит от того, над экраном помещен вибратор или под ним. Таким образом, усиление поля вибратором, помещенным на диэлектрический экран примерно эквивалентно вибратору, закрепленному над металлическим экраном на расстоянии к= Я / 4.
2,4 2,6
Частота, ГГц Рис.7. Зависимость напряженности Е поля в зазоре между вибратором и экраном из плексигласа от частоты СВЧ излучения при различных длинах вибратора Ь от 19 мм до 24.4.7
Для контроля численной методики были выполнены расчеты для условий, приведенных в [23,25,26], скорость потока 1 км/с, Я = 8.9 см. Мощность генератора составляла 1
МВт. Расчеты выполнялись для Е0 = 2 - 7 кВ/см. Инициатор располагался сначала в свободном пространстве, а затем на расстоянии к = Я/4 над металлическим экраном. Коэффициент усиления КЕ изменялся от 6 до 25 путем изменения длины резонатора.
Получены следующие результаты. В первые 2-3 мкс резко растет энергоподвод к разряду и температура перегретого плазменного канала. Величина температуры зависит от типа реализуемого разряда. При примерно одинаковом значении Е на резонаторе, более пространственно развитую структуру имеет разряд, реализуемый в поле с большим Е0 и меньшим К . Добавление металлического экрана приводит к удвоению пика температуры, которая может достигать 9000^.
Уже через несколько мкс поглощение разрядом мощности СВЧ излучения р начинает падать и через 30-40 мкс устанавливается асимптотическое значение, которое не зависит от типа разряда (см. рис.7). Температура газа в перегретом канале разряда остается примерно на постоянном уровне, который составляет 3500-3700А". Размеры плазменного облака (продольные и поперечные) примерно пропорциональны Е .
Рис.8. Зависимость поглощаемой разрядом мощности СВЧ излучения от времени для стримерного разряда с развитой пространственной структурой в начальный момент времени (1), для глубоко подкритического разряда (3), для разряда переходного между (1) и (3) типа (2).8
Поскольку время индукции топливной смеси составляет порядок десятков мс, то процесс подвода энергии можно считать квазистационарным.
Из рис.7 следует, что развитая пространственная стримерная структура играет роль
7,8Appendix 1
73
первые 5-10 мкс, далее нагрев газа от первоначальной структуры разряда не зависит.
Таким образом, имеет смысл использовать как можно меньшую величину Е и как можно больший коэффициент усиления КЕ . Предложенная в [26] простая оценка значения радиуса горячего следа Я в квазистационарном режиме для присоединенного разряда
Я =
Р
лмД
у
(15)
является вполне адекватной. В формуле (15) у - эффективный показатель адиабаты, -
скорость набегающего потока, Н - полная энтальпия, Т^0 - начальная температура торможения. Из формулы следует очевидный вывод - при постоянной температуре в перегретом канале, для сохранения постоянной площади поджигания с увеличением скорости потока нужно пропорционально увеличивать подводимую мощность.
Так, при порядке радиуса Я = 2 - 3 мм и типичной температуре Т = 3000К увеличение
скорости топливной смеси с единиц м/с до 1 км/с должно сопровождаться увеличением мощности, поглощаемой разрядом с 3-4 Вт до 2-3 кВт. Иначе разряд сдувается потоком.
Далее были выполнены расчеты для генератора с 1=12,3 см при мощности СВЧ излучения Р=1 -1.5 кВт.
Расчеты выполнялись для вибратора резонансной длины = 60 мм , размещенного на расстоянии Н = 1/4 от алюминиевого экрана, для для вибратора резонансной длины = 48 мм, размещенного на экране из плексигласа и с = 34,7 мм на экране из корунда. Скорость потока изменялась от 20 м/с до 500 м/с (20, 50, 100, 200, 500).
Результаты расчетов сравнивались с экспериментами, которые проводились на установке, схема которой приведена на рис.5. На рис.9 приведены фотографии разрядов над алюминиевым экраном при различной скорости потока.
Область разряда была накрыта сверху прозрачным листом плексигласа, как показано на рис.5. Порядок длины разряда: около 1 см при скорости 1 м/с и порядка 3 м/с при 500 м/с.
б)
Рис.9. Фото разряда в потоке воздуха со скоростью 1,5 м/с (а) и 500 м/с (б), расположенном на расстоянии h = Х/4 от алюминиевого экрана.9
Затем были выполнены эксперименты с инициаторами установленным под листом корунда. Скорость потока также была до 500 м/с. В этих экспериментах и последующих расчетах обеспечивалось значение напряженности поля на полюсах вибратора Е0 = 350 В/см. Измеренные в ходе эксперимента давления пробоя в обоих случаях подтвердили, что выбранный способ размещения вибраторов обеспечивает коэффициент усиления К около 150 (150 в расчетах и 148 по результатам эксперимента). Таким образом, при мощности источника СВЧ излучения 1 кВт численными расчетами и экспериментом показано, что зажигаемые на площади 9 х 9 см глубоко подкритические разряды не сдуваются потоком воздуха вплоть до скоростей 500 м/с.
Расчетные размеры разрядов: диаметр от 1-1.5 мм при скорости 1.5 м/с до 2.5-5 мм при скорости 500 м/с. Температура разряда - более 3000К. Поглощаемая разрядом мощность растет с увеличением скорости и, соответственно, размеров разрядов, в пределах от примерно 15 до 60 Вт.
В завершении были выполнены эксперименты по поджогу пропановой смеси серией глубоко подкритических разрядов размещенных под листом плексигласа над алюминиевым экраном (рис.10-а) и под корундовым экраном (рис.10-б). Разряд уверенно воспламенял поток при скорости до 500 м/с, но горение происходило только в непосредственно в ядре горячего следа разряда. Сам разряд имел характерный поперечный размер 5 мм. При небольших скоростях за разрядом образовывалось пламя, типичное для газовой горелки (рис.11). Из результатов расчетов разрядов в потоке получались размеры плазменного образования, распределение температуры и плотности подвода энергии.
9ЛррепШх 1
б)
Рис.10. Фото поджигания стехиометрической смеси пропана с воздухом глубоко подкритиче-ским присоединенным разрядом в потоке со скоростью 500 м/с, при размещении инициатора на расстоянии / длины волны от алюминиевого экрана (а) и непосредственно на поверхности экрана из корунда (б), на фото поток течет справа налево.10 Далее полагалось, что эти величины остаются постоянными, и они использовались для задания начальных условий для горения. Наличие в потоке электронов, влечение УФ излучения не учитывалось.
Рис.11. Фото поджигания стехиометрической смеси пропана с воздухом глубоко подкритиче-
ским присоединенным разрядом в потоке со скоростью 30 м/с, при размещении инициатора на расстоянии / длины волны от алюминиевого экрана, на фото поток течет слева.11
При умеренной скорости потока (20 м/с) и стехиометрическом соотношении топлива и воздуха получено очень хорошее совпадение результатов расчетов и эксперимента (см. рис.12).
Анализ результатов расчетов и их сравнение с результатами экспериментов позволил сделать следующие выводы.
Картина течения, образующаяся при поджигании топливной смеси присоединенным глубоко подкритическим разрядом, подобна идеальной горелке.
Моделирование разряда возможно заменой бесконечно тонким диском, к которой подводится энергия.
При нормальном соотношении топлива и окислителя, не слишком сильно отличающегося от стехиометрического, чисто термический механизм воспламенения хорошо объясняет получающуюся картину течения. Температура внутри разряда настолько велика, что внутри него горит даже очень бедная смесь с коэффициентом избытка пропана 0.2 (см. рис.13). При выходе из разряда горение прекращается. При бедных топливных смесях результаты эксперимента и расчетов существенно отличаются, следовательно, наличие УФ излучения и свободных электронов оказывают существенное влияние на процессы химического окисления топлива
а)
б)
г)
Рис.12. Сравнение результатов экспериментов горения разряда в воздухе (а) и в стехиометрической смеси пропана с воздухом (б) с результатами расчета поля температуры (в), поля концентрации СО2 (г). Скорость потока 20 м/с. Средний радиус разряда 2 мм, мощность тепловыделения разряда около 30 Вт.12
io,n,i2Appendix 1
75
Разряд (Discharge)
СО2
Рис.13. Концентрация продуктов горения в следе разряда, коэффициент избытка пропана 0.2, скорость 20 м/с.13
При концентрации топлива ниже «бедного» концентрационного предела мощности, подводимой от разряда к топливной смеси, не хватает для её поджигания, если не учитывать наличие плазмы и УФ излучения, т.е. чисто термический механизм не позволяет объяснять экспериментально наблюдаемое горение экстремально бедных смесей.
Чтобы разряд не сдувался и горение происходило устойчиво, необходимо увеличивать мощность, подводимую к резонатору пропорционально увеличению скорости (см. формулу 15).
III. ЗАКЛЮЧЕНИЕ
Исследования выявили ряд важных фактов. Наиболее существенный научный результат состоит в том, что глубокоподкрити-ческий разряд, присоединенный к инициатору, расположенному на диэлектрическом экране, не сдувается скоростным потоком и уверенно поджигает смесь. Разряд, зажженный на резонаторе вблизи диэлектрического или металлического экрана не сдувается потоком со скоростью в пределах от 50 м/с до 500 м/с. Разряд при этих скоростях уверенно поджигает как богатую, так и бедную топливную смесь. При дозвуковых скоростях факел распространяется в строгом соответствии с законами распространения ламинарного фронта пламени. При сверхзвуковых скоростях горение происходит в пределах конуса Маха, опирающегося на область разряда.
Вблизи «бедного» концентрационного предела горения поджигается и горит только топливо, непосредственно прошедшее через разряд. Результаты расчетов подтверждены экспериментом. Характер горения зависит от плотности энергии в области разряда.
Работа выполнена при финансовой' поддержке Министерства образования и науки РФ, Соглашение № 14.577.21.0277 (ЭБ 07515-2019-9 15), уникальный идентификатор
прикладных научных исследовании -RFMEFI57717X0277.
APPENDIX 1 (ПРИЛОЖЕНИЕ 1)
1Fig.1. Streamer microwave discharge in free space (a), on the surface of textolite (b) and glass (c), on a nylon mesh (d).
2Fig.2. Scheme of the initiator of the microwave discharge in the flowing combustion chamber above the metal screen (a), on the dielectric screen (b), under the dielectric screen (c).
3Fig.3. Spatial streamer discharge ignited on the initiator in still air (a) and in a supersonic air stream with M = 2 at a pressure p = 133 kPa, with an initial electric field strength of E = 1.7 kV / cm, 2 = 8.9 cm.
4Fig.4. Streamer discharge on the dielectric surface in still air (a) and in a stream with M = 2, E = 3.7 kV / cm, 2 = 8.9 cm.
5Fig.5. Experimental setup.
6Fig.6. The function of the gain on the length of the vibrator L, 2 = 12.3 cm, at the ends of the vibrator in free space (1), at the ends of the vibrator (2) and in the gap between the vibrator and the screen (3) when placing the vibrator on the surface of the screen made of plexiglass.
7Fig.7. Relation of the field strength E in the gap between the vibrator and the plexiglas screen on the microwave radiation frequency for various vibrator lengths L from 19 to 24.4 mm.
8Fig.8. Relation of the microwave radiation power absorbed by the discharge on time for a streamer discharge with a developed spatial structure at the initial instant of time (1), for a deeply subcritical discharge (3), for a transitional discharge between (1) and (3) type (2).
9Fig.9. Photo of a discharge in an air stream at a speed of 1.5 m / s (a) and 500 m / s (b) located at a distance of h = 2/4 from the aluminum screen.
10Fig.10. Photo of burning a stoichiometric propane-air mixture with a deeply subcritical attached discharge at a speed of 500 m / s, when the initiator is placed at a distance % of the microwave radiation wavelength from the aluminum screen (a) and directly on the surface of the corundum screen (b), in the photo the mixture flows from right to left.
13 Appendix 1
76
11Fig.11. Photo of burning a stoichiometric propane-air mixture with a deeply subcritical attached discharge in a stream at a speed of 30 m / s, when the initiator is placed at a distance % of the wavelength from the aluminum screen, the photo flows from left to right.
12Fig.12. Comparison of experimental results and calculations of discharge burning in air (a), in a stoichiometric mixture of propane with air (b) with the results of calculation of the temperature field (c), the concentration field of CO2 (d). The flow velocity is 20 m / s. The average radius of the discharge is 2 mm, the heat dissipation power of the discharge is about 30 watts.
12Fig.13. The concentration of combustion products in the Wake of discharge, coefficient of excess propane to 0.2, the speed of 20 m/s.
ЛИТЕРАТУРА (REFERENCES)
1. Chernyshev S.L., Skvortsov V.V., Ivanov V.V., Troschinenko G.A. A concept for generation and application of body- centered non-equilibrium discharge for initiation and intensification of fuel combustion in high-speed flows. Aviatsionnaya Promyshlennost', 2013, no. 2, pp. 19-25.
2. Rapp V.H., DeFilippo A., Saxena S., Chen J.Y., Dibble R.W., Nishiyama A., Moon A., Ikeda Y. Extending Lean Operating Limit and Reducing Emissions of Methane Spark-Ignited Engines Using a Microwave-Assisted Spark Plug. Journal of Combustion, Volume 2012, no 5, pp.1-8. Article ID 927081. doi:1 0.1155/2012/927081.
3. Starikovskiy A., Aleksandrov N. Plasma-assisted ignition and combustion. Progress in Energy and Combustion Science, 2013, vol. 39., pp. 331-368.
4. Ikeda Y., Nishiyama A., Wachi Y., Kaneko M. Research and development of a microwave plasma combustion engine (Part I: Concept of plasma combustion and plasma generation technique), SAE Technical Paper, 2009. 8 p. doi: https://doi.org/10.4271/2009-01-1050.
5. Ikeda Y., Hiroki K., Jeonj M., Kaneko H. Research and development of microwave plasma combustion engine (Part II: Engine performance of plasma combustion engine), SAE Technical Paper, 2009. doi: https://doi.org/10.4271/2009-01-1050.
6. Khodataev K.V. The Ignition of the Combustion and Detonation by the Undercritical Microwave Discharge. AIAA Plasmadynamics and Laser Conference, 32nd, and Weakly Ionized Gases Workshop, 4th, Anaheim, CA. Jun 11-14, 2001. AIAA Paper 2001-2941.
7. Khodataev K.V. The power effectivity of a microwave undercritical attached discharge, initiated
by resonant vibrator. 43-rd AIAA Aerospace Sciences Meeting and Exhibit. 10-13 January 2005, Reno, Nevada. AIAA-2005-0596. p.6.
8. Khosravyel M. H. Review of the New Combustion Technologies in Modern Gas Turbines, Progress in Gas Turbine Performance, 2013, 268 p., ISBN 978-953-51-1166-5. doi: 10.5772/2797.
9. Nishiyama A., Ikeda Y. Improvement of Lean Limit and Fuel Consumption Using Microwave Plasma Ignition Technology, SAE Technical Paper, 2012. doi: https://doi.org/10.4271/2012-01-1139.
10. Bulat, M.P., Bulat, P.V. The analysis centric isen-tropic compression waves. World Applied Sciences Journal, 2013, 27(8), c. 1023-1026.
11. Uskov, V.N., Bulat, P.V., Arkhipova, L.P. Classification of gas-dynamic discontinuities and their interference problems. Research Journal of Applied Sciences, Engineering and Technology, 2014, 8(22), c. 2248-2254.
12. Viktorovich, B.P. About the detonation engine. American Journal of Applied Sciences, 2014. 11(8), c. 1357-1364.
13. Uskov, V.N., Bulat, P.V. Shock and detonation wave in terms of view of the theory of interaction gasdynamic discontinuities. Life Science Journal, 2014. 11(SPEC. ISSUE 8), c. 307-310.
14. Bulat P.V., Denissenko P.V., Volkov K.N. Trends in the development of detonation engines for high-speed aerospace aircrafts and the problem of triple configurations of shock waves. Part I. Research of detonation engines. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2016, vol. 16, no. 1, pp. 1-21.
15. Bulat P.V., Denissenko P.V. Trends in the development of detonation engines for high-speed aerospace aircrafts and the problem of triple configurations of shock waves. Part II - Research of counterpropagating shock waves and triple shock wave configurations. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2016, vol. 16, no. 2, pp. 199-223. doi:10.17586/2226-1494-2016-16-2-199-223.
16. Bulat, P.V., Chernyshev, M.V. Existence regions of shock wave triple configurations. International Journal of Environmental and Science Education, 2016, 11(11), c. 4844-4854.
17. Kirill V. Khodataev. Weak detonation wave ignition and sustaining in over CJ-speed flow by means of undercritical microwave discharge. Symposium on «Thermo-chemical and plasma processes in aerodynamics». Sankt-Petersburg. 19-21 June, 2006.
18. Westbrook, C.K., Dryer, F.L. Simplified Reaction Mechanisms for the Oxidation of Fuels in Flames. Combustion Science and Technology, 1981, Vol.27, pp. 31-43.
19. Katta, V.R., Roquemore, W.M. Simulation of Unsteady Flows in an Axisymmetric Research Com-bustor Using Detailed-Chemical Kinetics.
AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, 34th, Cleveland, OH, July 1315, 199.
20. V. A. Levin, G. D. Smekhov, A. N. Khmelevsky. Simulation of equilibrium combustion of kerosene by a mixture of gaseous hydrocarbons. Physics of combustion and explosion, 2007, vol. 43, No 4, p. 3-7. 2.
21. Bulat P.V., Upyrev V.V. Calculation of the shock-wave structures in an imperfect gas // Journal of Engineering and Applied Sciences - 2016, Vol. 11, No. 13, pp.
22. Khodataev K. V. The Nature of Surface MW Discharges //48th AIAA Aerospace Sciences Meeting and Exhibition. 4-8 January 2010, Orlando, Florida. Paper AIAA 2010-1378.
23. Khodataev K.V. The initial phase of initiated un-dercritical microwave discharge. 43-rd AIAA
Aerospace Sciences Meeting and Exhibit. 10-13 January 2005, Reno, Nevada. AIAA-2005-0598.
24. Khodataev K.V., GorelikB.R. Diffusion and drift regimes of a plane ionization wave propagation in UHF field. //Physika plasmy, 1997, v.23, 3 pp. 236-245.
25. Pocklington H.C. Camb. Phil. Soc. Proc., 1897, v.9, p.324.
26. Khodataev K. V. The physical basis of the high ability of the streamer MW discharge to a resonant absorption of MW radiation. 42nd AIAA Aerospace Sciences Meeting 5-8 January 2004, Reno, NV. Paper AIAA-2004-0180.
27. Volkov K. N., Emelyanov V. N. Gas flows with mass flow in channels and paths of power plants. M: FIZMATLIT, 2011. - 455 p.
Сведения об авторах ft.,*
Есаков Игорь Иванович. Доктор физико-математических наук, заместитель генерального директора АО «МРТИ РАН» по научной работе. Область научных интересов: физика плазмы, СВЧ разряды.
E-mail: esa-kov@mrtiran.ru
Грачев Лев Петрович.
Руководитель экспериментальных исследований в НТЦ «Плазменных технологий» АО «МРТИ РАН». Область научных интересов: электрические разряды в газах E-mail: lqrachev@mrtiran.ru
Волобуев Игорь Алексеевич.
Университет ИТМО.
Область научных
интересов:
компрессорная
техника.
Е-таП:уо1о-
Ьие\ ig@mai1.ru
Раваев Александр Александрович.
Доктор физико-математических наук. Область научных интересов: физика плазмы, СВЧ разряды. E-mail:
aravaev@mrtiran.ru