авторефераты диссертаций БЕСПЛАТНАЯ  БИБЛИОТЕКА

АВТОРЕФЕРАТЫ КАНДИДАТСКИХ, ДОКТОРСКИХ ДИССЕРТАЦИЙ

<< ГЛАВНАЯ
АГРОИНЖЕНЕРИЯ
АСТРОНОМИЯ
БЕЗОПАСНОСТЬ
БИОЛОГИЯ
ЗЕМЛЯ
ИНФОРМАТИКА
ИСКУССТВОВЕДЕНИЕ
ИСТОРИЯ
КУЛЬТУРОЛОГИЯ
МАШИНОСТРОЕНИЕ
МЕДИЦИНА
МЕТАЛЛУРГИЯ
МЕХАНИКА
ПЕДАГОГИКА
ПОЛИТИКА
ПРИБОРОСТРОЕНИЕ
ПРОДОВОЛЬСТВИЕ
ПСИХОЛОГИЯ
РАДИОТЕХНИКА
СЕЛЬСКОЕ ХОЗЯЙСТВО
СОЦИОЛОГИЯ
СТРОИТЕЛЬСТВО
ТЕХНИЧЕСКИЕ НАУКИ
ТРАНСПОРТ
ФАРМАЦЕВТИКА
ФИЗИКА
ФИЗИОЛОГИЯ
ФИЛОЛОГИЯ
ФИЛОСОФИЯ
ХИМИЯ
ЭКОНОМИКА
ЭЛЕКТРОТЕХНИКА
ЭНЕРГЕТИКА
ЮРИСПРУДЕНЦИЯ
ЯЗЫКОЗНАНИЕ
РАЗНОЕ
КОНТАКТЫ

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

На правах рукописи

Каменщиков Сергей Александрович

ВЗАИМОДЕЙСТВИЕ ПЛАЗМЫ ПРОДОЛЬНО – ПОПЕРЕЧНОГО И

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

ВОЗДУШНО - ПРОПАНОВЫМ ПОТОКОМ

01.04.08 – физика плазмы.

Автореферат диссертации на соискание ученой степени

кандидата физико-математических наук

.

Москва - 2011.

Работа выполнена на кафедре физической электроники физического факультета Московского Государственного Университета им. М.В. Ломоносова.

Научный руководитель:

Кандидат физико-математических наук, доцент Черников В.А.

Официальные оппоненты:

Доктор физико-математических наук, профессор Коссый И.А.

Кандидат физико-математических наук Артамонов В.И.

Ведущая организация:

ФГУП «Московский Радиотехнический Институт» Российской Академии Наук.

Защита состоится 12 мая 2011 г. в 16:30 на заседании диссертационного совета Д 501.001.66. при Московском государственном университете им. М.В.Ломоносова, 119991, ГСП-1, Москва, Ленинские горы, МГУ им. М.В.Ломоносова, Дом 1, строение 2, Физический Факультет, Северная физическая аудитория.

С диссертацией можно ознакомиться в библиотеке физического факультета МГУ им. М.В.Ломоносова.

Автореферат разослан 11 апреля 2011 г.

Ученый секретарь диссертационного совета Д 501.001.66.

кандидат физико-математических наук Карташов И.Н.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Диссертация посвящена определению нестационарных свойств продольно – поперечного и плазмодинамического разрядов в условиях натекающего воздушно – пропанового потока. Исследовано влияние свойств разрядной плазмы на предварительное перемешивание смеси и характеристики плазменного воспламенения.

Актуальность темы Создание летательных аппаратов (ЛА), для которых натекающий поток воздуха является надкритическим относительно гидродинамических параметров подобия (число Маха, число Рейнольдса, число Струхаля и т.д.) связано с рядом проблем прикладного и фундаментального характера. В частности, необходимость новых подходов возникает в случае движения ЛА со сверхзвуковыми (число Маха M 1) и гиперзвуковыми (M 5) скоростями. Задачи, возникающие при этом, могут быть решены в рамках применения плазменных источников различных типов. Весь сложный комплекс подобных проблем можно условно разделить на три основных группы:

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

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

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

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

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

Цель работы В качестве основных целей работы необходимо выделить следующие:

Определение нестационарных свойств разряда в присутствии воздушно – пропанового и воздушного течения.

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

Достижение оптимального перемешивания натекающей топливно – воздушной смеси в условиях исследованных режимов существования разряда.

