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

Клементьева, Ирина Борисовна

Особенности взаимодействия электрических разрядов с газовыми потоками во внешних магнитных полях применительно к задачам интенсификации смешения тема автореферата и диссертации по физике, 01.04.08 ВАК РФ

Клементьева, Ирина Борисовна АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Москва МЕСТО ЗАЩИТЫ

2007 ГОД ЗАЩИТЫ

01.04.08 КОД ВАК РФ

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

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

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

КЛЕМЕНТЬЕВА Ирина Борисовна

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

Специальность 01 04 08 - физика плазмы

Автореферат диссертации на соискание учёной степени кандидата физико-математических наук

□ОЗ1742 Ю

Москва - 2007

003174210

Работа выполнена в Объединенном институте высоких температур РАН

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

доктор физико-математических наук В А Битюрин

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

доктор физико-математических наук В М Шибков,

доктор физико-математических наук А Б Ватажин

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

Физико-технический институт им Иоффе РАН

Защита состоится «

нлжйш 2007 г в^

часов на заседании

диссертационного совета Д 00^110 02 в Объединенном институте высоких температур Российской академии наук по адресу 125412, Москва, ул Ижорская, 13/19, Экспозиционный зал

С диссертацией можно ознакомиться в библиотеке ОИВТ РАН

Ваш отзыв на автореферат в 2-х экз , заверенный печатью организации, просим выслать по адресу 125412, Москва, ул Ижорская, 13/19, ОИВТ РАН, ученому секретарю диссертационного совета Д 002 110 02

Телефон для справок (495) 485-26-38

Учёный секретарь диссертационного совет Доктор физико-математических наук

Хомкин

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

Актуальность темы

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

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

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

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

