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

005001336

На правах рукописи УДК 537.52

Киндышева Светлана Викторовна

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

01.04.08 - физика плазмы

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

10 НОЯ 2011

Долгопрудный - 2011

005001336

Работа выполнена в Федеральном государственном бюджетном

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

«Московский физико-технический институт (государственный университет)»

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

доктор физико-математических наук, профессор МФТИ Александров Николай Леонидович

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

доктор физико-математических наук, профессор, нач. отдела ТРИНИТИ Напартович Анатолий Петрович

Доктор физико-математических наук, старший научный сотрудник ОИВТ РАН Найдис Георгий Вениаминович

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

Институт общей физики им. А.М. Прохорова РАН

СО

Защита состоится 30 ноября 2011 года в/У часов на заседании диссертационного совета Д 212.156.06 при Московском физико-техническом институте по адресу: 117393, г. Москва, ул. Профсоюзная, 84/32, корпус В-2.

Отзывы направлялись по адресу: 141700, Московская обл., г. Долгопрудный, Институтский пер., д. 9, МФТИ.

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

Автореферат разослан_октября 2011 года.

Ученый секретарь

диссертационного совета Д 212.156.06 кандидат технических наук

$.......Чубинский Н.П.

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

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

Работа посвящена теоретическому исследованию процессов в плазме высоковольтного наносекундного импульсного разряда, который развивается в объеме в форме волны ионизации или вдоль диэлектрической поверхности (поверхностный барьерный разряд). Указанные типы разряда поддерживаются при высоких значениях приведенного электрического поля Е/Ы (IV - концентрация нейтральных частиц), что обеспечивает значительный вклад энергии разряда в ионизацию молекул и возбуждение их электронных уровней. Характерные концентрации электронов в таких разрядах лежат в диапазоне 1012 - 1015 см'3. Кроме этого, для волн ионизации характерно наличие области параметров, в которой развитие разряда происходит однородно, что позволяет значительно упростить моделирование разрядных процессов за счет перехода к нульмерной постановке задачи.

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

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

з

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

Цель работы

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

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

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

- теоретическое изучение механизмов быстрого (на временах менее 1 мкс при 1 атм) нагрева в воздухе и других смесях N2:02 при высоких приведенных электрических полях £Ж(300 - 1000 Тд, где 1 Тд = 10"17 В см2).

Научная новизна работы

1. Выполнено моделирование распада плазмы высоковольтного наносекундного разряда в молекулярных газах при высоких (~1012 см"3)

4

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

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

3. Впервые теоретически исследован быстрый нагрев газа в разрядной плазме в воздухе и других смесях Ы2:02 при аномально высоких (~103 Тд) значениях приведенных электрических полей и определены основные каналы нагрева в этих условиях.

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

1. Результаты моделирования распада плазмы высоковольтного наносекундного разряда в молекулярных газах и их смесях для низких давлений и комнатной температуры, а также для давлений порядка атмосферного и температур газа 600 - 2000 К. Влияния на скорость распада термализации электронов, кластеризации положительных ионов и, тройной электрон-ионной рекомбинации с третьим телом -электроном. Результаты сравнения моделирования с измерениями плотности электронов при распаде плазмы.

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

5

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

3. Кинетическая модель и результаты моделирования быстрого нагрева в разрядной плазме в воздухе и других смесях N2:02 при аномально высоких (~103 Тд) приведенных электрических полях.

Практическая и научная ценность работы

Модели, развитые в данной работе могут быть использованы:

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

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

3. для определения доли энергии разряда, быстро переходящей в тепло, в воздухе и других смесях Ы2:02 при наличии больших приведенных электрических полей.

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

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

Список работ по теме диссертации содержит 35 публикаций, в том числе - 9 статей в реферируемых научных журналах, входящих в список ВАК РФ. Личный вклад автора в работы, вошедших в диссертацию, является

6

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

Основные результаты работы докладывались на конференциях: 48 - 52 Научные конференция МФТИ (Долгопрудный, 2005-2009), 45th AIAA Aerospace Sciences Meeting and Exhibit (Reno, Nevada, USA, 2007), 7th International Workshop on Magnetoplasma Aerodynamics (Moscow, Russia, 2007), 26th International Symposium on Shock Waves (Gôttingen, Germany, 2007), 28th ICPIG (Prague, Czech Republic, 2007), 21st ICDERS (Poitiers, France, 2007), Workshop on Frontiers in Low Temperature Plasma Diagnostics (Beverley, UK, 2007), 2nd European Conference for AeroSpace Sciences (Université Libre de Bruxelles, 2007), AEROSPACE THEMATIC WORKSHOP on the "Fundamentals of Aerodynamic Flow and Combustion control by Plasmas" (Villa Monastero, Varenna, Como lake, Italy, 2007), 32nd International Symposium on Combustion (Canada, 2008), XVII International Conference on Gas Discharges and Their Applications (Cardiff, UK, 2008), ISPC18 (Kyoto, Japan, 2008), 46th AIAA Aerospace Sciences Meeting and Exhibit (2008, Reno, Nevada, USA), XXXVI - XXXVIII Международная конференция по физике плазмы и управляемому термоядерному синтезу (Звенигород, 2009 - 2011), 19th International Symposium on Plasma Chemistry (Bochum, 2009), 8th International Workshop on Magneto-Plasma Aerodynamics (Moscow, 2009), 22nd International Colloquium on the Dynamics of Explosions and Reactive Systems (Minsk, Belarus, 2009), 47th AIAA Aerospace Sciences Meeting and Exhibit (Orlando, Florida, 2009), 10th International Workshop on Magneto-Plasma Aerodynamics (Moscow, 2011).

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

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

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

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

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

Вторая глава посвящена изучению распада плазмы, созданной высоковольтным наносекундным импульсным разрядом в молекулярных газах. Этот вопрос исследовался экспериментально при измерении плотности электронов СВЧ интерферометром в послесвечении разряда как при комнатной температуре [1], так и в условиях, когда начальная температура газа была повышенной и лежала в диапазоне 600 - 2400 К при нагреве газа за фронтом ударной волны [2]. Цель данной работы состоит в моделировании распада плазмы для условий этих экспериментов.

В первой части данной главы исследуется распад плазмы в пяти различных газах: чистых Ы2, 02, С02, воздухе и парах воды при комнатной температуре газа и давлениях 1-10 Тор.

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

где пе0 = псЦ=0) и а - коэффициент электрон-ионной рекомбинации. В более сложном случае это уравнение уже не справедливо, но с его помощью

8

можно определить эффективный коэффициент рекомбинации а.ф = М(пе0т), где т - время, за которое начальная плотность электронов уменьшается в 2 раза. Величина ае/ в реальной ситуации может зависеть не только от температуры электронов Те, но и от давления газа и начальной концентрации электронов.

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

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

\

...... . -п

¡¡т. ¿г.

с1к, Ык, п,п{ +-п,

где Т - температура газа, \'Е - частота релаксации энергии электронов при столкновении с нейтральными частицами, п, — плотность ионов и к2 и к3 ~ константы скорости парной и тройной (с третьим телом — электроном) электрон-ионной рекомбинации.

Сравнение расчетного коэффициента рекомбинации ае/ с экспериментальными данными в N2, С02 и парах воды показывает, что «е/

зависит от времени во всех исследуемых газах: за счёт влияния термализации энергии электронов (при низких давлениях в N2), а также за счет образования комплексных положительных ионов в трехтельных ионно-молекулярных процессах (в парах воды, Ы2 и С02). Константа скорости рекомбинации с комплексными ионами на порядок выше, чем с простыми молекулярными ионами. В отличие от N2 и С02, в парах воды «е/значительно возрастает из-за существования нескольких типов устойчивых комплексных ионов. Сравнение результатов моделирования динамики концентрации электронов с данными эксперимента [1] для Ы2 и С02 приведены на рис. 1.

а) б)