Достижение плазменного воспламенения в условиях натекающего воздушно – пропанового надкритического потока. Определение влияния нестационарных свойств разряда на инициирование плазменного воспламенения.

Методы исследований Результаты диссертации получены с использованием следующего диагностического оборудования:

Оптические детекторы – фотоэлектронный умножитель, рефракционный датчик, эмиссионный спектрометр.

Электромеханические детекторы – датчики давления тензометрического типа, трубки Пито – Прандтля, термопарный датчик температуры.

Методы определения электрических параметров разрядов: делитель напряжения, шунт.

Научная новизна Все полученные результаты являются новыми, основная их часть получена автором самостоятельно. Научная новизна работы состоит в следующем:

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

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



Выяснено влияние неустойчивого продольно – поперечного разряда на усиление хаотических свойств течения за областью разряда. Установлена связь между средними характеристиками разряда и спектральными свойствами течения.

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

Реализовано плазменное воспламенение потока при помощи неустойчивого продольно – поперечного и плазмодинамического разрядов. Показано влияние средних характеристик разряда на оптимальность воспламенения.

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

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

Апробация работы Результаты диссертации докладывались:

На «Звенигородской конференции по физике плазмы и управляемому термоядерному синтезу», Звенигород, 2007,2008 гг.;

на конференции «Ломоносовские чтения» в МГУ, г. Москва, 2007,2008 гг.;

на семинаре «Школа – семинар по магнитоплазменной аэродинамике», ИВТАН, г. Москва, 2007-2009 гг.;

на конференции «AIAA Aerospace Sciences Meeting», США, 2009,2010 гг.;

на конференции «6th International Workshop and Exhibition on Plasma Assisted Combustion, Germany», Германия, 2010 г.;

на семинарах кафедры физической электроники физического факультета МГУ, 2008-2009 гг.

Публикации.

Основные результаты по теме диссертации опубликованы в семи работах, список которых приводится в конце автореферата. Шесть работ [2-7] опубликованы в изданиях, входящих в утвержденный ВАК перечень ведущих рецензируемых научных журналов и изданий, в которых должны быть опубликованы основные научные результаты диссертации на соискание ученой степени доктора и кандидата наук.

Структура и объем диссертации.

Диссертация состоит из введения, пяти глав и заключения. Первая глава носит вводный характер и содержит обзор литературы. Вторая глава содержит описание экспериментальной установки, методов исследования и результаты предварительного исследования течения без присутствия источников плазмы и инжекции пропана. Третья, четвертая и пятая главы оригинальны. В целом диссертация содержит 125 страниц, включая 62 рисунка, 1 таблицу и библиографию из 109 наименований.

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

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

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

Рис. 1. Схема первой экспериментальной установки. 1 – вакуумная камера, 2 – компрессор, 3 – кран, 4 – воздухопроводы, 5 – сборочный коллектор, 6 – баллонная рампа, 7 – манометры, 8 – управляемый пневмоклапан, 9 – омический нагреватель, 10 – блок питания нагревателя, 11 – сверхзвуковой канал, 12 – датчики давления, 13 – соединение с вакуумной системой, 14 – иллюминаторы, 15 – топливораздаточные фланцы, 16 – фланец для крепления плазменных источников воспламенения, 17 – система подачи топлива.

Основным элементом экспериментальной установки является вакуумная камера (1) объемом 3 м3, с которой через гибкое соединение (13) соединен канал (11). Принцип работы установки следующий: с помощью компрессора (2) воздух по трубопроводу (4) через сборочный коллектор (5) нагнетается в дополнительную баллонную систему (6). Ее наличие позволяет довести общий объем системы высокого давления до 0.5 м3.

Рабочее давление в системе высокого давления (системе подачи воздуха) контролируется с помощью манометров (7) и может достигать 6 Бар. Из системы подачи воздуха поток через управляемый пневматический клапан (8) поступает в омический нагреватель (9), и далее в аэродинамический канал (11). В настоящую работу результаты исследований влияния предварительного нагрева не были включены.

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

Формирующая секция предусмотрена с целью создания однородного потока. Для осуществления оптической диагностики в двух секциях вмонтированы кварцевые иллюминаторы. Внутреннее сечение рабочих секций составляет 2,5x4 см2. Внутренний диаметр формирующей секции D=2.5 см. Таким образом, на границе между формирующей секцией и следующей за ней рабочей секцией существует обратный уступ, создающий зону внезапного расширения (застойную зону). Канал оборудован системой подачи топлива (17) и системой контроля давления (12) в виде датчиков давления ( единиц). Во всех экспериментах в качестве топлива использовался пропан. Поступление пропана осуществляется при помощи двух фланцев (15). Первый фланец размещен между секцией сопла и формирующей секцией;