Применение МГД-метода интенсификации базируется на концепции реакционного объема Горение как элементарная химическая реакция может произойти только в объеме, где имеет место столкновение молекул топлива и окислителя Этот объем носит название "Реакционный объем" и формируется в результате диффузии на контактной поверхности струй топлива и окислителя V(t)=S(f) ö(t), S(f) - площадь контактной поверхности струй топлива и окислителя, 6(0 - толщина диффузионного слоя Концепция реакционного объема сформулирована, описана и апробирована в работах [1, 2]

Молекулярная диффузия достаточно медленный процесс на масштабах доступных для смешения времен, лимитированных временем пребывания топлива и окислителя в камере сгорания 5 &-J~Dt, D и f - коэффициент диффузии и время Однако ускорить процесс создания реакционного объема можно увеличивая площадь контактной поверхности, например, генерацией вихревых возмущений в пото-

ках топлива и окислителя (кинематическое перемешивание) Как показано в работе [3] и продемонстрировано в [4], в завихренном турбулентном потоке длины материальных линий и площади материальных поверхностей, в нашем случае площадь контактной поверхности топлива и окислителя, растут со временем по экспоненциальному закону 5~5оехр {/£2}, Йо - начальное значение площади контактной поверхности, О - масштаб характерной неоднородности завихренного потока

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

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

Процессы и явления, обусловленные применением МГД-метода интенсификации смешения, изучены недостаточно полно, и их экспериментальное и расчетно-теоретическое исследование актуально для физики плазмы и плазменной аэродинамики

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

Научная новизна работы состоит в следующем

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

2 Получены экспериментальные результаты по эволюции и динамике электрического разряда в газовой среде/потоках газа во внешнем магнитном поле в широком диапазоне параметров

3 Проведено численное исследование эволюции и динамики электрического разряда в газовой среде/потоках газа во внешнем магнитном поле и влияния взаимодействия электрического разряда с внешним магнитным полем на структуру и смешение потоков газа

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

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

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

7 Продемонстрирована эффективность применения МГД-метода интенсификации смешения предварительно неперемешанных газовых струй, через характер и особенности эволюции и динамики электрического разряда в газовой среде/потоках газа во внешнем магнитном поле, структуру потоков и оценку величины реакционного объема

Достоверность результатов

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

Научная и практическая значимость

В результате экспериментальных и расчетно-теоретических исследований получены новые сведения об особенностях взаимодействия контрагированных электрических разрядов с газовыми потоками во внешнем магнитном поле и о влиянии особенностей этого взаимодействия на смешение предварительно неперемешанных газовых струй Результаты исследований демонстрируют эффективность применения МГД-метода интенсификации смешения и необходимы при его практической реализации

Положения, вынесенные на защиту

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

2 Квазитрехмерная трехтемпературная модель контрагированного электрического разряда в потоках газа во внешнем магнитном поле

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

4 Анализ результатов исследований параметров электрического разряда, его эволюции и динамики в газовой среде/потоках газа во внешнем магнитном поле,

влияния взаимодействия электрического разряда с внешним магнитным полем на структуру и смешение потоков газа, полученных с использованием пакета программ двумерного моделирования PLASM AERO, разработанного к ф -м н АН Бочаровым

Апробация работы

Основные результаты работы докладывались и обсуждались на следующих конференциях 4-th Workshop on Magnetoplasma Aerodynamics in Aerospace Applications, 9-11 April, 2002, 14th International Conference on MHD Power Generation and High Temperature Technologies, Maui, Hawaii, 2002, 41st, 43rd, 44th, 45th Aerospace Sciences Meeting and Exhibit, Reno, NV, 2003 - 2007, 15th International Conference on MHD Energy Conversion and 6th International Workshop on Magnetoplasma Aerodynamics, Moscow, May 24-27, 2005, 27-th International Conference on Phenomena in Ionized Gases (ICPIG), Eindhoven, the Netherlands, 18-22 July, 2005, Всероссийская научно-техническая конференция молодых ученых и специалистов «проблемы создания перспективных авиационных двигателей», Москва, ЦИАМ, 2730 сентября, 2005, Минский международный коллоквиум по физике ударных волн, горения, взрыва, детонации и неравновесным процессам MIC 2005, Минск, 14-19 ноября 2005, Школа-семинар по магнитоплазменной аэродинамике, 26-27 апреля, Москва, 2006, 31st International symposium on combustion, Heidelberg, Germany, August 611, 2006, 16th International Conference on Gas Discharges and their Applications, Xian, China, Sept 11-15, 2006, 7th International Workshop on Magnetoplasma Aerodynamics, Moscow, May 24-27, 2007, 16th International Conference on MHD Energy Conversion, 25-28 June 2007, Miami, Florida

Публикации

Результаты диссертационной работы опубликованы в 14 научных работах, список которых приводится в конце автореферата

Объем работы

Диссертационная работа состоит из введения, четырех глав, заключения и списка литературы Объём диссертации составляет 130 страниц, включая 95 рисунков

СОДЕРЖАНИЕ РАБОТЫ

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

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

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

Во второй главе приведена принципиальная схема экспериментов и описана экспериментальная установка, созданная автором для исследования эволюции и динамики электрического разряда в газовой среде/потоках газа во внешнем магнитном поле и МГД стимулированного смешения и поджига. Представлены результаты экспериментальных исследований эволюции, динамики и электрических характеристик разряда в широком диапазоне параметров для нескольких конфигураций экспериментальной секции.

Принципиальная схема экспериментов показана на рис. 2.1.

Основной поток 1, моделирующий окислитель, поступает в экспериментальную секцию слева направо. Встречно/спутно основному потоку через расположенную на оси аэродинамического тракта систему инжекции подается газ 2, моделирующий топливо. Электрический разряд 6, инициируемый в магнитном поле 3 между проволочным 4 и кольцевым 5 электродами, под действием электромагнитной силы Р деформируется, перемещается и создаёт завихренности в потоке, что является интенсифицирующим фактором перемешивания топлива и окислителя.

Схема экспериментальной установки для изучения эволюции и динамики разряда в газовой среде/потоках газа во внешнем магнитном поле и МГД стимулированного смешения и поджига представлена на рис. 2.2.

Эволюция и динамика электрического разряда атмосферного давления в магнитном поле в газовой среде в конфигурации рис. 2.1. Параметры проведения экспериментов: ток короткого замыкания цепи /кз = 2,8; 1,4; 0,7; 0,35; 0,175 А; величина индукции магнитного поля 6=0,01 - 0,05 Тл; напряжение на разряднике {А* =7,02 кВ; давление Р=105 Па.

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

Рис. 2.1. Принципиальная схема экспериментов:

1 - основной поток; 2 - инжектируемый поток; 3 — линии индукции магнитного поля; 4 - проволочный электрод; 5 - кольцевой электрод; 6 -электрический разряд

Рис. 2.2. Схема экспериментальной установки:

I - проволочный электрод; 2 - кольцевой электрод; 3 - генератор; 4 - магнитная катушка; 5 - система включения магнитного поля и подачи сигнала на генератор импульсов; 6 -аккумуляторная батарея; 7- реостат; 8 - кварцевая труба; 9 - ресивер; 10 - компрессор;

II - система инжекции; 12 - баллон с газом; 13 - скоростной фоторегистратор; 14 - компьютер; 15-зеркало; 16-схема измерения тока в цепи /агс и напряжения на дуге (Лгс; 17 - токовый шунт; 18, 19- расходомеры; 20 - осциллограф; 21 - генератор импульсов

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

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

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

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

Рис. 2.3. Форма разряда в магнитном поле и перезамыкание на внешний электрод, /а,с=2,8 А, 8=0,05 Тл

разряда, движущейся преимущественно в радиальном направлении к кольцевому электроду, зафиксированного на фотоснимке рис 2 3

Видимое вращение разряда есть результат его эволюции под действием электромагнитной силы в рассматриваемой конфигурации и перезамыкания внешней части токового канала на кольцевой электрод (см рис 2 3)

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

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

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

Установлено, что при низких значениях электрического тока и -индукции магнитного поля (в экспериментах с /«э=0,35, 0,175 А, В=0,01 Тл) на динамику разряда существенное влияние оказывает сила Архимеда

В качестве особенности поведения разряда в магнитном поле отмечено возможное явление повторного пробоя на внешний электрод - катод, обусловленное

распределением электрического поля в рассматриваемой области. Положение повторного пробоя, как правило, локализовано в месте привязки катодного пятна, что непосредственно видно на фотоснимках рис. 2.4.

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

Экспериментальные результаты эволюции и динамики разряда представлены на рис. 2.5, 2.6, где частоты вращения разрядов для двух режимов построены как функции тока /кз_ и индукции магнитного поля В.

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

Рис. 2.5. Зависимость частоты вращения разряда от тока /кз. для нескольких значений величины магнитной индукции В

Рис. 2.4. Перепробой разряда в процессе его эволюции на внешний катод для /„.3=2,8; 0,7; 0,175 А и В=0,05, 0,01 Тл соответственно

1000 800 £ 600

Ф 400

£ 200

Anode

600 1.4А

500 0.7А 0.35А

N 400 —•—0.175А

—2.8А

g 300 ф 1 200 /

iL 100 /

Cathode

0,00 0,01 0,02 0,03 0,04 0,05 0.0

в.т

0,01 0,02 0,03 0.04 0,05 0,06

в,т

Рис. 2.6. Зависимость частоты вращения разряда от величины магнитной индукции В для нескольких значений тока /„.3.

Вольт-амперные характеристики разряда для нескольких значений величины магнитной индукции В приведены на рис. 2.7. Ток и напряжение дуги измерены при

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

Рис. 2.7. ВАХ разряда для нескольких значений величины магнитной индукции В

Эволюция и динамика электрического разряда давления ниже атмосферного Р=1,5 10* - 2,5-10* Па в магнитном поле в газовой среде в конфигурации рис. 2.1. Для исследования эволюции и динамики электрического разряда в магнитном поле при низких давлениях экспериментальная секция монтируется в вакуумной камере. Параметры проведения экспериментов: ток короткого замыкания цепи /к.3=1 А; величина индукции магнитного поля В=0,1 - 0,24 Тл; напряжение на разряднике (Л.х.= 1 кВ; давление в барокамере Р=1,5Ю4 - 2,5-Ю4 Па. Как и в случае дуги атмосферного давления, в экспериментах зафиксированы эволюция раз-

20 i i i i

СО

Е *

о 10

о *

0.01 0.02

в, т

Рис. 2.8. Скорость перемещения

элемента разряда /агс=0,9 А

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

Обнаружено, что при величине индукции магнитного поля £¡=0,03 Тл вращение разряда не происходит, катодное пятно привязано к определенной точке кольцевого электрода в течение всего времени наблюдения. Для магнитного поля В=0,01 Тл отмечено ярко выраженное попятное движения катодного пятна.

Экспериментально определенны скорости перемещения приэлектродной части разряда у центрального анода для двух значений величины индукции магнитного поля В=0,1 и 0,2 Тл для давления Р=0,2-105 Па. Показано, что скорость движения дуги хорошо коррелирует с квадратичным (по скорости) законом сопротивления движения тел в сплошной среде (рис. 2.8).

Эволюция и динамика электрического разряда атмосферного давления в магнитном поле в потоке газа скорости 0,2 - 10 м/с в конфигурации рис. 2.1. Одной из особенностей поведения разряда в рассматриваемых условиях является его снос газовым потоком вдоль электродов и разрыв токовой цепи. Параметры проведения экспериментов: ток короткого замыкания цепи /к.3=1; 0,7 А; величина индукции магнитного поля 8=0,01 - 0,24 Тл; напряжение на разряднике (Л.х.= 1; 7,02 кВ; давление Р=105 Па, скорость воздушного потока - 0,2 - 10 м/с.

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

Рис. 2.9. Кадры скоростной видеосъемки взаимодействия электрического разряда с внешним магнитным полем В=0,05 Тл в воздушном потоке 10 м/с (кадры представлены через каждые 0,42 мс). Внешний электрод - катод, /к.3=1 А. Как и в случае экспериментов без потока, наблюдается привязка катодного пятна к определенной точке кольцевого электрода в течение некоторого времени

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

Динамика электрического разряда атмосферного давления в потоке газа на рельсовом электроде в конфигурации рис. 2.10.

Рис. 2.10. Экспериментальная секция Рис 211 интегральные фотографии ди-

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

вый электрод но. Скорость потока 1,2 м/с

В рассматриваемой конфигурации разряд, будучи прикрепленным одним концом к центральному электроду, скользит другим концом вдоль рельсового в воздушном потоке. Параметры проведения экспериментов: ток короткого замыкания цепи 'к.з.=1 А; напряжение на разряднике ¿Ухх.=1 кВ; давление Р=105 Па, скорость воздушного потока - 1,2 - 10,4 м/с (на входе в тестовую секцию).

В случае конфигурации рельсовый электрод - катод разряд непрерывно вытягивается вдоль потока рис. 2.11, а при достижении некоторой длины (-10 см в рассматриваемых условиях), определяемой мощностью источника, перепробивается в зазоре между электродами, и процесс его вытяжения повторяется. Причиной вторичного пробоя в этом случае является потеря устойчивости разряда, обусловленная внешней электрической цепью, которая перестает работать в режиме генератора тока с ростом сопротивления вытягивающейся дуги. Для катодного пятна характерно перемещение вдоль электрода со скоростью меньшей скорости потока: 0,8-Ю"3 - 1,7-10"3 м/с для скоростей потока 1,2 - 10,4 м/с соответственно. То, что скорость катодного пятна ниже, чем скорость потока, по-видимому, говорит о тепловом механизме его перемещения вдоль рельсового катода.

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

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

Эволюция электрического разряда атмосферного давления во встречных потоках газа во внешнем магнитном поле в конфигурации рис. 2.1. Параметры проведения экспериментов: ток короткого замыкания цепи /хз=1 А; напряжение на разряднике ихх = 1 кВ; величина индукции магнитного поля В=0,1 Тл; давление Р=105 Па; скорость основного воздушного потока - 1 - 3,5 м/с; в качестве газа встречного потока используется азот, скорость встречного потока - до 30 м/с.

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

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

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

Квазитрехмерная трехтемпературная модель электрического разряда в газовой среде/потоках газа во внешнем магнитном поле. Особенность квазитрехмерной модели заключается в том, что она строится на соотношениях макробалланса [2]. Источниками основных характеристик и зависимостей для построения модели служат результаты физических и двумерных численных экспериментов. Не претендуя на описание начальной стадии разряда, предполагается, что дальнейшее его развитие может быть описано в рамках модели трехтемпературной плазмы. Модель предполагает, что тяжелые атомы, молекулы и ионы имеют одина-

Magnetic field on

J \/riHnno

j-1 Voltage

Рис. 2.12. Эволюция разряда в магнитном поле

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

При построении физической модели электрического разряда используются следующие положения

1) масса вещества разряда считается постоянной во времени величиной,

2) рассматриваются осредненные по сечению разряда характеристики,

3) учитывается сильная неоднородность характеристик разряда по длине,

4) описание потоков через разрядную поверхность осуществляется на основе

гипотезы диффузионного слоя,

5) ускорение разряда определяется как собственной массой канала, так и присоединенной массой обтекающего его внешнего потока