Рис. 1. Зависимость концентрации электронов от времени: а) в N1 для р=3 Тор и б) в С02 для р=3.5 Тор. Сплошные кривые соответствуют измерениям, штриховые - вычислениям в предположении мгновенной термализации электронов, а цприхиунктирые -вычислениям с учетом термшшзации электронов за конечное время.

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

пренебречь. При моделировании все константы скорости брались из

ю

экспериментов, за исключением константы скорости тройной рекомбинации с третьим телом - электроном. Константа скорости этого процесса экспериментально и теоретически исследовалась только для атомарных ионов. Именно эта константа обычно и приводится в обзорах по скоростям рекомбинации и используется при моделировании рекомбинации и атомарных, и молекулярных ионов. В то же время имеются расчеты [3] для модельной системы, где показано, что константа тройной электрон-ионной рекомбинации к3 с третьим телом электроном, к3 = аТе~", может существенно (на порядок величины и более) увеличиваться при переходе от атомарных ионов к молекулярным. Показатель степени, равный п = 4.5 в случае атомарных ионов, может также уменьшиться до п ~ 2 - 2.5 в случае молекулярных ионов. Эти различия между рекомбинацией атомарных и молекулярных ионов связаны с возможностью образования при рекомбинации молекулярных ионов сильновозбужденной нейтральной молекулы, способной разваливаться за счет предиссоциации [3]. Учет этого эффекта позволил достичь согласия между расчетом и измерениями и в случае распада плазмы в воздухе и 02. На рис. 2 показано сравнение экспериментальной и вычисленной динамики электронов во времени в этих газах. Для обнаружения указанного эффекта важно, что в рассматриваемых условиях во время распада плазмы доминирующим ионом является простой ион О/, скорость диссоциативной рекомбинации для которого существенно меньше аналогичной величины для комплексных ионов. При распаде же плазмы в N2, С02 и парах воды, как показывают дополнительные расчеты, ионы являются комплексными, и роль тройной рекомбинации относительно мала.

Во второй половине главы исследуется распад плазмы, созданной высоковольтным наносекундным разрядом при повышенных температурах (600 - 2400 К) в воздухе и смесях Ы2:02:С02 в диапазоне начальных давлений 0.05- 1.2 атм.

и

а) б)

Рис. 2. Зависимость концентрации электронов от времени в послесвечении разряда а) в кислороде при давлении 8 Тор и б) в воздухе при давлении 2.5 Тор. Кривая 1 - измерения, кривая 2 - расчет с коэффициентом тройной рекомбинации кг = аТе'" из [4], 3 - расчет с увеличенным значением а в 10 раз, 4 - расчет с увеличенным значением а в 10 раз и и = 2.

Сравнение эффективного коэффициента рекомбинации, вычисленного для воздуха и смесей, содержащих углекислый газ, с экспериментальными измерениями [2] показывает, что наблюдается хорошее согласие расчетов с экспериментальными данными при Т>1300 К, но при более низких температурах вычисленные коэффициенты рекомбинации значительно превышают экспериментальные значения. Расчёты для воздуха приведены на рис. 3.

В данных условиях главными каналами гибели электронов являются диссоциативная рекомбинация с 0{ и тройное прилипание электронов к молекулам 02. Обнаруженное различие, как показывают дополнительные исследования, связаны с наработкой атомов О во время высоковольтного наносекундного разряда. В результате процессы прилипания компенсируются отлипанием электронов от ионов 02~ в столкновениях с атомами О, и гибель электронов, в соответствии с экспериментом, определяется только рекомбинацией с ионом О/. В отсутствии атомов О распад плазмы может носить рекомбинационный характер с эффективным 1

коэффициентом рекомбинации, не зависящим от скоростей прилипания и отлипания электронов.

0.09 атм ▲ 0.13 атм

* 0.4 атм я 0.165 атм • 0.2 атм

♦ 0.22 атм * 0.39 атм

4 \ ▼ 0.4 атм в 0.616 атм

о.К 1.04 агм * 1.17 а тм

я

•

500 1000 1500 2000 2500

т,к

Рис. 3. Зависимость эффективного коэффициента электрон-ионной рекомбинации (символы), измеренного в воздухе, от температуры газа. Сплошная кривая соответствует коэффициенту диссоциативной рекомбинации для иона Ог , измеренному в работе [5], а штриховые - вычислениям при различных давлениях газа в'отсутствии атомов О.)

В третьей главе (в первой части) исследуется наработка активных частиц под действием высоковольтного наносекундного импульсного разряда в стехиометрических смесях насыщенных углеводородов (от метана до пентана) с молекулярным кислородом, на 90% разбавленных аргоном, при давлении порядка атмосферного и температурах в диапазоне 800 - 3000 К. Под активными частицами здесь понимаются электронно-возбужденные молекулы 02 и атомы Аг, атомы О, Н, радикалы углеводородных молекул, а также электроны и положительные ионы.

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

13

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

параметров [6 - 8].

Во второй части исследуется влияние добавки азота на наработку активных частиц и воспламенение под действием высоковольтного наносекундного разряда смеси СН4:02^2:Аг = 1:4:15:80 для диапазона давлений 0.4 - 2.0 атм и температур 1200 - 1800 К. Вычисления показывают, что добавление с одной стороны приводит к заметному уменьшению средней энергии электронов из-за эффективного возбуждения колебательных уровней азота, а с другой - за счет тушения электронно-возбужденных состояний Ы2 на Ог появляется эффективный канал наработки атомов О, которые и в этом случае остаются доминирующим сортом активных частиц.

Далее на основе данных по концентрациям активных частиц, наработанных в разряде, моделируется изменение во времени концентраций нейтральных частиц и газовой температуры на фазе воспламенения на основе численного решения уравнений баланса этих частиц и уравнения, описывающего изменение температуры за счет химических реакций. Как видно из рис. 4, полученные расчетные времена индукции воспламенения хорошо согласуются с измерениями [9] при низких (-0.5 атм) давлениях газа. Различие же между расчетными и измеренными временами индукции воспламенения при высоких (-2 атм) давлениях может быть обусловлено влиянием неоднородности развития разряда и воспламенения, которые были зарегистрированы в эксперименте [9].

Рис. 4. Время индукции самовоспламенения (сплошные кружки - эксперимент, полые -расчет) и воспламенения после разряда (треугольники и звездочки - эксперимент, вертикальные отрезки - расчет) в смеси CH4:02:N2:Ar = 1:4:15:80 в зависимости от обратной газовой температуры при различных давлениях газа. Сплошная линия -аппроксимация измеренного времени индукции самовоспламенения, штриховая кривая -расчет времени индукции самовоспламенения при давлении 0.5 атм. Звездочки соответствуют эксперименту с одновременным измерением временных зависимостей разрядного тока и электрического поля. Для этих условий и проводились расчеты.

В четвертой главе исследуется быстрый нагрев воздуха и других смесей N2:02 при комнатной температуре в присутствии высоких приведенных электрических полей. За основу взята модель [10] быстрого нагрева в плазме воздуха при умеренных (-100 Тд) приведенных электрических полях FJN и дополнена процессами, важными при высоких значениях E/N. Экспериментальное исследование эффективности быстрого нагрева газов в плазме наносекундных разрядов в воздухе можно найти в работе [11] для 20 Topp и в [12, 13] для 1 атм. Их анализ был выполнен в работе [14], где было показано, что доля энергии, быстро передаваемая от электронов в тепло, может составлять 36 - 40 % при E/N- 600 Тд в случае малых давлений (р = 20 Topp) и 56 - 66 % при E/N ~ 900 Тд в случае больших давлений (р = 1 атм). Развитые ранее кинетические модели относились к меньшим значениям E/N и не могли объяснить эти данные.

а)