второй фланец – между формирующей и рабочей секциями. В центре обоих фланцев созданы цилиндрические отверстия диаметром 25 мм.

Внутри первого фланца размещены четыре канала ввода пропана, позволяющих осуществлять инжекцию перпендикулярно направлению потока. Диаметр каналов d=2 мм.

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

При проведении исследований, результаты которых представлены в Главе 3 и Главе 4, был использован модифицированный экспериментальный стенд (ЭУ2). В этом случае канал прямоугольного сечения был непосредственно соединен с камерой низкого давления (рис.1, позиция 13).. Электродный узел был размещен во второй секции канала.

Основные параметры представленных экспериментальных стендов следующие:

число Маха сопла Лаваля М = 2, давление системы создания низкого давления рст = 0.1 Бар, давление системы создания высокого давления p = 2.6 Бар;

давление в системе подачи пропана pf =4 Бара, время процесса взаимодействия разряда с воздушно – пропановым потоком 1-2 секунды. При помощи датчиков давления и сопел Вентури, расположенных в системах подачи компонент смеси определялся расход обоих компонент и входной коэффициент эквивалентности смеси, который был равен E=1.1. Для определения локализации сверхзвуковой области течения в канале модифицированной установки было проведено предварительное исследование. При помощи трубки Пито, область измерения которой располагалась на оси симметрии канала, было определено локальное число Маха, усредненное по времени существования разряда. Исследование проводилось без включения источников плазмы и системы подачи пропана. Трубка Пито состояла из двух датчиков статического и полного давления потока, расположенных в непосредственной близости. Отношение полного и статического давлений в области позволяло определить число Маха натекающего потока. М, используя измерения следующее соотношение [1], выражающее отношение полного и статического давлений.

166,7 M pf (1) = ( 7 M 2 1 )2, ps Распределение локального числа Маха в центре симметрии прямоугольного сечения рабочего канала, представлено в Таблице 1. Координата z соответствует расположению среза трубки, предназначенной для измерения полного давления.

Статистическая ошибка метода M=0.1 при времени усреднения показаний T=1 c. Левый столбец Таблицы 1 описывает качественно местоположение трубки Пито – Прандтля. За начало координат принята координата центра выходного сечения сопла.

Таблица. 1. Распределение числа Маха на оси модифицированного канала.

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

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

Продольно – поперечный разряд постоянного тока создавался при помощи электродного узла, размещаемого через нижний фланец второй секции канала (рис.2).

Рис.2. Схема размещения электродов во второй секции канала. 1 – катод, 2 – анод, 3 – люк крепления электродов, 4 – нижняя стенка канала. Стрелка – направление потока.

Координата центра межэлектродного зазора могла быть изменены при помощи фиксации высоты электродов. В настоящем исследовании она соответствовала координате центра сечения канала. Длина горизонтальной части верхнего электрода L1=10 мм, а горизонтальной части нижнего электрода L2=40 мм. Источник, питающий разряд, работал в режиме генератора тока, при этом выходное напряжение источника питания U5 кВ, а выходной ток I20 А. Вкладываемая в разряд средняя мощность не превышала P=3 кВт.

В качестве определяющего стационарного параметра, характеризующего разряд, был выбран средний ток разряда [2].

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

Статическое давление в разряде в неподвижном воздухе было выбрано идентичным давлению в разрядной камере в воздушном потоке. Измерение давления было при этом реализовано при помощи тензометрического датчика, расположенного на середине верхней стенки разрядной секции. Съемка разряда в неподвижном воздухе при помощи видеокамеры показала наличие осцилляций святящегося образования, представленного на рис.3 а. Осцилляции с осью вращения, расположенной в конце анода представляют собой блуждание катодного пятна, что позволяет сделать вывод о частичной эрозии нижнего электрода и его сопутствующем нагреве. Данная конфигурация разряда состоит исключительно из приэлектродных областей. Положительный столб не наблюдается. В эмиссионном спектре области свечения разряда в неподвижном воздухе при этом были отмечены интенсивные линии меди с длиной волны 1=510 нм и 2=578 нм.