Для определения характеристик разряда и его положения в пространстве в каждый момент времени решается система уравнений, записанная для каждого малого элемента объема разряда V

— \р (IV = М = сог^, ей

(1)

+ V,

М^ + а..АЗ = ду V, V,

М -=^хВ V,

с1е.

(4)

(3)

(2)

(5)

(6)

(7)

Здесь р, М - плотность и масса газа элемента электрической дуги, е0. ее, е„ -осредненные по массе элемента полная энергия, энергия электронного газа и колебательная энергия возбужденных двухатомных молекул Слагаемые до, де, <?У - полные потоки тепла через единичную поверхность элемента разряда Оа, Qe, имеют смысл осредненных по объему элемента источниковых членов - полный энерговклад, энерговкпад в электронный газ и в колебательные степени свободы тяжелого газа д - плотность ./-той компоненты В уравнении (7) г и и - радиус-векторы и скорости центров масс элементов разряда, при этом и=\/в+\/Яо„(г), где У в - скорость движения, обусловленная действием электромагнитной силы на проводник с током в магнитном поле, УтмАт) - скоростное поле невозмущенного потока В уравнении (8) Р - давление газа, представляемое суммой парциальных давлений компонент

, Я и к - универсальная газовая постоянная и постоянная

7=1

Больцмана п, - концентрации соответствующих компонент

Связь между элементами разряда осуществляется из условия протекания общего тока, определяемого уравнением внешней цепи Е=1(Я0+Яагс), где Е - ЭДС источника и /?0 - внешнее сопротивление цепи - заданные параметры, а сопротивление электрического разряда Шагс =-———

<т(Те)лг2

Постановка задачи Для исследования поведения электрического разряда в потоке во внешнем магнитном поле решается модельная задача в постановке рис 2 1 Расчетной областью является пространство, ограниченное цилиндрической поверхностью Скоростное поле газового потока УлаЛт) представлено как суперпозиция двух полей Уяои^Л+УгМ, где У^сопв! и \/2=01Ля2 - величины скоростей основного (течет слева направо) и инжектируемого потоков О - расход жидкости через сферическую поверхность радиуса г Электрический разряд инициируется между проволочным электродом, расположенным на оси, и электродом внутренней поверхности цилиндра, расстояние между которыми 30 мм Под действием электромагнитной силы и газового потока разряд перемещается и деформируется Данная схема обеспечивает локализацию разряда вблизи контактной поверхности основного и встречного потоков Можно ожидать, что локализация разряда вблизи их контактной поверхности решит задачу интенсификации смешения и горения в спутных струях

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

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

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

Time, s X, cm

Рис. 3.1. Температуры в разряде /ка.= 5 А Рис. 3.2 Эволюция разряда в поле

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

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

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

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

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

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

В четвертой главе представлен анализ результатов исследования параметров электрического разряда, его эволюции и динамики в газовой среде/потоках газа во внешнем магнитном поле, влияния взаимодействия электрического разряда с внешним магнитным полем на структуру и смешение потоков газа. Результаты получены с использованием пакета программ двумерного моделирования PLASM AERO, разработанного к.ф.-м.н. А Н. Бочаровым.

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

Рис. 4.1. Температурное поле из разряда не исчезает, то каждый период разряд

Ту

рассекающей его цилиндрической поверхности, расположенной достаточно далеко от электродов.

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

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

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