б)

Рис. 5. Изменение во времени доли электронной энергии, быстро передаваемой в тепло в сухом воздухе а) при 20 Topp и EIN - 600 Тд и б) при 1 атм и EIN = 10 Тд. Расчеты выполнены при различных значениях ле/. Вертикальные отрезки соответствуют экспериментальным данным [14].

На рис. 5 показана расчетная динамика быстрого нагрева в воздухе для условий эксперимента. Расчеты, дающие оценку эффекта сверху, показывают, что электронная энергия, быстро переходящая в тепло, растет от 10% при EIN = 102 Тд до 30 - 55% при EIN = 103 Тд. Результат существенно зависит от давления газа и слабо зависит от концентрации электронов nef в момент окончания разряда. Этот рост с давлением обусловлен влиянием тройного прилипания электронов к 02 и тройной ион-ионной рекомбинации.

Известно, что при умеренных значениях E/N эффективность быстрого нагрева в смесях N2:02 снижается с уменьшением доли 02 [10]. В настоящей работе рассмотрена эффективность быстрого нагрева в различных смесях N2:02 при высоких значениях E/N. Здесь доля энергии, идущая в быстрый нагрев, также меньше в смесях с меньшим содержанием кислорода.

Сравнение результатов настоящего расчёта с расчетом [15] и с расчетом [10] даёт хорошее согласие при EIN = 100 - 200 Тд. При более высоких значениях EIN результаты данной работы лежат заметно выше, чем результаты расчетов, в которых учитываются только каналы быстрого нагрева, доминирующие при умеренных приведенных электрических полях.

Это различие связано с учетом каналов нагрева через сильновозбужденные состояния N2 и каналов, связанных с заряженными частицами.

Расчёты для более высоких полей качественно согласуются с результатами анализа наблюдений поверхностного наносекундного высоковольтного импульсного разряда в воздухе при 100 Topp [16], где было получено, что в быстрый нагрев воздуха уходит 40±10% вложенной в разряд энергии при EIN в диапазоне 100 - 500 Тд.

После выполнения данной работы и опубликования ее результатов появилась работа [17], в которой также теоретически изучается быстрый нагрев в воздухе и других смесях N2:02 при высоких приведенных электрических полях, но при этом используется несколько другая кинетическая схема для описания быстрого перехода энергии свободных электронов в тепло. Одним из главных отличий этой модели является предположение о том, что при диссоциативной рекомбинации электронов с ионами 04+ и при тройной ион-ионной рекомбинации кислородных ионов одна из молекул 02 диссоциирует. Следовательно, в этих процессах энергия, переходящая быстро в тепло, должна быть меньше той, что предполагается в данной работе. В расчетах [17] доля этой энергии не превышает 40 %, что несколько меньше результатов данной работы.

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ ДИССЕРТАЦИИ:

1. Выполнено численное моделирование распада плазмы высоковольтного наносекундного разряда, инициированного в N2, С02 и парах Н20 при давлениях газа 1-10 Topp и комнатной температуре. Получено хорошее согласие с измерениями концентрации электронов в зависимости от времени при распаде плазмы. Показано, что наблюдаемая зависимость эффективного коэффициента рекомбинации от давления и времени связана с динамикой остывания электронов (для N2), образованием кластерных ионов в тройных столкновениях (для всех исследованных

17

газов) и тройной электрон-ионной рекомбинацией с третьим телом -молекулой (для Н20).

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

Теоретически изучен распад плазмы высоковольтного наносекундного разряда в воздухе и смесях Ы2:02:С02 при повышенных газовых температурах (500 - 2500 К). Показано, что при наличии быстрого прилипания и отлипания электронов распад плазмы в изучаемых смесях может описываться рекомбинационным законом с независящим от времени эффективным коэффициентом рекомбинации, который значительно больше (до порядка величины и более) коэффициента диссоциативной рекомбинации. При этом, в отличие от ранее изученных случаев, эффективный коэффициент рекомбинации не зависит от скоростей прилипания и отлипания электронов. Указанный режим распада плазмы не был реализован в конкретных экспериментах из-за высокой плотности атомов О, наработанных во время высоковольтного наносекундного разряда. Учет разрушения отрицательных ионов на атомах О позволил объяснить полученную в эксперименте динамику распада плазмы в рассматриваемых условиях.

Развита кинетическая схема и выполнено численное моделирование наработки активных частиц в плазме высоковольтного наносекундного разряда в стехиометрических смесях углеводородов СкН2ы2 (к = 1 - 5) с кислородом или воздухом при их разбавлении аргоном применительно к конкретным ударноволновым экспериментам, в которых измерялись электрическое поле в разряде и его ток в зависимости от времени. Показано, что среди наработанных разрядом активных частиц преобладают атомы О. Добавление N2 к газовой смеси с одной стороны приводит к уменьшению эффективности образования активных частиц в разряде из-за снижения средней энергии электронов при заданном электрическом поле, а с другой - к появлению новых эффективных каналов образования атомов О через тушение электронно-возбужденных молекул N2 на молекулах 02.

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

Изучен быстрый (на временах менее I мкс при 1 атм) нагрев в воздухе и других азот-кислородных смесях под действием наносекундного разряда при аномально высоких (до 103 Тд) приведенных электрических полях Е/И, где N - плотность нейтральных частиц. Показано, что при

высоких значениях E/N становится важным нагрев через высоколежащие электронно-возбужденные состояния молекулярного азота и при электрон-ионной и ион-ионной рекомбинации. Выполнена оценка сверху доли разрядной энергии, быстро переходящей в тепло, в зависимости от EJN и давления газа. Показано, что при высоких значениях E/N и атмосферном давлении до 50% энергии разряда может быстро переходить в тепло, что находится в согласии с результатами экспериментов.

Список литературы

[1] Aleksandrov N.L., Kindysheva S.V., Kirpichnikov А.А., Kosarev I.N., Starikovskaia S.M., Starikovskii A.Yu. Plasma decay in N2, C02 and H20 excited by high-voltage nanosecond discharge // J. Phys. D: Appl. Phys. 2007. Vol. 40. No. 15. P. 4493 - 4502.

[2] Aleksandrov N.L., Kindysheva S.V., Kosarev I.N., Starikovskii A.Yu. Plasma decay in air and N2:02:C02 mixtures at elevated gas temperatures // J. Phys. D: Appl. Phys. 2008. Vol. 41. No. 21. 215207.

[3] Collins C.B. Collisional - dissociative recombination of electrons with molecular ions // Phys. Rev. 1965. Vol. 140. P. A1850 - A1857.