Интенсивность данных линий превосходит интенсивность линий, соответствующих излучению иона N2+, H() и радикала CN(0,0). Анализ эмиссионного спектра разрядной области в воздушном потоке показал появление линии атомарного кислорода 1=777,2 нм, интенсивность которой соответствует средней интенсивности линий меди. Таким образом, появление сноса разрядного канала (рис.3 б) приводит к усилению процессов диссоциации и возрастанию роли положительного столба.

а) б) в) Рис. 3. Кадры видеосъемки разряда для случая разряда в неподвижном воздухе (а) с экспозицией 20 мкс, разряда в сверхзвуковом потоке воздуха с экспозицией 20 мкс (б) и в сверхзвуковом потоке воздуха с экспозицией 1 с (в). Стрелка – направление потока.

Кадр (а) соответствует более подробному масштабу сравнительно с б) и в).

Данные скоростной видеосъемки позволили заметить существование нескольких фаз разряда: реализация разряда в области межэлектродного зазора (фаза I, рис.3 (а)) и конвективная фаза (фаза II, рис.3 (б)). Наличие первой фазы связано с локальным торможением потока вблизи края катода в связи с наличием косых скачков и нагрева газа.

Существование разряда в области пониженной скорости приводит к преобладанию скорости дрейфа заряженных частиц над скоростью сноса. В связи с этим разряд сохраняет стабильную фазу. Вторая фаза развития разряда (конвективная фаза), обусловленна нарушением указанного выше условия. Данная фаза обладает отличительной особенностью относительно поперечного разряда постоянного тока в потоке. Присутствие катода приводит к наличию пограничного слоя и как следствие – к созданию циркуляции скорости. Данный эффект создает условия для деформации положительного столба. При этом макроскопическим масштабом деформации является межэлектродное расстояние. Обработка последовательности кадров скоростной видеосъемки позволила определить относительное время существования фаз разряда в воздушном потоке. Был проанализирован квадратный участок изображения, на диагонали которого лежал отрезок, соединяющий заостренные края электродов. Длина отрезка соответствовала минимальному межэлектродному расстоянию. Область анализа была ассоциирована с существованием разряда в первой фазе. При помощи программы, написанной на языке Mat lab 7, изображение было переведено в массив, элементами которого являлась яркость пикселей участка в формате RGB. Пример сигналов интенсивности в окне анализа представлен на рис. 4.

а) б) Рис. 4. Сигнал интенсивности свечения в окне анализа.

Римскими цифрами обозначены фазы разряда. а) Средний разрядный ток I=19 А.

б) Средний разрядный ток I=7 А.

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

При помощи калиброванной термопары ХА-68, позволяющей измерять температуру в пределе от -50 °С до 1300 °С было проведено исследование температуры в следе за разрядом. Для этого термопара было размещена в начале третьей секции канала (дозвуковая область), за разрядом (Таблица 1). Зависимость максимального значения э.д.с.

сигнала термопары в области насыщения от среднего разрядного тока представлена на рис. 5. Тенденция роста температуры в спутном следе газа за разрядом позволяет сделать вывод о том, что определяющую роль в нагреве газа играет ток.

Рис.5. Зависимость максимального потенциала термопары в зависимости от среднего тока разряда.

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

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

Первый путь связан с использованием части положительного столба предшественника и приводит к вторичному пробою в идентичной фазе разряда: фазе – II. Второй путь сопряжен с полной рекомбинацией плазмы плазменного канала – предшественника и пробою в минимальном межэлектродном зазоре, т.е. в фазе – I. Механизм вторичного пробоя, реализованный в данном типе разряда, обусловлен фактором увеличения длины положительного столба [2] и сопутствующим ростом напряжения на разряде. В первом случае рост напряжения может привести к достижению пробойного значения и реализации пробоя в фазе – I или фазе – II (с участием части канала предшественника).

Средняя частота вторичных пробоев при этом определяется отношением средней напряженности электрического поля в разряде и пробойного электрического поля [2].

2 E u (2) f1 min E пр L В приведенном соотношении E – напряженность электрического поля разряда, u – скорость движения заряженных частиц, Eпр – пробойная напряженность электрического поля, являющаяся функцией геометрии электродов, L0 – межэлектродное расстояние. В цитируемой работе используется предположение о малости скорости дрейфа заряженных частиц относительно скорости потока и наличие фиксированного пространственного масштаба (межэлектродного расстояния). Для условий настоящей работы, отмеченная закономерность не может полностью описать среднюю частоту пульсаций, связанных с разрядом ввиду существования разрядных фаз и флуктуаций потока. Заметим, что в случае вторичного пробоя во второй фазе пробойное напряжение отличается от пробойного напряжения в минимальном межэлектродном зазоре и поэтому указанное соотношение для частоты является нижним пределом. Прямое использование соотношения для f1 может быть реализовано для состояний разряда с преобладанием фазы - I (низких токов).