Численное моделирование эволюции и динамики разряда в газовой среде во внешнем магнитном поле в конфигурации коаксиапьно расположенных электродов. Для исследования особенностей эволюции электрического разряда в газовой среде во внешнем магнитном поле решается задача в постановке рис. 2.1. На рис. 4.3, 4.4 для нескольких последовательных моментов времени представлена эволюция электрического разряда /к.3 = 5 А в поле Р= 1 Тл через температурное поле, мольную долю электронов и показаны профили скорости газа в кольцевом зазоре между электродами. На начальной стадии (<30 мкс) преобладает механизм теплового расширения столба разряда вследствие джоулева тепловыделения. В то же самое время электромагнитная сила, действуя по часовой стрелке (ток течет от центрального электрода к кольцевому), вызывает движение газа в том же направлении. Условие прилипания на поверхности электродов приводит к искривлению токового канала у центрального электрода, что в свою очередь определяет изменение направления действия силы. Как и в квазитрехмерном расчете, вращение токового канала вокруг центрального электрода происходит из-за постоянного изменения направления действия электромагнитной силы. Точка привязки токового канала (пятно размером порядка 0.5 мм, через которое протекает большая часть тока) к центральному электроду также вращается. В пограничном слое у поверхности центрального электрода, где конвективное движение невелико, единственным механизмом, обеспечивающим перемещение пятна, может являться лишь прогрев газа перед токовым каналом. Также на рис. 4.3, 4.4 можно заметить явление перезамыкания разряда на центральный электрод. Скорости вращения разряда вокруг центрального и внешнего электродов отличаются примерно в 25 раз. Механизм вращения точки привязки то-

0.046 -

0.0Е+0 2.5Е-4 5.0Е-4 7.5Е-4 1.0Е-3 Time, sec

Рис. 4.2. Масса разряда и координата вдоль оси X, /кз = 5 А, В=1 Тл

кового канала на внешнем электроде в общем тот же самый, что и на внутреннем. Большое различие в скоростях перемещения обусловлено действием двух факторов: электромагнитной силой и полем давления в кольцевом зазоре. Важным выводом из результатов является то, что в движение приводится не только область разряда, но и холодный невозмущенный газ. Хорошо виден спиралевидный характер разряда, приобретаемый им в процессе движения в магнитном поле. Внешняя часть разрядной спирали перемещается преимущественно в направлении кольцевого электрода под действием электромагнитной силы, в то время как остальная часть вращается вокруг центрального. Область у центрального электрода с течением времени становится все более однородной, можно предположить, что именно здесь про-исходит наиболее интенсивное смешение. В центральной части характерная ско-рость ~500 м/с, в то время как у внешнего электрода среда почти неподвижна. В про-цессе раскручивания разряда средняя скорость жидкости в зазоре возрастает. И в конце наблюдения весь первоначально холодный и неподвижный объем вовлечен в интенсивное движение, вращение.

Деформация контактной поверхности воздуха и водорода электрическим разрядом, эволюционирующим в конфигурации коаксиально расположенных электродов в момент времени (=0,4мс и распределение электрического потенциала представлены на рис. 4.5, 4.6. Видно, что электрический ток растекается в двух направлениях вокруг центрального электрода и снова концентрируется у кольцевого.

©о

Рис. 4.3. Поле температуры и скорости Рис. 4.4. Мольная доля электронов

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

Рис. 4.5. Деформация контактной поверхнос- Рис. 4.6. Температура и мольная доля ти и распределение потенциала воды

Исследование влияния взаимодействия электрического разряда с внешним магнитным полем на структуру и смешение потоков газа. Постановка задачи схематично изображена на рис. 4.7. На рис. 4.8 приведены поля температуры и скорости в расчетной области для нескольких последовательных моментов времени эволюции разряда тока /агс=2 А в магнитном поле 8=1 Тп. Результат смешения представлен на рис. 4.9. Показано, что движение электрического разряда под действием электромагнитной силы сопровождается генерацией системы вихрей вокруг проводящего ядра, которые интенсифицируют перенос тепла и компонентов из проводящей области в окружающий поток. С другой стороны холодный и почти невозмущенный газ поступает в зону проводящего ядра. Таким образом, проводящее ядро работает как насос, который вовлекает свежую неперемешанную жидкость в процесс смешения.

На рис. 4.10 приведены зависимости координаты разряда от времени для токов /агс=1; 2; 4; 8 А соответственно.

Рис. 4.7.Схема расчета

•ч;

Рис. 4.8. Поля температуры и скорости Рис. 4.9. Зона смешения

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

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

0.00 0.20 0.40 0.60 0.80 time, ms

Рис. 4.10. Положение разряда во времени

0.20 0.40 0.60 time, ms

Рис. 4.11. Зависимость величины реакционного объема от времени

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

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

1. Создана экспериментальная установка для исследования эволюции и динамики электрического разряда в газовой среде/потоках газа во внешнем магнитном поле и исследования эффективности МГД-метода интенсификации смешения предварительно неперемешанных газовых струй. В этой связи апробированы различные варианты экспериментальной секции для исследования эволюции и динамики электрического разряда в газовой среде/потоках газа во внешнем магнитном поле.

2. Разработана квазитрехмерная трехтемпературная модель контрагирован-ного электрического разряда в потоках газа во внешнем магнитном поле.

3. Выполнен теоретический анализ эволюции и динамики разряда в магнитном поле в газовой среде и влияния взаимодействия электрического разряда с внешним магнитным полем на структуру потока.

4. Проведено экспериментальное исследование эволюции и динамики электрического разряда в газовой среде/потоках газа во внешнем магнитном поле.

5. На базе разработанной квазитрехмерной модели проведено расчетно-теоретическое исследование эволюции, динамики и параметров электрического разряда в газовой среде/потоках газа во внешнем магнитном поле.

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

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

РЕЗУЛЬТАТЫ ДИССЕРТАЦИИ ИЗЛОЖЕНЫ В СЛЕДУЮЩИХ РАБОТАХ

1 ANBocharov, SBLeonov, DSBaranov, IВ Klement'eva, VABityurin MHD Enhanced Mixing and Combustion in Cocurrent Streams // In 4th Workshop on Magnetoplasma Aerodynamics for Aerospace Applications, Moscow, Russia, 9-11 April, 2002, pp 220 - 230

2 A Bocharov, V Bityunn, I Klement'eva, and S Leonov Experimental and Theoretical Study of MHD Assisted Mixing and Ignition in Co-Flow Streams // In 14th International Conference on MHD Power Generation and High Temperature Technologies, Maui, Hawaii, 2002, AIAA Paper 2002-2248

3 ABocharov, V Bityunn, I Klement'eva, and S Leonov A Study of MHD Assisted Mixing and Combustion II In 41st Aerospace Sciences Meeting and Exhibit, Reno, NV, 2003, AIAA Paper 2003-5878

4 A Bocharov, I Klement'eva, A Klimov, V Bityunn A Study of MHD Assisted Mixing and Combustion in Couter-Flow Streams // In 43rd Aerospace Sciences Meeting and Exhibit, Reno, NV, 2005, AIAA Paper 2005-0600

5 I Klement'eva, A Bocharov, V Bityunn, A Klimov Experimental and Numencal Study of MHD Assisted Mixing // In 15th International Conference on MHD Energy Conversion and 6 International Workshop on Magnetoplasma Aerodynamics, Moscow, May 24-27, 2005, pp 365 - 374

6 I Klement'eva, ABocharov, V Bityunn, A Klimov MHD Assisted Mixing through the Experimental and Numerical Study // In 27th ICPIG, Eindhoven, the Netherlands, 1822 July, 2005, CD

7 И Б Клементьева, A H Бочаров, В А Битюрин, А И Климов Взаимодействие электрического разряда с газовым потоком во внешнем магнитном поле // Сб докладов Всероссийской научно-технической конференции молодых ученых и специалистов «проблемы создания перспективных авиационных двигателей», Москва, ЦИАМ, 27-30 сентября, 2005

8 A Bocharov, V Bityunn, I Klement'eva, A Klimov, Numencal and Expenmental Study of MHD Assisted Mixing and Combustion // In 44th Aerospace Sciences Meeting and Exhibit, Reno, NV, 2006, AIAA Paper 2006-1009

9 V Bityunn, I Klement'eva, A Klimov Mixing and combustion intensification of high speed reacting flows by the electrical discharges in magnetic field // In 31st International symposium on combustion, Heidelberg, Germany, August 6-11, 2006, CD

10 VBityunn, IKlementyeva, ABocharov Investigation of electrical discharges interaction with gas flows in external magnetic fields for problems of mixing and combustion

intensification // In 16th International Conference on Gas Discharges and their Applications, Xian, China, Sept 11-155, 2006, pp 425 - 428

11 A Bocharov, V Bityunn, I Klement'eva, A Klimov Numerical and Experimental Study of MHD Assisted Mixing and Combustion // In 45th Aerospace Sciences Meeting and Exhibit, Reno, NV, 2006, AIAA Paper 2007-1024

12 I Klementyeva, A Bocharov, V Bityunn Electrical Discharges - Gas Flows Interaction in External Magnetic Fields // In 7th International Workshop on Magnetoplasma Aerodynamics, Moscow, May 24-27, 2007, CD

13 I Klementyeva, V Bityunn and A Bocharov Constricted Discharge Interaction With High Speed Gas Flows // In 38th AIAA Plasmadynamics/16th MHD Energy Conversion Conference, Miami, FL, June 25-28, 2007, AIAA Paper 2007-4141

14 И Б Клементьева, АН Бочаров, В А Битюрин Особенности взаимодействия электрического разряда с газовым потоком во внешнем магнитном поле // Письма в ЖТФ Т 33 Вып 22, 2007, с 16 - 22

1 V A Bityunn, А N Bocharov Advanced MHD Assisted Mixing of Reacting Streams // In 39th Aerospace Sciences Meeting and Exhibit, Reno, NV, 2001, Paper 2001-0793

3 А С Монин, А МЯглом Статистическая гидромеханика, механика турбулентности - М Наука, 1967, с 513 - 517

4 V A Bityunn, AN Bocharov, V G Potebnya, N Yu Babaeva, V G Naidis, ALTseskis Modeling of Processes in GCB around Current Zero - M IVTAN - ANRA, 1999, p 61

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

КЛЕМЕНТЬЕВА Ирина Борисовна

Автореферат

Подписано в печать 11 09 07

Печать офсетная

Тираж 100 экз_

Уч -изд л 1 5 Заказ № У2Ч

Формат 60x84/16 Уел -печ л 1 39 Бесплатно

ОИВТ РАН 125412, Москва, ул Ижорская, 13/19

Содержание диссертации автор исследовательской работы: кандидата физико-математических наук, Клементьева, Ирина Борисовна

Введение

Глава 1. Обзор работ.

1.1. Разряды в потоках газа.

1.2. Электрический разряд в магнитном поле.

1.3 Влияние разрядов в магнитном поле на характеристики потоков.

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

Исследование взаимодействия поток - плазменное образование является перспективным и быстро развивающимся направлением современной науки. Необходимость более тщательного изучения такого взаимодействия возникла в связи с созданием высокоскоростных летательных аппаратов [1].

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

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

Одним из методов, позволяющим интенсифицировать смешение и горение в газовых потоках разделенных компонентов топлива и окислителя, является МГД-метод, который заключается: 1) в создании импульсного электрического разряда в газовом потоке; 2) в использовании взаимодействия электрического разряда с внешним магнитным полем для организации вторичных вихревых течений (рис.1).