[4] Kossyi I.A., Kostinsky A.Yu., Matveev A.A., Silakov V.P. Kinetic scheme of the non-equilibrium discharge in nitrogen-oxygen mixtures // Plasma Sourse Sci. Technol. 1992. Vol. 1. No. 3. P. 207-220.

[5] Cunningham A.J., Hobson R.M. Dissociative recombination at elevated temperatures III. 02+dominated afterglows // J. Phys. B: Atom. Molec. Phys. 1972. Vol. 5. No. 12. P. 2320 - 2327.

[6] Kosarev I.N., Aleksandrov N.L., Kindysheva S.V., Starikovskaia S.M., Starikovskii A.Yu. Kinetic mechanism of plasma-assisted ignition of hydrocarbons Hi. Phys.D: Appl. Phys. 2008. Vol. 41. No. 3.032002.

[7] Kosarev I.N., Aleksandrov N.L., Kindysheva S.V., Starikovskaia S.M., Starikovskii A.Yu. Kinetics of ignition of saturated hydrocarbons by

20

nonequilibrium plasma: C2H6 - to C5Hi2 - containing mixtures // Combustion and Flame. 2009. Vol. 156. No. 1. P. 221 -233.

[8] Косарев И.Н. Кинетика неравновесного воспламенения газовых смесей под воздействием лазерного флэш-фотолиза и наносекундного разряда : диссертация на соискание ученой степени канд. физ.-мат. наук. / Московский физико-технический институт. Москва, 2008.154 с.

[9] Starikovskaia S.M., Kukaev E.N., Kuksin A.Yu., Nudnova M.M., Starikovskii A.Yu. Analysis of the spatial uniformity of the combustion of a gaseous mixture initiated by a nanosecond discharge // Combust. Flame. 2004. Vol. 139 . No. 3. P. 177- 187.

[10] Попов H.A. Исследование механизма быстрого нагрева аота и воздуха в газовых разрядах // Физика плазмы. 2001. Т. 27. Вып. 10. С. 940 - 949.

[И] Anokhin Е.М., Starikovskaia S.M., Starikovskii A.Yu. Energy Transfer in Hypersonic Plasma Flow and Flow Structure Control by Low Temperature Nonequilibrium Plasma // 42nd AIAA Aerospace Sciences Meeting and Exhibit (Reno, NV, 5-8 January 2004) paper AIAA -2004 - 674 .

[12] Roupassov D.V., Nikipelov A.A., Nudnova M.M., Starikovskii A.Yu. AIAA J. 2009. Vol.47. P. 168.

[13] Starikovskii A.Yu., Nikipelov A.A., Nudnova M.M., Roupassov D.V. SDBD plasma actuator with nanosecond pulse-periodic discharge // Plasma Sources Sci. Technol. 2009. Vol. 18. No. 3.034015.

[14] Aleksandrov N.L., Kindysheva S.V., Nudnova M.M., Starikovskiy A.Yu. Mechanism of ultra-fast heating in a non-equilibrium weakly ionized air discharge plasma in high electric fields // J. Phys. D: Appl. Phys. 2010. Vol. 43. No. 25. 255201.

[15] Flitti A., Pancheshnyi S. Gas heating in fast pulsed discharges in N2-02 mixtures // Eur. Phys. J. Appl. Phys. 2009. Vol. 45.21001.

[16] Знаменская И.А., Латфулин Д.Ф., Лутский A.E., Мурсенкова И.В., Сысоев Н.Н. Развитие газодинамических возмущений из зоны

распределенного поверхностного скользящего разряда // ЖТФ. 2007. Т. 77. Вып. 5. С. 10.

[17] Popov N.A. Fast gas heating in a nitrogen - oxygen discharge plasma: 1. Kinenic mechanism // J. Phys. D: Appl. Phys. 2011. Vol. 44.285201.

Список основных публикаций по тёме диссертации

1. Aleksandrov N.L., Kindysheva S.V., Kirpichnikov А.А., Kosarev I.N., Starikovskaia S.M., Starikovskii A.Yu. Plasma decay in N2, C02 and H20 excited by high-voltage nanosecond discharge // J. Phys. D: Appl. Phys.

2007. Vol. 40. No. 15. P. 4493 - 4502.

2. Aleksandrov N.L., Kindysheva S.V., Kosarev I.N., Starikovskii A.Yu. Plasma decay in air and N2:02:C02 mixtures at elevated gas temperatures // J. Phys. D: Appl. Phys. 2008. Vol. 41. No. 21.215207.

3. Kosarev I.N., Aleksandrov N.L., Kindysheva S.V., Starikovskaia S.M., Starikovskii A.Yu. Kinetic mechanism of plasma-assisted ignition of hydrocarbons // J. Phys.D: Appl. Phys. 2008. Vol. 41. No. 3.032002.

4. Kosarev I.N., Aleksandrov N.L., Kindysheva S.V., Starikovskaia S.M., Starikovskii A.Yu. Kinetics of ignition of saturated hydrocarbons by nonequilibrium plasma: CH4-containing mixtures // Combustion and Flame.

2008. Vol. 154. No. 3. P. 569 - 586.

5. Kosarev I.N., Aleksandrov N.L., Kindysheva S.V., Starikovskaia S.M., Starikovskii A.Yu. Kinetics of ignition of saturated hydrocarbons by nonequilibrium plasma: C2H6 - to C5H12 - containing mixtures // Combustion and Flame. 2009. Vol. 156. No. l.P. 221-233.

6. Александров H.JI., Киндышева C.B., Кунаев E.H., Стариковская С.М., Стариковский А.Ю. Моделирование динамики воспламенения метано-воздушной смеси высоковольтным наносекундным разрядом // Физика плазмы. 2009. Т. 35. Вып. 10. С. 941 - 956.

7. Стариковская С.М., Александров H.JL, Косарев И.Н., Киндышева С.В., Стариковский А.Ю. Воспламенение с помощью низкотемпературной плазмы: кинетический механизм и экспериментальная проверка// Химия высоких энергий. 2009. Т. 43. Вып. 3. С. 259 - 265.

8. Aleksandrov N.L., Kindysheva S.V., Kosarev I.N., Strarikovskaia S.M., Starikovskii A.Yu. Mechanism of ignition by non-equilibrium plasma // Proceedings of the Combustion Institute. 2009. Vol. 32. P. 205 - 212.

9. Aleksandrov N.L., Kindysheva S.V., Nudnova M.M., Starikovskiy A.Yu. Mechanism of ultra-fast heating in a nonequilibrium weakly-ionized air discharge plasma on high electric fields // J. Phys. D: Appl. Phys. 2010. Vol. 43.255201.

Также содержание диссертации опубликовано в тезисах и докладах на 26-ти

отечественных и международных научных конференциях.

Киндышева Светлана Викторовна МОДЕЛИРОВАНИЕ КИНЕТИЧЕСКИЗ ПРОЦЕССОВ В ПЛАЗМЕ ВЫСОКОВОЛЬТНОГО НАНОСЕКУНДНОГО ИМПУЛЬСНОГО РАЗРЯДА В МОЛЕКУЛЯРНЫХ ГАЗАХ Автореферат

Подписано в печать 18.10.2011.Формат 60x84 усл.печ.л.1.0 Тираж 80 экз Заказ № 91 Федеральное государственное бюджетное образовательное учреждение Высшего профессионального образования «Московский физико-технический институт (государственный университет)» Отдел оперативной полиграфии «Физтех-полиграф» 141700,Московская обл.г. Долгопрудный,Институтский пер.9