Второй механизм, ответственный за возникновение вторичного пробоя – распад плазмы за счет локализации положительного столба в области более низкого электрического поля, например, вне плоскости электродов. На основании принципа размерности, можно выразить среднюю частоту распада f2 через скорость движения заряженных частиц и пространственный масштаб_единственным образом:

u f 2 const (3) _ L Здесь u – скорость движения заряженных частиц, L – пространственный масштаб. Верхнее подчеркивание соответствует процедуре усреднения по времени и в дальнейшем будет опущено. Запишем стандартные соотношения для плотности тока разряда [3]:

2.82 10 4 n U (4) j q nu E I /s L В соотношении (3) j – плотность тока, q – заряд компоненты плазмы, n – концентрация заряженных частиц, – проводимость плазменного канала, U – напряжение разряда, - эффективная частота столкновений электронов, I – ток, протекающий через плазменный канал, s – сечение плазменного канала. Выразим длину плазменного канала из соотношения (4). Тогда на основании полученной падающей вольт – амперной характеристики линейного типа может быть записано выражение для второй частоты:

( ( p ) e Z ) u ' (5) f 2 const ( k 0 k1 I ) Выразив длину плазменного канала в соотношении (2) можно записать конечное выражение для частоты (6), в котором const и const – безразмерные константы.

Преобладание одного из факторов гибели плазменного канала является ограничивающим в рассматриваемой системе.

( ( p ) e Z ) u 2 e ( p ) u f max const ' (6), const ( k 0 k1 I ) E пр L Как видно из соотношения, в пределе высоких токов, определяющим становится эффект f2. В настоящей работе было распада плазменного канала, определяемый частотой исследовано влияние среднего разрядного тока на максимальную частоту пульсаций. Для этого в начале третьей секции канала была расположена трубка Пито (см. Таблица 1).

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

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

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

Тогда анализ временного хода характеризующей поток величины может позволить определить рост энтропии в определенном диапазоне временных масштабов, или частичную энтропию. При этом соотношение для частичной энтропии может быть выражено следующим образом [4]:

ln(N ()) S p lim ( D ) (7) D lim 0 ln()) Здесь Sp – частичная энтропия, D – обобщенная размерность фазового пространства, N – число элементов с временным масштабом, покрывающих данное пространство.

Процедура, позволяющая по временному ходу заданной величины определить число временных масштабов, была предложена Х.Харстом [5]. Здесь Sp – частичная энтропия, D – обобщенная размерность фазового пространства, N – число элементов с временным масштабом, покрывающих данное пространство. Процедура, позволяющая по временному ходу заданной величины определить число временных масштабов, была предложена Х.Харстом [5] и получила название RS анализа временного сигнала.

Изложим расширенный алгоритм RS анализа Харста на примере временного ряда XN:

1). Составление переменного ряда объемом Nmin из начального временного ряда.

t 2). Составление кумулятивного ряда X t, N u 1 (eu M N ), где eu - член переменного ряда, MN – среднее ряда, N принимает значения от Nmin до максимального значения, t принимает значения от 1 до N-1. 3). Определение размаха переменного ряда R=max(Xt,N) – min(Xt,N). 4). Определение стандартного отклонения ряда R=max(Xt,N) – min(Xt,N).;

5). Построение зависимости ln(R/S) от ln(N). 6). Определение показателя Харста как тангенса угла наклона данной зависимости. 7). Определение обобщенной размерности фазового пространства. 8). Определение частичной энтропии Sp.

RS H, Итогом процедуры анализа Харста является показатель Харста характеризующий изменение статистических свойств ряда XN в зависимости от масштаба времени. При этом обобщенная размерность D связана с показателем Харста соотношением D=2-H.

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

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

Рис. 7. Зависимость приведенного приращения частичной энтропии в следе за разрядом от среднего тока разряда.

За масштаб принята энтропия невозмущенного потока.