Преимуществом МГД-метода интенсификации смешения является быстродействие и бесконтактное воздействие на потоки.

Применение МГД-метода интенсификации базируется на концепции реакционного объема. Горение как элементарная химическая реакция может произойти только в объеме, где имеет место столкновение молекул топлива и окислителя [5-10].

Этот объем носит название "Реакционного объёма" и формируется в результате диффузии на контактной поверхности струй топлива и окислителя: V(t)=S(t)-S(t), S(t) -площадь контактной поверхности струй топлива и окислителя, S(t) - толщина диффузионного слоя. Концепция реакционного объема сформулирована и описана в работах [11 -14] и апробирована в [15 -18].

Рис.1. Схема МГД-метода.

Рис.2. Деформация контактной поверхности в поле вторичных вихревых течений.

Молекулярная диффузия достаточно медленный процесс на масштабах доступных для смешения времен, лимитированных временем пребывания топлива и окислителя в камере сгорания: SvjDt, D и t - коэффициент диффузии и время. Однако, ускорить процесс создания реакционного объема можно увеличивая площадь контактной поверхности, например, генерацией вихревых возмущений в потоках топлива и окислителя (кинематическое перемешивание) [19 - 22]. Как показано в работах [23 - 25] и продемонстрировано в [11 - 14], в завихренном турбулентном потоке длины материальных линий (рис.2) и площади материальных поверхностей, в нашем случае площадь контактной поверхности топлива и окислителя, растут со временем по экспоненциальному закону: So - начальное значение площади контактной поверхности, - масштаб характерной неоднородности завихренного потока.

Для организации вторичных вихревых течений и используется взаимодействие тока электрических разрядов с внешним магнитным полем. Можно показать [26 - 28], что в среде с неоднородной проводимостью во внешнем магнитном поле непотенциальная, вихревая электромагнитная сила не может быть уравновешена градиентом давления, и в системе возникают вихревые течения (рис.3).

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

Цель работы

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

В соответствии с целями работы решены следующие задачи:

2. Разработана квазитрехмерная трехтемпературная модель контрагированного электрического разряда в потоках газа во внешнем магнитном поле.

6. С использованием пакета программ двумерного моделирования PlasmAero, разработанного к.т.н. А.Н. Бочаровым, проведены исследования параметров электрического разряда, его эволюции и динамики в газовой среде/потоках газа во внешнем магнитном поле, рассмотрено влияние взаимодействия электрического разряда с внешним магнитным полем на структуру и смешение потоков газа.

Научная новизна работы состоит в следующем:

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

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

3. Проведено численное исследование эволюции и динамики электрического разряда в газовой среде/потоках газа во внешнем магнитном поле и влияния взаимодействия электрического разряда с внешним магнитным полем на структуру и смешение потоков газа.

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

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

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

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

Достоверность результатов

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

Научная и практическая значимость

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

Адробадия работы