ВВЕДЕНИЕ.

ГЛАВА 1. ОБЗОР ЛИТЕРАТУРЫ.

1.1. Кинетика распада неравновесной газоразрядной плазмы.

1.2. Моделирование стимулированного плазмой воспламенения газовых смесей.

1.3. Быстрый нагрев в плазме молекулярных газов при высоких электрических полях.

ГЛАВА 2. РАСПАД ПЛАЗМЫ ВЫСОКОВОЛЬТНОГО НАНОСЕКУНДНОГО РАЗРЯДА.

2.1. Введение.

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

2.2.1. Анализ экспериментальных данных.

2.2.2. Кинетическая модель.

2.2.3. Результаты вычислений.

2.2.4. Неоднородность высоковольтного наносекундного разряда.

2.2.5. Модификация кинетической схемы.

2.3. Распад плазмы высоковольтного наносекундного импульсного разряда при высоких газовых температурах.

2.3.1. Экспериментальные данные.

2.3.2. Кинетическая модель.

2.3.3. Результаты вычислений.

2.4. Выводы главы 2.

ГЛАВА 3. КИНЕТИКА НАРАБОТКИ АКТИВНЫХ ЧАСТИЦ И

ВОСПЛАМЕНЕНИЯ ПРЕДЕЛЬНЫХ УГЛЕВОДОРОДОВ СН4-С5Н12 ПОД

ДЕЙСТВИЕМ ВЫСОКОВОЛЬТНОГО НАНОСЕКУНДНОГО РАЗРЯДА.

3.1. Введение.

3.2. Смеси с кислородом: СпН2П+2:02: Аг.

3.2.1. Экспериментальные данные.

3.2.2. Кинетическая модель плазмы в разряде и его послесвечении.

3.2.3. Результаты вычислений для стадии разряда и его послесвечения.

3.3. Смесь с воздухом: СН4:воздух:Аг.

3.3.1. Экспериментальные данные.

3.3.2. Кинетическая модель наработки активных частиц в разряде и процессов воспламенения.

3.3.3. Результаты моделирования разрядных процессов.

3.3.4. Результаты моделирования воспламенения смеси.

3.4. Выводы главы 3.

ГЛАВА 4. БЫСТРЫЙ НАГРЕВ ГАЗА В СЛАБОИОНИЗОВАННОЙ НЕРАВНОВЕСНОЙ ПЛАЗМЕ В СИЛЬНОМ ЭЛЕКТРИЧЕСКОМ ПОЛЕ.

4.1. Введение.

4.2. Расчетная модель.

4.3. Результаты расчета.

4.4. Сравнение с экспериментом и другими расчетами.

4.5. Выводы главы 4.

ВЫВОДЫ И ОСНОВНЫЕ РЕЗУЛЬТАТЫ.

Актуальное гь темы. Работа посвящена теоретическому исследованию процессов в плазме высоковольтного наносекундного импульсного разряда, который развивался в объеме в форме волны ионизации или вдоль диэлектрической поверхности (поверхностный барьерный разряд). Указанные типы разряда поддерживаются при высоких значениях приведенного электрического поля Е/Ы (И — плотность нейтральных частиц), что обеспечивает значительный вклад энергии разряда в ионизацию молекул и возбуждение их электронных уровней. Характерные начальные концентрации электронов лежат в диапазоне Ю12- 1015 см"3. Кроме этого, для волн ионизации характерно наличие области параметров, в которой развитие разряда происходит однородно, что позволяет значительно упросшть моделирование разрядных процессов за счет перехода к нульмерной постановке задачи.

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

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

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

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

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

- теоретическое изучение механизмов быстрого (на временах менее 1 мкс при 1 атм) нагрева в воздухе и других смесях N2:02 при высоких приведенных электрических полях £/¥(300 - 1000 Тд, где 1 Тд = Ю-17 В см2).

Научная новизна. 1. Выполнено моделирование распада плазмы высоковольтного 3 наносекундного разряда в молекулярных газах при высоких (~10~см ) начальных концентрациях электронов, а также произведена верификация модели на основе экспериментальных данных. Впервые моделирование выполнено в предположении о том, что электрон - ионная рекомбинация с третьим телом — электроном для молекулярных ионов может идти гораздо быстрее, чем для атомарных.

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

3. Впервые теоретически исследован быстрый нагрев газа в разрядной плазме в воздухе и других смесях N2:02 при аномально высоких (~103Тд) значениях приведенных электрических полей, и определены основные каналы нагрева в этих условиях.

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

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

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

3. Для определения доли энергии разряда, быстро переходящей в тепло, в воздухе и других смесях N2:02 при наличии больших приведенных электрических полей.

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

Апробация работы. Основные результаты работы докладывались и обсуждались на 48 — й научной конференции МФТИ (Долгопрудный, 2005), 49 — й научной конференции МФТИ (Долгопрудный, 2006), 50-й научной конференции МФТИ (Долгопрудный, 2007), 45th AIAA Aerospace Sciences Meeting and Exhibit (Reno, Nevada, USA, 2007), 7th International Workshopv on Magnetoplasma Aerodynamics (Moscow, Russia, 2007), 26th International Symposium on Shock Waves (Gottingen, Germany, 2007), 28th ICPIG (Prague, Czech Republic, 2007), 21st ICDERS (Poitiers, France, 2007), Workshop on Frontiers in Low Temperature Plasma Diagnostics (Beverley, UK , 2007), 2nd European Conference for AeroSpace Sciences (Université Libre de Bruxelles, 2007), AEROSPACE THEMATIC WORKSHOP on the "Fundamentals of Aerodynamic Flow and Combustion control by Plasmas" (Villa Monastero, Varenna, Como lake, Italy, 2007), 51 — й научной конференции МФТИ (Долгопрудный, 2008), 32nd International Symposium on Combustion (Canada, 2008), XVII International Conference on Gas Discharges-and Their Applications (Cardiff, UK, 2008), 18';' International Symposium on Plasma Chemistry (Kyoto, Japan, 2008), 46th AIAA Aerospace Sciences Meeting and Exhibit (2008, Reno, Nevada, USA), 52-й научной конференции МФТИ (Долгопрудный, 2009), XXXVI Международной конференции по физике плазмы и у управляемому термоядерному синтезу (Звенигород, 2009), 19th International Symposium on Plasma Chemistry (Bochum, 2009), 8th International Workshop on Magneto - Plasma Aerodynamics (Moscow, 2009), 22nd International Colloquium on the Dynamics of Explosions and Reactive Systems (Minsk, Belarus, 2009), 47th AIAA Aerospace Sciences Meeting and Exhibit (Orlando, Florida, 2009), XXXVII Международной конференции по физике плазмы и управляемому термоядерному синтезу (Звенигород, 2010), XXXVIII Международной 6 конференции по фи зике плазмы и управляемому термоядерном у синтезу (Звенигород, 2011), 10th International Workshop on Magneto - Plasma Aerodynamics (Moscow, 2011).

Публикации. Список работ по теме диссертации содержит 35 публикаций, в том числе -9 статей в реферируемых научных журналах, входящих в список ВАК РФ [1—9]. Личный вклад автора в работы, вошедшие в диссертацию, является определяющим (выполнение моделирования, изучение влияния различных параметров на получаемые результаты, анализ результатов и сравнение с экспериментом).