В четвертой главе настоящей работы приведено исследование влияния инжекции пропана на нестационарные свойства продольно – поперечного разряда. Показана связь между средними параметрами разряда и перемешиванием воздушно – пропановой натекающей смеси. Влияние добавления пропана в воздушный поток может проявляться в результате следующих процессов: газодинамическое влияние, влияние пропана на скорость элементарных процессов в плазме и его влияние на параметры химических реакций, протекающих вне плазменной среды. Газодинамическое влияние сводится к изменению средних безразмерных характеристик потока и амплитуды его возмущений без присутствия разряда. Инжекция струи пропана в область воздушного потока приводит к возникновению дополнительного возмущения и росту осцилляции газодинамических параметров. При помощи трубки Пито – Прандтля, размещенной в конце второй секции, были определены статическое и полное давление. На основании сигналов, усредненных по времени T=0,2 c были получены стандартные отклонения. При добавлении пропана в поток обнаружен рост стандартного отклонения статического давления в st=1.5 раз и полного давления в f =22,9 для стандартной схемы инжекции, рассмотренной во второй главе. Измерения проводились без включения разряда, что позволило сделать вывод об интенсификации возмущений, связанных непосредственно с инжекцией пропана.

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

Изучение данных скоростной видеосъемки показало, что при инжекции пропана качественно сохраняются все нестационарные свойства разряда в воздушном потоке:

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

Обнаружено, что путь вторичного пробоя не соответствует минимальному пути между частью разрядного канала - предшественника и катодом как в случае воздушного течения.

Вместо этого вторичный пробой реализуется по криволинейной трехмерной траектории (рис. 8). Был выполнен анализ кадров видеосъемки при помощи алгоритма, идентичного тому, что рассмотрен в Главе 3. Было проведено суммирование яркости всех пикселей для каждого кадра из анализируемого окна и определена величина I I min 2.

N Структура этого соотношения соответствует структуре стандартного отклонения, в котором вместо среднего значения интенсивности взято его минимальное значение Imin.

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

Зависимость от среднего тока разряда представлена на рис. 9. Величина представлена в относительных единицах.

Рис. 8. Вторичный пробой. Два последовательных кадра скоростной видеосъемки, выполненные с экспозицией 20 мкс. Вторичный пробой с переходом в фазу - II.

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

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

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

Графическая зависимость, соответствующая данным измерениям представлена на рис. 10.

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

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

Рис. 9. Зависимость меры присутствия разряда в фазе – I от среднего тока разряда.

Рис. 10. Зависимость меры присутствия разряда в фазе – I от среднего тока разряда.

Обработка сигналов статического давления позволила получить ряд значений метрической энтропии от токовых режимов (рис. 11).

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

Рис. 11. Зависимость приведенной частичной энтропии от среднего тока разряда.

Треугольник – режим, соответствующий максимальному току.

В качестве масштаба выбрана энтропия невозмущенного разрядом потока.

Римскими цифрами на рисунке отмечены стадии возмущения давления, где I – область, соответствующая неподвижному воздуху, II – область невозмущенного разрядом течении, III – область линейного возмущения плазмой, IV – область сильного возмущении.

Существование последней области определяется началом экзотермических реакций, индуцированных плазмой (плазменное воспламенение). Время существования каждой из стадий соответствует TII=0.1 с, TIII=0.2 с, TIV=0.65 с. Отметим, что величина обратная времени существования каждой из стадий меньше максимальной частоты макроскопических пульсаций. При этом на основании рис. 11 можно сделать вывод о влиянии разряда, существующего в режиме максимального тока на свойства натекающего потока. Как показывает закономерность, представленная на данном рисунке, число элементарных масштабов возмущений (степеней свободы) сокращается, т.е. поток поступающей в область разряда воздушно – пропановой смеси приближается к ламинарному состоянию - стабилизируется. При этом частичная энтропия падает до значения S/S0=10.55. меньшего значения S/S0=22.56, соответствующего минимальному разрядному току I=5.44 А.

Рис. 12. Сигнал датчика статического давления трубки Пито – Прандтля, расположенной в конце первой секции канала. Средний ток I=15.6 А.

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

Изменения продольной координаты электродов позволило усилить влияние обратного течения, дестабилизирующего поток и неоднородности состава смеси. В середине верхней стенки разрядной камеры было измерено статическое давление. Пример сигнала давления представлен на рис. 13. Рассчитано отношение K совокупного времени стадии воспламенения II к длительности разряда. Данная величина может быть интерпретирована как стабильность воспламенения.

Рис.13. Сигнал датчика статического давления датчика, расположенного в середине верхней стенки разрядной секции. II – плазменное воспламенение, I – режим нормального взаимодействия разряда с потоком.

Рис.14. Зависимость коэффициента стабильности плазменного воспламенения от среднего разрядного тока.