Основные результаты работы докладывались и обсуждались на следующих конференциях: 4-th Workshop on Magnetoplasma Aerodynamics in Aerospace Applications, 911 April, 2002; 14th International Conference on MHD Power Generation and High Temperature Technologies, Maui, Hawaii, 2002; 41st, 43rd, 44th, 45th Aerospace Sciences Meeting and Exhibit, Reno, NV, 2003 - 2007, 15th International Conference on MHD Energy Conversion and 6th International Workshop on Magnetoplasma Aerodynamics, Moscow, May 24-27, 2005; 27-th International Conference on Phenomena in Ionized Gases (ICPIG), Eindhoven, the Netherlands, 18-22 July, 2005; Всероссийская научно-техническая конференция молодых ученых и специалистов «проблемы создания перспективных авиационных двигателей», Москва, ЦИАМ, 27-30 сентября, 2005; Минский международный коллоквиум по физике ударных волн, горения, взрыва, детонации и неравновесным процессам MIC 2005, Минск, 14-19 но-ября 2005; Школа-семинар по магнитоплазменной аэродинамике, 26-27 апреля, Мо-сква, 2006; 31st International symposium on combustion, Heidelberg, Germany, August 6-11, 2006; 16th International Conference on Gas Discharges and their Applications, Xian, China, Sept. 11-15, 2006; 7th International Workshop on Magnetoplasma Aerodynamics, Moscow, May 24-27, 2007; 16th International Conference on MHD Energy Conversion, 25-28 June 2007, Miami, Florida.

Публикации

Результаты диссертационной работы опубликованы в 14 научных работах.

Положения, вынесенные на защиту:

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

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

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

4. Анализ результатов исследований параметров электрического разряда, его эволюции и динамики в газовой среде/потоках газа во внешнем магнитном поле, влияния взаимодействия электрического разряда с внешним магнитным полем на структуру и смешение потоков газа, полученных с использованием пакета программ двумерного моделирования PlasmAero, разработанного к.т.н. А.Н. Бочаровым.

Краткое содержание работы

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

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

В четвертой главе представлен анализ результатов исследования параметров электрического разряда, его эволюции и динамики в газовой среде/потоках газа во внешнем магнитном поле, влияния взаимодействия электрического разряда с внешним магнитным полем на структуру и смешение потоков газа. Результаты получены с использованием пакета программ двумерного моделирования PlasmAero, разработанного к.т.н. А.Н. Бочаровым.

Заключение диссертации по теме "Физика плазмы"

Общие выводы

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

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

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

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

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

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

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

3. С использованием пакета программ двумерного моделирования PlasmAero, разработанного к.т.н. А.Н. Бочаровым, проведены исследования параметров электрического разряда, его эволюции и динамики в газовой среде/потоках газа во внешнем магнитном поле, рассмотрено влияние взаимодействия электрического разряда с внешним магнитным полем на структуру и смешение потоков газа

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

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

4.5. Заключение

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

Список источников диссертации и автореферата по физике, кандидата физико-математических наук, Клементьева, Ирина Борисовна, Москва

1. W.H.Heiser, D.T.Pratt. Hypersonic Airbreathing Propulsion. Washington DC.: A1.A Education Series, 1994,588 p.

2. В.Е.Алемасов, А.Ф.Дрегалин, А.С.Черенков. Основы теории физико-химических процессов в тепловых двигателях и энергетических установках. М.: Химия, 2000.

3. Ф.Бартльме. Газодинамика горения. М.: Энергоиздат, 1981,280 с.

4. Ф.А.Вильямс. Теория горения. М.: Наука, 1971,616 с.

5. П.К.Третьяков. Газодинамика горения в сверхзвуковом потоке. Новосибирск: ИТПМ СО АН СССР, 1979,101 с.

6. Ю.М.Пчелкин. Камеры сгорания газотурбинных двигателей. М.: Машиностроение, 1973,392 с.

7. Е.С.Щетинков. Физика горения газов. М.: Наука, 1965.

8. V.A.Bityurin, V.G.Potebnya, A.L.Tseskis. Evolution of a Current Plasma Clot in Turbulent Flow // In: Proc. of 33rd SEAM Conf., Tennessee, June 12-15,1995, p.IV.7.

9. В.А.Битюрин, В.Г.Потебня, А.Л.Цескис. Об эволюции токонесущего плазменного сгустка в среде со случайным полем скоростей // ПЖТФ, 1996, т.22, No.2, с.80-83.

10. В.А.Битюрин, В.Г.Потебня, А.Л.Цескис. Перенос энергии в газоплазменном течении при наличии токопроводящих кластеров // Магнитная гидродинамика, 1997, т.ЗЗ, No.3, с.297-305.

11. VABityurin, A.N.Bocharov, V.G.Potebnya, N.Yu.Babaeva, V.G.Naidis, AX.Tseskis. Modeling of Processes in GCB around Current Zero // Research Report IVTAN ANRA No.99/26, June 1999,61 p.

12. V.A.Bityurin, A.N.Bocharov. Advanced MHD Assisted Mixing of Reacting Streams // In: 39th AIAA Aerospace Sciences Meeting and Exibit, Reno, Nevada, 8-11 January, 2001, AIAA Paper 2001-0793.

13. V.A.Bityurin, A.N.Bocharov, D.S.Baranov. Parametric Study of Advanced Mixing of Fuel/Oxidant System in High Speed Gaseous Flows and Experimental Validation Planning // In: Final Report ISTC No.l809p, August 2001,57 p.

14. A.N.Bocharov, S.B.Leonov, D.S.Baranov, I.B.Klement'eva, V.A. Bityurin, MHD-Enhanced Mixing and Combustion in Co-current Streams // In: Proc. of 4th workshop on magnetoplasma aerodynamics for aerospace applications, Moscow, 9-11 April, 2002, p.220 -230.

15. Т.Н.Абрамович, Т.А.Гиршович, С.Ю.Краменинииков, А.Н.Секундов, И.П.Смирнова. Теория турбулентных струй. М.: Наука, 1984,720 с.

16. Г.Н Абрамович. Прикладная газовая динамика. М.: Наука, 1969,824 с.

17. Т.Н. Абрамович. Турбулентные свободные струи жидкостей и газов. М.: Государственное энергетическое издательство, 1948,288 с.

18. Бай Ши-И. Турбулентное течение жидкостей и газов. М.: Издательство иностранной литературы, 1962,244 с.

19. J.M.Ottino. The kinematics of mixing: stretching, chaos, and transport. Cambridge University Press, 1997,364 p.

20. А.С.Монин, А.М.Яглом. Статистическая гидромеханика, механика турбулентности, часть 1. М.: Наука, 1966, с.505 529.

21. А.С.Монин, А.М.Яглом. Статистическая гидромеханика, механика турбулентности, часть 2. М.: Наука, 1967, с.513-517.

22. W.F.H.Merck et. al. Constricted Discharges in Ar-Cs MHD Generators. Metallurgical Technologies, Energy Conversion and MHD Flows. V.148 Progress in Astronautics and Aeronautics, 1990, p.373 397.

23. Бай Ши-И. Магнитная газодинамика и динамика плазмы. М.: Мир, 1964,304 с.

24. Т.Каулинг. Магнитная гидродинамика. М.: Издательство иностранной литературы, 1959,132 с.

25. В.И.Алфёров. К вопросу определения поля плотности потока при визуализации вихревых жгутов методом высоковольтного разряда. // Труды ЦАГИ, 1972, Вып. 1421, с.13-21.