Структура и объем диссертации. Диссертация состоит из введения, четырех глав и заключения. Работа изложена на 115 страницах, включая 43 рисунка, 7 таблиц и список литературы, насчитывающей 114 наименований. Во введение приводится актуальность выбранной темы, научная новизна, а также цели и краткое содержание данной работы.

ВЫВОДЫ И ОСНОВНЫЕ РЕЗУЛЬТАТЫ

1. Выполнено численное моделирование распада плазмы высоковольтного наносекундного разряда, инициированного в N2, С02 и парах Н2О при давлениях газа 1-10 Тор и комнатной температуре. Получено хорошее согласие с измерениями концентрации электронов в зависимости от времени при распаде плазмы. Показано, что наблюдаемая зависимость эффективного коэффициента рекомбинации от давления и времени связана с динамикой остывания электронов (для N2), образованием кластерных ионов в тройных столкновениях (для всех исследованных газов) и тройной электрон-ионной рекомбинацией с третьим телом — молекулой (для Н2О).

2. Моделирование распада плазмы в воздухе и кислороде при аналогичных условиях дало существенно меньшие скорости распада плазмы по сравнению с наблюдаемыми в эксперименте. Высказана гипотеза о том, что в этом случае существенную роль играет тройная рекомбинация электронов с ионами 02+ (третье тело — электрон), и что'константа скорости этого процесса на порядок величины больше аналогичной величины для атомарных ионов, а зависимость скорости рекомбинации от температуры электронов более слабая. Эти предположения, согласующиеся с модельным расчетом сорокалетней давности, позволил» получить согласие между измерениями концентрации электронов и расчетом в случае распада плазмы в воздухе и кислороде.

3. Теоретически изучен распад плазмы высоковольтного наносекундного разряда в воздухе и смесях №:02:С02 при повышенных газовых температурах (500 — 2500 К). Показано, что при наличии быстрого прилипания и отлипания электронов распад плазмы в изучаемых смесях может описываться рекомбинационным законом с независящим от времени эффективным коэффициентом рекомбинации, который значительно больше (до порядка величины и более) коэффициента диссоциативной рекомбинации. При этом, в отличие от ранее изученных случаев, эффективный коэффициент рекомбинации не зависит от скоростей прилипания и отлипания электронов. Указанный режим распада плазмы не был реализован в конкретных экспериментах из—за высокой плотности атомов О, наработанных во время высоковольтного наносекундного разряда. Учет разрушения отрицательных иоиов на атомах О позволил объяснить полученную в эксперименте динамику распада плазмы в рассматриваемых условиях.

4. Развита кинетическая схема и выполнено численное моделирование наработки активных частиц в плазме высоковольтного наносекундного разряда в стехиометрических смесях углеводородов С^Нгк+г (к = 1 - 5) с кислородом или воздухом при их разбавлении

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

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

6. Изучен быстрый (на временах менее 1 мкс при 1 атм) нагрев в воздухе и других азот-кислородных смесях под действием наносекундного разряда при аномально высоких (до 10 Тд) приведенных электрических полях ЕМ, где N - плотность нейтральных частиц. Показано, что при высоких значениях Е/Ы становится важным нагрев через высоколежащие электронно-возбужденные состояния молекулярного азота и при электрон-ионной и ион-ионной рекомбинации. Выполнена оценка сверху доли разрядной энергии, быстро переходящей в тепло, в зависимости от Е/Ы и давления газа. Показано, что при высоких значениях Е/.N и атмосферном давлении до 50 % энергии разряда может быстро переходить в тепло, что находится в согласии с результатами экспериментов.

1. Aleksandrov N.L., Kindysheva S.V., Kirpichnikov A.A., Kosarev I.N., Starikovskaia S.M., Starikovskii A.Yu. // J. Phys. D: Appl. Phys. 2007. Vol. 40. No. 15. P. 4493 4502.

2. Aleksandrov N.L., Kindysheva S.V., Kosarev I.N., Starikovskii A.Yu. // J. Phys. D: Appl. Phys. 2008. Vol. 41. No. 21. 215207.

3. Kosarev I.N., Aleksandrov N.L., Kindysheva S.V., Starikovskaia S.M., Starikovskii A.Yu. Kinctic mechanism of plasma-assisted ignition of hydrocarbons // J. Phys.D: Appl. Phys. 2008. Vol. 41. No. 3. 032002.

4. Kosarev I.N., Aleksandrov N.L., Kindysheva S.V., Starikovskaia S.M., Starikovskii A.Yu. // Combustion and Flame. 2008. Vol. 154. No. 3. P. 569 586.

5. Kosarev I.N., Aleksandrov N.L., Kindysheva S.V., Starikovskaia S.M., Starikovskii A.Yu. // Combustion and Flame. 2009. Vol. 156. No. 1. P. 221 -233.

6. Александров ПЛ., Киндышева С.В., Кукаев Е.Н., Стариковская С.М., Стариковский АЛО. // Физика плазмы. 2009. Т. 35. Вып. 10. С. 941 956.

7. Стариковская С.М., Александров Н.Л., Косарев И.Н., Киндышева С.В., Стариковский А.Ю. // Химия высоких энергий. 2009. Т. 43. Вып. 3. С. 259 265.

8. Aleksandrov N.L., Kindysheva S.V., Kosarev I.N., Strarikovskaia S.M., Starikovskii A.Yu. // Proceedings of the Combustion Institute. 2009. Vol. 32. P. 205 212.

9. Aleksandrov N.L., Kindysheva S.V., Nudnova M.M., Starikovskiy A.Yu. // J. Phys. D: Appl. Phys. 2010. Vol. 43. No. 25. 255201.

10. Dyatko N.A., Kochetov I.V., Napartovich A.P. // J. Phys. D: Appl. Phys. 1993. Vol. 26. No. 3. P. 418-423.

11. Дятко H.A., Кочетов И.В., Напартович А.П. // Физика плазмы. 2002. Т. 28. Вып. 11. 1046- 1053.

12. Dyatko N.A., Ionikh Y.Z., Kolokolov N.B., Meshchanov A.V., Napartovich A.P. // IEEE Transactions on Plasma Science. 2003. Vol. 31 No. 4. P. 553 563.

13. Sa P.A., Guerra V., Loureiro J., Sadeghi N. // J. Phys. D: Appl. Phys. 2004. Vol. 37. No. 2. P. 221-231.

14. Вихарсв А.Л., Горбачев A.M., Иванов О.И., Исаев В.А., Колганов Н.Г., Колыско A.JL, Офицеров М.М. // ЖТФ. 1996. Т. 66. Вып. 7. С. 56 72.

15. Chang Jen Shin, Hobson R.M., Laframboise J.G. and Organ G.L. // J. Phys. B: Atom. Molec. Phys. 1978. Vol. 11. No. 9. P 1675- 1679.

16. Хаксли JI. и Кромптон Р. Диффузия и дрейф электронов в газах. М: Мир, 1977. 672 с.

17. Бычков В.Л., Елецкий A.B. // Химия плазмы. Москва: Энергоатомиздат. 1985. Вып. 12, С. 119.

18. Смирнов Б.М. Ионы и возбужденные атомы в плазме. М.: Атомиздат, 1974. 456 с.

19. Вихарев А.Л. Иванов О. А., Степанов А.Н. //ЖТФ. 1984. Т. 54. Вып. 8. С. 1617-1619.