Зависимость стабильности от разрядного тока представлено на рис.14. Диаграмма иллюстрирует наличие двух предельных режимов, в области между которыми зафиксирован рост стабильности плазменного воспламенения от среднего тока разряда.

В пятой главе диссертации приведены результаты исследования взаимодействия плазмы импульсного эрозионного плазмотрона с натекающим воздушно – пропановым потоком в экспериментальном канале типа ЭУ2. Схема конструкции и питания плазматрона, представлена на рис. 15. Напряжение источника питания составляло 3,5 кВ, емкость накопительного конденсатора (С) была 50 µФ. Квазипериод разрядного тока плазмотрона, определенный при помощи пояса Роговского составлял T=50 мкс при полном времени выделения энергии =100 мкс. Схема размещения плазмотрона в канале представлена на рис. 15 б. Плазменная струя, созданная плазмотроном, распространялась навстречу потоку. Угол между осью симметрии плазматрона и направлением потоком составлял 30°. Интегральные фотографии распространения плазменной струи в первой секции канала, полученные с различной экспозицией, позволили определить динамические параметры области свечения плазмы (рис. 16 а). Предполагается, что пространственный профиль концентрации возбужденных частиц остается неизменным.

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

а) б) Рис. 15. а) Электрическая схема разряда. 1 –анод, 2 – диэлектрик, 3 – медная шайба, 4 – катод, 5 – источник питания (ИП) постоянного тока, 6 – разрядник. б) Схема размещения разряда в потоке. 1 –плазматрон, 2 – ось симметрии, 3 – зона циркуляции, 4 – ступень.

Размеры указаны в миллиметрах.

а) б) Рис.16. а) Фотография плазменной струи в первой секции канала, выполненная с экспозицией 45 мкс. б) Зависимость скорости передней границы плазмы (сплошная линия) и максимальной ширины плазмы (треугольники от времени).

Анализ интегральных фотографий позволил построить зависимость проекции скорости передней границы плазмы на ось плазматрона от времени (рис. 16 а, сплошная линия). Максимальная скорость оказывается равной V=3,1 км/с и соответствует времени t=20 мкс. Скорость боковой границы достигает при этом V=5 км/c и имеет максимум в тот же момент времени.

Показанные на рисунке 16 б экспериментальные точки, обозначенные как треугольники, соответствуют максимальной ширине плазмы. Как видно из диаграммы процесс расширения завершается в момент времени t=30 мкс, после чего начинается компрессия плазменного формирования потоком и его дальнейший конвективный снос.

На рис. 17 показаны сигналы с датчиков давления при импульсном разряде в пропан – воздушной смеси. В этом случае регистрируются сигналы давления, связанные с распространением плазменного образования (I и II). В определенный момент времени датчики регистрируют резкое возрастание давления в канале. При этом импульс повышенного давления распространяется как по потоку (положительное направление), так и против потока (отрицательное направление, прямая III на рис. 17-б).

а) б) Рис. 17. Сигналы с датчиков давления при импульсном разряде в потоке топливной смеси, а - в панорамном, б - в подробном временных масштабах. Плазматрон расположен в области датчика № 15, датчики с номерами меньше 15 расположены против потока, а с номерами больше 15 – по потоку.

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

На рис. 18 а, приведен временной ход сигналов датчиков давления, установленных в различных секциях канала в середине верхней стенки. Амплитуда сигналов выражена в условных единицах давления. Скачки давления сопровождаются возмущениями плотности. Для анализа данных возмущений были проведены рефракционные измерения в 4-ой секции канала (луч лазер проходил через центр симметрии прямоугольного иллюминатора), а также измерения, выполненные при помощи фотоумножителя (ФЭУ).

На рис. рис.18 б. и рис.18 в. представлены сигналы с рефракционного датчика и с ФЭУ, соответственно.

а) б) в) Рис. 18. а) Сигнал возмущения давления в четырех секциях канала;

б) Пример сигнала возмущения плотности, полученного рефракционным методом;

в) Пример сигнала ФЭУ, расположенного перед иллюминатором четвертой секции канала. Прямоугольниками выделены области неоднородности плотности и яркости.

Сравнение трех диаграмм рис. 18 позволяет сделать вывод о том, что по направлению первоначального потока топливной смеси распространяется плазменная область, на границе которой происходит плазменное воспламенение воздушно – пропановой смеси. Анализ динамики импульса давления, вызванного распространением плазмы, показывает, что амплитуда импульса не является затухающей только в присутствии пропана, при этом средняя продольная скорость распространения возмущения I (рис. 17 б) в четырех секциях канала составляет 681 м/c.