26. В.И.Алфёров, А.С.Бушмин. Электрический разряд в сверхзвуковом потоке воздуха. // ЖЭТФ, 1963, т.44, No 6, с.1775.

27. В.И.Алфёров, А.С.Бушмин, Б.В.Калачев. Экспериментальное исследование свойств электрического разряда в потоке воздуха. // ЖЭТФ, 1966, т.51, No 5, с.1281.

28. В.И.Алфёров, Л.М.Дмитриев. Электрический разряд в потоке газа при наличии градиентов плотности. // ТВТ, 1985, т.23, No 4 с.677 682.

29. В.И. Алфёров. Исследование структуры электрического разряда большой мощности в высокоскоростном потоке воздуха. // МЖГ, 2004.

30. А.П.Ершов, И.Б.Тимофеев, В.АЛерников, В.М. Шибков. Газовые разряды в сверхзвуковых потоках воздуха. // Прикл. Физика, 1999, No 6, с. 12.

31. А.П.Ершов и др. Поперечные электрические разряды в сверхзвуковых потоках воздуха Механизмы распространения и неустойчивости разряда // ТВТ, 2004, т.42, No 4, с.516-522.

32. С.А.Двинин и др. Моделирование разряда постоянного тока в поперечном сверхзвуковом потоке газа // ТВТ, 2004, т.42, No 2, с.181 —191.

33. А.П.Ершов. Взаимодействие электрических разрядов со сверхзвуковыми газодинамическими возмущениями. // Автореф. дисс. на соискание уч. ст. доктора физ.-мат. наук, М.: МГУ, 2006,36 с.

34. А.П.Ершов, В.А.Черников, В.М. Шибков. Поперечные электрические разряды в сверхзвуковых потоках воздуха М.: МГУ им. М.В.Ломоносова, Физический факультет 2006,94 с.

35. S.Leonov, V.Bityurin, A.Klimov, Yu.Kolesnichenko, A.Yuriev. Influence of Structural Electric Discharges on Parameters of Streamlined Bodies in Airflow. // In: 32nd AIAA Plasmodynamic and Laser Conf., 11-14 June 2001, Anaheim, CA, AIAA Paper 2001-3057.

36. V.G.Brovkin et. al. Study of Microwave Plasma-Body Interaction in Supersonic Airflow. // In: 30th AIAA Plasmodynamic and Laser Conf., June 1999, Norfolk, Virginia, AIAA Paper 1999-3740.

37. T.Cain et al. Influence of a HF Corona Plasma Structure on Drag of an Axial-Symmetric Body in a Supersonic Airflow. // In: 3rd Workshop on WIG, November 1999, Norfolk, Virginia, AIAA Paper 1999-4856.

38. С.Б.Леонов, Д.А.Яранцев. Управление отрывными явлениями в высокоскоростном потоке с помощью приповерхностного электрического разряда. // АМГД, No 4,2006.

39. S.Leonov, A.Kuryachii, D.Yarantsev, A.Yuriev. Supersonic/Transonic Flow Control by Electro-Discharge Plasma Technique. In: Proc.of ICAS-2006, Hamburg, Germany, 4-5 Sept, 2006.

40. S.Leonov, A.Kuryachii, V.Soloviev, D.Yarantsev. Study of Friction and Separation Control by Surface Plasma. In: 42nd Aerospace Sciences Meeting and Exhibit, Reno, NV, 2004, AIAA Paper 2004-0512.

41. S.Leonov et. al. Hydrocarbon Fuel Ignition in Separation Zone of High Speed Duct by Discharge Plasma. // In: Proc.of 4th Workshop PA and MHD in Aerospace Application, Moscow, IVTAN, April 2002.

42. С.Б.Леонов. Воздействие электрических разрядов на структуру и параметры высокоскоростного воздушного потока. // Автореф. дисс. на соискание уч. ст. доктора физ.-мат. наук, М.: ОИВТРАН, 2006,41 с.

43. A.A.Sakharov, N.P.Mende, S.V.Bobashev, D.M.Van Wie. MHD Control of Supersonic Flow about Body. // Pis'ma v Zhurnal Tekhnicheskoi Fiziki, 2006, v.32, No 14, p.40 45.

44. Б.Брон. Движение дуги в магнитном поле. М.: Госэнергоиздат, 1944.

45. Б.Брон. Электрическая дуга в аппаратах управления. М.: Госэнергоиздат, 1954.

46. Г.А.Кукеков. Исследование дуги постоянного тока в магнитном поле. // ЖТФ, 1941, т.11, вып.Ю, с.972 978.

47. М.Ф.Жуков, В.Я.Смоляков, Б.А.Урюков. Электродуговые нагреватели газа М.: Наука, 1973,232 с.

48. Б.А.Урюков. Теория идеальной электрической дуги в коаксиальном плазмотроне. // ПМТФ, 1969, No 1.

49. В.А.Кисель, Б.А.Урюков, В.ИДцров. Обобщение вольтамперных характеристик коаксиального плазмотрона с магнитной стабилизацией электрической дуги. // Изв. СО АН СССР, Серия техн. наук, 1967, No 3, вып.2.

50. L.MShpanin et. al. The Interaction of High Current Electric Arcs with Spatially Varying Cross Magnetic Fields. // In: Proc. Of 15th Int. Conf. On Gas Discharges, Toulouse, 5-10 Sept. 2004, p.183 -186.

51. L.M.Shpanin et. al. The Use of Spatially Distributed Magnetic Fields for Arc Control and Current Interruption at Atmospheric Pressure. In: 16th International Conference on Gas Discharges and their Applications, Xian, China, Sept. 11-155,2006, p.41 44.

52. L.M.Shpanin et. al. Formation and Propulsion of an Atmospheric Pressure Plasma Ring. In: 16th International Conference on Gas Discharges and their Applications, Xian, China, Sept. 11-155,2006, p.5-8.

53. T.Ombrello et. al. Enhancement of Combustion and Flame Stabilization Using Stabilized Non-Equilibrium Plasma. // In: 43 rd Aerospace Sciences Meeting and Exhibit, Reno, NV,2005, AIAA Paper 2005-1194.

54. T.Ombrello et. al. Ignition Enhancement Using Magnetic Gliding Arc. // In: 45th Aerospace Sciences Meeting and Exhibit, Reno, NV, 2007, AIAA Paper 2007-1025.

55. S.O.Macheret, M.N.Shneider, R.B.Miles. MHD and EHD Control of Hypersonic Flows of Weakly Ionized Plasmas. // In: 14th International Conference on MHD Power Generation and High Temperature Technologies, Maui, Hawaii, 2002, AIAA Paper 2002-2249.

56. S.O.Macheret. Ptysics of Magnetically Accelerated Nonequilibrium Surface discharge in High-Speed Flow. // In: 44th Aerospace Sciences Meeting and Exhibit, Reno, NV, 2006, AIAA Paper 2006-1005.

57. A.V.Likhanskii et. al. Modeling of Interaction Between Weakly Ionized Near-Surface Plasmas and Gas Flow. // In: 44th Aerospace Sciences Meeting and Exhibit, Reno, NV,2006, AIAA Paper 2006-1204.