20. Александров П.Л., Киндышева C.B., Кирпичников A.A., Косарев И.Н., Стариковская С.М., Стариковский А.Ю. // XXXVIII Международных конференции по физике плазмы и управляемому термоядерному синтезу (Звенигород 2011, 14 — 18 февраля).

21. Aleksandrov N.L., Anokhin Е.М., Kindysheva S.V., Kirpichnikov A.A., Kosarev I.N., Nudnova M.N., Starikovskaia S.M., Starikovskiy A.Yu. // 10th International Workshop on Magneto Plasma Aerodynamics (Moscow, 21-24 March, 2011).

22. Starikovskaia S.M. // J. Phys. D: Appl. Phys. 2006. Vol. 39. No. 16. P. R265 -R299.

23. Попов H.A. // TBT. 2007. T. 45. C. 296 315.

24. Leonov S.B., Yarantsev D.A., Napartovich A.P., Kochetov I.V. // IEEE Transactions on plasma Science. 2006. Vol. 34. No. 6. P. 2514 2525.

25. Кочетов И.В., Напартович А.П., Леонов С.Б. // Химия высоких энергий. 2006. Т. 40. Вып. 2. С. 1 8.

26. Naidis GN.il J. Phys. D: Appl. Phys. 2007. Vol. 40. No. 15. P. 4525-4531.

27. Shibkov V. M., Aleksandrov A.F., Chernikov V.A., Ershov A.P., Shibkova L.V. // Journal of Propulsion and Power. 2009. Vol. 25. No. 1. P. 123 137.

28. Напартович A.П., Кочетов И.В., Леонов С.Б. // Теплофизика Высоких Температур. 2010. Т. 48. Вып. 1.С. 60-66.

29. Adamovich I.V., Choi I., Jiang N., Kim J H., Keshav S., Lempert W.R., Mintusov E., Nishihara M., Samimy M., Uddi M. // Plasma Sources Sei. Technol. 2009. Vol. 18. No. 3. 034018.

30. Starik A.M., Kozlov V.E., Titova N.S. // Combustion and Flame. 2010. Vol. 157. P.313 -327.

31. Williams S., Popovic S., Vuskovic L., Carter C., Jacobson L., Kuo S., Bivolaru D., Corera S., Kahandawala M., Sidhu S. // Proc. 42nd AIAA Aerospace Sciences Meeting and Exhibit. (Reno, NV, 5 8 January 2004) paper AIAA - 2004 - 1012.

32. Bozhenkov S., Starikovskaia S., Starikovskii A.Y. // Combustion and Flame. 2003. Vol. 133. P. 133-146.

33. Starikovskaia S.M., Kukaev E.N., Kuksin A.Yu., Nudnova M.M., Starikovskii A.Yu. // Combust. Flame. 2004. Vol. 139. No. 3. P. 177- 187.

34. Fridman A. Plasma Chemistry. New York: Cambridge University Press, 2008. 978 p.

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

36. Marode Е., Bastien F., Bakker М. // J. Appl. Phys. 1979. Vol. 50. P. 140.

37. Bletzinger P., Ganguly B.N., Van Wie D., Garscadden A. // J. Phys. D: Appl. Phys. 2005. Vol. 38. No. 4. P. R33 R58.

38. Moreau E. // J. Phys. D: Appl. Phys. 2007. Vol. 40. No. 3. P. 605 636.

39. Попов H.A. Исследование механизма быстрого нагрева аота и воздуха в газовых разрядах // Физика плазмы. 2001. Т. 27. Вып. 10. С. 940 949.

40. Boeuf J.P., Kunhardt Е.Е. // J. App. Phys. 1986. Vol. 60. P. 915.

41. Александров А.Ф., Кузовников A.A., Шибков В.М. // Инж. физ. журн. 1992. Т. 62. С. 726.

42. Бердышев A.B., Вихарев А.О., Гитлин М.С., Дерюгин A.A., Иванов O.A., Кочетов И.В., Литвак А.Г., Напартович А.П., Полушкин И.Н., Степанов А.Н., Щербаков А.И. // Теплофизика Высоких Температур. 1988. Т. 26. Вып. 4. С. 661 666.

43. Anokhin Е.М., Starikovskaia S.M., Starikovskii A.Yu. // 42nd AIAA Aerospace Sciences Meeting and Exhibit (Reno, NV, 5-8 January 2004) paper AIAA -2004 674 .

44. Знаменская И.А., Латфулин Д.Ф., Лутский А.Е., Мурсенкова И.В., Сысоев Н.Н. // ЖТФ. 2007. Т. 77. Вып. 5. С. 10.

45. Roupassov D., Nudnova М., Nikipelov A., Starikovskii А. // 46th AIAA Aerospace Sciences Meeting and Exhibit (Reno, NV, 7-10 January 2008) paper AIAA 2008 - 1367.

46. Roupassov D.V., Nikipelov A.A., Nudnova M.M., Starikovskii A.Yu. AIAA J. 2009. Vol. 47. P. 168.

47. Starikovskii A.Yu., Nikipelov A.A., Nudnova M.M., Roupassov D.V. // Plasma Sources Sci. Technol. 2009. Vol. 18. No. 3. 034015.

48. Popov N.A. // J. Phys. D: Appl. Phys. 2011. Vol. 44. 285201.

49. Mintoussov E.I., Pendleton S.J., Gerbault F.G., Popov N.A., Starikovskaia S.M. // J. Phys. D: Appl. Phys. 2011. Vol. 44. 285202.

50. Flitti A., Pancheshnyi S. // Eur. Phys. J. Appl. Phys. 2009. Vol. 45. 21001.

51. Aleksandrov N.L., Kirpichnikov A.A., Kindysheva S.V., Kosarev I.N., Starikovskii A.Yu. // 45th AIAA Aerospace Science Meeting and Exhibit (Reno, NV, 8-11 January 2007) paper AIAA 2007 - 0997.

52. МакИвен M., Филлипс Л. // Химия атмосферы. М.: Мир, 1978. 204 с.

53. Смирнов Б.М. // Комплексные ионы. Москва: Наука, 1983. 149 с.

54. Kossyi I.A., Kostinsky A.Yu., Matveev A.A., Silakov V.P. // Plasma Sourse Sci. Technol. 1992. Vol. 1. No. 3. P. 207 220.

55. Mitchell J.B.A. // Phys. Rep. 1990. Vol. 186. P. 215-248.

56. Александров Н.Л., Кончаков A.M., Шачкин Л.В., Шашков B.M. // Физика плазмы. 1986. Т. 12. С. 1218- 1224.

57. Maier H.N., Fessenden R.W. // J. Chem. Phys. 1975. Vol. 62. No. 12. P. 4790 4794.

58. Warman J.M., Sennhauser E.S., Armstrong D.A. // J. Chem. Phys. 1979. Vol. 70. No. 2. P. 995 999.

59. Senhauser E.S., Armstrong D.A. and Warman J.M. // Radiat. Phys. Chem. 1980. Vol. 15. P. 479-483.

60. Johnsen R. // J. Chem. Phys. 1993. Vol. 98. No. 7. P. 5390 5395.

61. Florescu Mitchell A.I., Mitchell J.B.A. // Phys. Rep. 2006. Vol. 430. P. 277 - 374.

62. Смирнов Б.М. Отрицательные ионы. M.: Атомиздат, 1978. 176 с.

63. Hagelaar G.J.M., Pitchford L.C. // Plasma Sources Sci. Technol.2005. Vol. 14. P. 722 733.