При помощи трубки Пито – Прандтля было определено число Маха перед фронтом возмущения давления I. Трубка была расположена на оси симметрии цилиндрического переходного фланца, подключенного к камере низкого давления. Было, получено, что скачок давления распространяется с дозвуковой скоростью, M1=0,54. Усреднение сигнала было проведено за 1 секунду. Статистическая ошибка определения числа Маха составляла M=0,1. На основании полученных параметров можно сделать вывод о том, что по направлению потока распространяется дозвуковая волна давления, интенсивность которой определяется сопутствующим фронтом реакции плазменного воспламенения.

В заключение работы сформулированы основные результаты диссертации:

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

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

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

Показана связь между нестационарными свойствами разряда и потока.

4. Исследование взаимодействия продольно – поперечного разряда постоянного тока со сверхзвуковым воздушно - пропановым потоком позволило определить дополнительные нестационарные свойства разряда при различных режимах его существования.

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

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

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

7. Исследование взаимодействия плазмодинамического импульсного разряда с натекающей воздушно – пропановой смесью позволило определить две стадии распространения разряда: стадия, при которой разряд не испытывает существенного влияния со стороны потока и стадия переноса плазменного следа. Зафиксировано плазменное воспламенение на границе плазмы, которое существует в течении всего времени сноса плазменного образования потоком.

СПИСОК ЦИТИРУЕМОЙ ЛИТЕРАТУРЫ.

1. Г. Н. Абрамович. Прикладная газовая динамика//Наука, 1976.

2. Ершов А.П., Колесников Е. Б., Тимофеев И. Б., Черников В. А., Чувашев С. Н., Шибков В. М., Плазмодинамические разряды в поперечных сверхзвуковых потоках воздуха // ТВТ, Т.44., № 4, 2006.

3. Райзер Ю.П., Физика газового разряда//Наука, 1992.

4.В. Бутаков, А. Граковский, Оценка уровня стохастичности временных рядов произвольного происхождения при помощи показателя Херста//Computer modeling and new technologies, Vol.9, No.2, 2005.

5.Г. М. Заславский, Р. З. Сагдеев, Введение в нелинейную физику//Наука, 1988.

ПУБЛИКАЦИИ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ.

1. V.Chernikov, E.Kolesnikov, S.Kamenshchikov, Explosive combustion, initiated by plasmodynamic discharge in propane - air mixture//43rd AIAA Aerospace Sciences Meeting and Exhibit, USA, Orlando, 2010.

2. А.П. Ершов, С.А. Каменщиков, Е.Б. Колесников, А.А. Логунов, А.А. Фирсов, В.А.

Черников, Измерение скорости потока с помощью поперечного разряда//Вестник МГУ. Серия 3. Физика, Астрономия, № 3, 2008.

3. А.Ф. Александров, А.П. Ершов, С.А. Каменщиков, А.А. Логунов, В.А. Черников, Воспламенение сверхзвуковой пропан-воздушной смеси с использованием импульсной плазмы//Вестник МГУ. Серия 3. Физика, Астрономия, № 2, 2008.

4. А.П. Ершов, С.А. Каменщиков, Е.Б. Колесников, А.А. Логунов, А.А. Фирсов, В.А.

Черников, О возможности измерения скорости потока с помощью маломощного импульсно-периодического разряда//Изв. РАН. Механика жидкости и газа, № 4, 2008.

5. А. П. Ершов, С. А. Каменщиков, А. А. Логунов, В. А. Черников, Инициация горения сверхзвукового пропан - воздушного потока разрядом магнитоплазменного компрессора//ТВТ, том 47, № 6, 2009.

6. А. П. Ершов, С. А. Каменщиков, А. А. Логунов, В. А. Черников, Воспламенение высокоскоростного воздушно-пропанового потока, инициируемое продольно поперечным разрядом постоянного тока//ТВТ, том 47, № 6, 2009.

7. Афонина Н.Е., Громов В.Г., Ершов А.П., Каменщиков С.А., Черников В.А., Воспламенение высокоскоростного воздушно-пропанового потока, инициируемое разрядом постоянного тока: эксперимент и численное моделирование // Нелинейный мир, № 11, 2009.



 

Похожие работы:





 
2013 www.netess.ru - «Бесплатная библиотека авторефератов кандидатских и докторских диссертаций»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.