58. S.H.Zaidi et. al. Snowplow Surface Discharge in Magnetic Field for High Speed Boundary Layer Control. // In: 44th Aerospace Sciences Meeting and Exhibit, Reno, NV, 2006, AIAA Paper 2006-1006.

59. Ю.П.Райзер. Физика газового разряда. М.: Наука, 1987,536 с.

60. В.Финкедьнбург, Г.Меккер. Электрические дуги и термическая плазма. М.: Издательство Иностранной Литературы, 1961,370 с.

61. Э.М.Базелян, Ю.П.Райзер. Искровой разряд. М.: Издательство МФТИ, 1997,320 с.

62. Э.М.Базелян, Ю.П.Райзер. Физика молнии и молниезащиты. М.: Физматлит, 2001, 320 с.

63. A.Bocharov, I.Klement'eva, A.Klimov, V.Bityurin. A Study of MHD Assisted Mixing and Combustion in Couter-Flow Streams // In: 43rd Aerospace Sciences Meeting and Exhibit, Reno, NV, 2005, AIAA Paper 2005-0600.

64. I.Klement'eva, A.Bocharov, V.Bityurin, A.Klimov. MHD Assisted Mixing through the Experimental and Numerical Study // In: 27th ICPIG, Eindhoven, the Netherlands, 18-22 July, 2005, CD.

65. A.Bocharov, V.Bityurin, I.Klement'eva, A.Klimov, Numerical and Experimental Study of MHD Assisted Mixing and Combustion // In: 44th Aerospace Sciences Meeting and Exhibit, Reno, NV, 2006, AIAA Paper 2006-1009.

66. V.Bityurin, I.Klement'eva, A.Klimov. Mixing and combustion intensification of high speed reacting flows by the electrical discharges in magnetic field // In: 31st Interna-tional symposium on combustion, Heidelberg, Germany, August 6-11,2006, CD.

67. A.Bocharov, V.Bityurin, I.Klement'eva, A.Klimov. Numerical and Experimental Study of MHD Assisted Mixing and Combustion // In: 45th Aerospace Sciences Meeting and Exhibit, Reno, NV, 2006, AIAA Paper 2007-1024.

68. Klementyeva, A. Bocharov, V. Bityurin. Electrical Discharges Gas Flows Interac-tion in External Magnetic Fields // In: 7th International Workshop on Magnetoplasma Aerodynamics, Moscow, May 24-27,2007, CD.

69. I. Klementyeva, V. Bityurin and A. Bocharov. Constricted Discharge Interaction With High Speed Gas Flows // In: 38th AIAA Plasmadynamics/16lh MHD Energy Conver-sion Conference, Miami, FL, June 25-28,2007, AIAA Paper 2007-4141.

70. М.Я.Выгодский. Справочник по высшей математике. M.: Физматлит, 1963, с.751 833.

71. В.М.Юренев, П.Д.Лебедев. Теплотехнический справочник. Т.2. М.: Энергия, 1976, 396с.

72. Т.М.Бапгга. Гидравлика, гидравлические машины и гидравлические приводы. М.: Машиностроение, 1970,504 с.

73. И.Е.Идельчик. Справочник по гидравлическим сопротивлениям. М.: Машиностроение, 1975,560 с.

74. A.Bocharov, V.Bityurin, I.Klement'eva, and S.Leonov. A Study of MHD Assisted Mixing and Combustion // In: 41st Aerospace Sciences Meeting and Exhibit, Reno, NV, 2003, AIAA Paper 2003-5878.

75. И.Б. Клементьева, A.H. Бочаров, B.A. Битюрин. Особенности взаимодействия электрического разряда с газовым потоком во внешнем магнитном поле // Письма в ЖТФ Т.ЗЗ. Вып.22,2007, с.16-22.

76. Л.М.Милн-Томсон. Теоретическая гидродинамика. М.: Мир, 1964,656 с.

77. Л.Г.Лойцянский. Механика жидкости и газа. М.: Дрофа, 2003,840 с.

78. Г.Биркгоф. Гидродинамика. М.: Издательство иностранной литературы, 1963,244с.

79. ЛДЛандау, В.М.Лифпшц. Гидродинамика. М.: Наука, 1988,736 с.

80. Н.А.Слезкин. Динамика вязкой несжимаемой жидкости. М.: Изд. Технико-теоретической литературы, 1955,520 с.

81. А.Н.Патрашев. Гидромеханика. М.: Военно-морское издательство, 1953,720 с.

82. Дж.Э.Андерсон. Явления переноса в термической плазме. М.: Энергия, 1972.

83. Александров Н.Л., Базелян А.Э., Базелян Э.М., Кочетов И.В., Моделирование длинных стримеров в газе атмосферного давления, Физика плазмы, 1995, Т 21, № 1, с. 60-80.

84. Александров Н.Л., Базелян А.Э., Кончаков А.М., Параметры плазмы в канале длинного лидера в воздухе, Физика плазмы, 2001, Т 27, № 10, с. 928-939.

85. В.С.Имшенник, Н.А.Боброва. Динамика столкновительной плазмы. М.: Энергоатомиздат, 1997.

86. В.А.Рожанский, Л.Д.Цендин, Столкновительный перенос в частично ионизованной плазме. Москва. М.: Энергоатомиздат, 1988.

87. В.И.Артемов, Ю.СЛевитан, ОЛ.Синкевич. Неустойчивости и турбулентность в низкотемпературной плазме. М.: Издательство МЭИ, 1994,412 с.

88. Д.А.Франк-Каменецкий. Лекции по физике плазмы. М.: Атомиздат, 1964,284 с.

89. Н.С.Бахвалов, Н.П.Жидков, Г.М.Кобельков. Численные методы. М.: Наука, 1987.

90. И.С.Березин, Н.П.Жидков. Методы вычислений. М.: Ф.М., 1962,464 с.

91. А.А.Амосов, Ю.А.Дубинский, Н.В.Копченова. Вычислительные методы для инженеров. М.: Высшая школа, 1994,596 с.

92. L.SJacobsen, J.A.Schetz, G.S.Elliot, J.W.Crafton. An experimental investigation of a DC plasma-torch igniter // In: 11th AIAA Conf., Orleans, France, 29 Sept 4 Oct, 2002.

93. И.Б.Клементьева. Интенсификация смешения и горения в спутных струях во внешних электромагнитных полях, физическая и математическая модели дугового разряда. Бакалаврская работа, 2002.

94. И.Б.Клементьева. Разработка физической и численной моделей электрического разряда для исследования процессов интенсификации смешения и зажигания в высокоскоростных потоках. Магистерская диссертация, 2004.

95. A.Bocharov, V.Bityurin, E.Filimonova and A.Klimov, Numerical Study of Plasma Assisted Mixing and Combustion in Non-Premixed Flows. // In: 42nd Aerospace Sciences Meeting and Exhibit, Reno, NV, 2004, AIAA Paper 2004-1017.

96. Andriatis A.V., Zhluktov S.A., Sokolova I.A. Transport Coefficients for Chemical Nonequilibrium Components of Air Mixtures. // J. Mathematical Modeling, v.4, No.l, 1992.