64. Dyatko N.A., Kochetov I.V., Napartovich A.P., Sukharev A.G. EEDF: the software package for calculations of the electron energy distribution function in gas mixtures. http://www.lxcat.laplace.univ tlse.fr/software/EEDF/

65. Биберман JT.M., Воробьев B.C., Якубов И.Т. Кинетика неравновесной плазмы. М.: Наука, 1982. 378 с.

66. Collins С.В. // Phys. Rev. 1965. Vol. 140. P. A1850 A1857.

67. Jungen С., Pratt S.T. // Phys. Rev. Lett. 2009. Vol. 102. 023201.

68. Johnsen R., Guberman S.L. // Adv. At. Mol. Opt. Phys. 2010. Vol. 59. P. 75.

69. Cunningham A.J., Hobson R.M. // J. Phys. B: Atom. Molec. Phys. 1972. Vol. 5. No. 12. P. 2320-2327.

70. Радциг A.A. Смирнов Б.М. Справочник по атомной и молекулярной физике. Москва: Атомиздат, 1980. 240 с.

71. Bates D.R. // Adv. Atom. Mol. Phys. New York: Academic. 1979. No. 15. P. 235 262.

72. Напартович А.П., Наумов В.Г., Шашков B.M. // Физика плазмы. 1975. Т. 1. С. 821 -829.

73. Pack J.L., Phelps A.V. // J. Chem. Phys. 1966. Vol. 45. No. 11. P. 4316 4329.

74. Pack J.L., Phelps A.V. // J. Chem. Phys. 1966. Vol. 44. No. 5. P. 1870 1883.

75. Eliasson В., Kogelschatz U. Basic data for modeling of electric discharge in gases: oxygen, Report KLR-86-11С // Report Brown Boveri Forschungszentrum CH 5405 Baden, 1986.

76. Braginskiy O.V., Vasilieva A.N., Klopovskiy K.S., Kovalcv A.S., Lopaev D.V., Proshina O.V., Rakhimova T.V., Rakhimov A.T. // J. Phys. D Appl. Phys. 2005. Vol. 38. No. 19. P. 36093625.

77. Ionin A.A., Kochetov I.V., Napartovich A.P., Yuryshev N.N. // J. Phys. D Appl. Phys. 2007. Vol.40. No. 2. P. R25-R61.

78. Delcroix J.L., Ferreira C.M., Ricard A. // Principles of Laser Plasmas. New York: Wiley, 1976. P. 176 — 244.

79. Справочник констант элементарных процессов с участием атомов, ионов, электронов, фотонов. / Под ред. А.Г. Жиглинского. Санкт-Петербург: Издательство Санкт-Петербургского университета , 1994. 335 с.

80. Balamuta J., Golde M.F., Moyle A.M. // J. Chem. Phys. 1985. Vol. 82. No. 7. P. 3169 3178.

81. Velazco J. E., Kolts J.H., Setser D. W. // J. Chem. Phys. 1978. Vol. 69. No. 10. P. 4357 -4373.

82. Velazco J. E., Kolts J.H., Setser D. W. // J. Chem. Phys. 1976. Vol. 65. No. 9. 3468 3480.

83. Praxmarer C., Hansel A., Lindinger W., Herman Z. // J. Chem. Phys. 1998. Vol. 109. No. 11. P. 4246-4252.

84. Bates D.R. // Astrophys. J. 1986. Vol. 306. P. L45 L47.

85. Janev R. K., Reiter D. // Phys. Plasmas. 2004. Vol. 11. P. 780 829.

86. Шкаровский И., Джонстон Т., Бачинский М. Кинетика частиц плазмы. М.: Атомиздат, 1969. 396 с.

87. Tachibana К. // Phys. Rev. А. 1989. Vol. 34. P. 1007 1015.

88. Bychkov V.L., Kochetov I.V., Bychkov D.L., Volkov S.A. // IEEE Trans. Plasma Sei. 2009. Vol. 37. 2280 2285.

89. Iannuzzi M.P., Jeffries J.B., Kaufman F. // Chem. Phys. Lett. 1982. Vol. 87. P. 570.

90. DeSousa A.R., Touzeau M., Petitdidier M. // Chem. Phys. Lett. 1985. Vol. 121. P. 423.

91. Piper L.G. // J. Chem. Phys. 1988. Vol. 88. No. 11. P. 6911 6921.

92. Piper L.G.//J. Chem. Phys. 1988. Vol. 88. No. l.P. 231 -239.

93. Словецкий Д.И. Механизмы химических реакций в неравновесной плазме. М.: Наука, 1980. 310 с.

94. Piper L.G., Caledonia G.E., Kennealy J.P. // J. Chem. Phys. 1981. Vol. 75. No. 6. P. 2847 -2852.

95. Piper L.G. // J. Chem. Phys. 1987. Vol. 87. No. 3. Р. 1625 1629.

96. Herron J.T.// J. Phys. Chem. Ref. Data. 1999. Vol. 28. P. 1453.

97. Pintassilgo C.D., Loureiro J., Cernogora G., Touzeau M. // Plasma Sources Sei. Technol. 1999. Vol. 8. P. 463.

98. Piper L.G. // J. Chem. Phys. 1992. Vol. 97. No. 1. P. 270 275.

99. Marinelli W.J., Kessler W.J., Green B.D., Blumberg W.A.M. // J. Chem. Phys. 1989. Vol. 90. P. 2167-2173.

100. Konnov A.A. // Proc. of the Combustion Institute. Pittsburg. 2000. Vol. 28. P. 317.103. www.me.berkeley.edu/gri mech.

101. Magne L., Pasquiers S., Gadonna K., Jeanney Р., Blin Simiand N., Jorand F., Postel C. // J. Phys. D: Appl. Phys. 2009. Vol. 42. No. 16. 165203.

102. Capitelly M., Ferreira C.M., Gordiets B.F., Osipov A.I. Plasma Kinetics in Atmospheric Gases. Berlin: Springer, 2000. 301 p.

103. Cao Y.S., Johnsen R. // J. Chem. Phys. 1991. Vol. 95. No. 10. P. 7356.

104. Bates D.R. //J. Phys. B: Atom. Molec. Phys. 1981. Vol. 14. No. 18. P. 3525 3534.

105. Александров Н.Л.//УФН. 1988. Т. 154. С. 177.

106. Phelps A.V., Pitchford L C. // Phys. Rev. A. 1985. T. 31. C. 2932.

107. Bates D.R. // J. Phys. B: Atom. Molec. Phys. 1980. Vol. 13. No. 2. P. L51 L56.

108. Bohringer II., Durup Ferguson M., Fahey D.W., Fehsenfeid F.C., Ferguson E.E. // J. Chem. Phys. 1983. Vol. 79. No. 9. P. 4201 - 4213.

109. Dobler W„ Ramler H., Villinger H., Howorka F., Lindinger W. // Chem. Phys. Lett. 1983. Vol. 97. P. 553.

110. Durup Ferguson M., Bohringer H., Fahey D.W., Fehsenfeid F.C., Ferguson E.E. // J. Chem. Phys. 1984. Vol. 81. No. 6. P. 2657 - 2666.

111. Federer W., Dobler W., Howorka F., Lindinger W., Durup Ferguson M., Ferguson E.E. // J. Chem. Phys. 1985. Vol. 83. No. 3. P. 1032 - 1038.