Корреляционные методы измерения скорости и температуры высокоскоростного воздушного потока и реагирующих смесей

Фирсов, Александр Александрович

Корреляционные методы измерения скорости и температуры высокоскоростного воздушного потока и реагирующих смесей тема автореферата и диссертации по физике, 01.04.01 ВАК РФ

Фирсов, Александр Александрович АВТОР

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

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

2011 ГОД ЗАЩИТЫ

01.04.01 КОД ВАК РФ

Диссертация по физике на тему «Корреляционные методы измерения скорости и температуры высокоскоростного воздушного потока и реагирующих смесей»

Автореферат диссертации на тему "Корреляционные методы измерения скорости и температуры высокоскоростного воздушного потока и реагирующих смесей"

005002016

Фирсов Александр Александрович

КОРРЕЛЯЦИОННЫЕ МЕТОДЫ ИЗМЕРЕНИЯ СКОРОСТИ И ТЕМПЕРАТУРЫ ВЫСОКОСКОРОСТНОГО ВОЗДУШНОГО ПОТОКА И РЕАГИРУЮЩИХ СМЕСЕЙ

01.04.01 - Приборы и методы экспериментальной физики

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

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

2 4 НОЯ 2011

Москва-2011

005002016

Работа выполнена в Учреждении Российской академии наук Объединенном институте высоких температур РАН (ОИВТ РАН) и на кафедре Физической электроники Физического факультета МГУ имени М.В. Ломоносова.

Научный руководитель: доктор физико-математических наук

Леонов Сергей Борисович

Научный консультант: доктор физико-математических наук, профессор

Александров Андрей Федорович

Официальные оппоненты: доктор физико-математических наук, профессор

Голуб Виктор Владимирович

кандидат технических наук Гольдфельд Марат Абрамович

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

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

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

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

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

Диссертационный совет Д 501.001.66 Ученому секретарю Карташову И.Н.

119991, ГСП-1, Москва, Ленинские горы, МГУ имени М.В. Ломоносова, дом 1, строение 2, Физический факультет.

Автореферат разослан «15» ноября 2011 г.

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

диссертационного совета Д 501.001.66 кандидат физико-математических наук

Карташов И.Н.

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

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

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

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

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

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

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

• Исследование динамики температуры и скорости предварительно нагретого высокоскоростного потока газа.

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

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

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

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

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

• Впервые проведено экспериментальное определение скорости сверхзвукового потока при плазмо-инициированном взаимодействии воздуха с водородом и этиленом в области, частично занятой плазмой разряда постоянного тока. Показано, что при скорости набегающего сверхзвукового потока 495 м/с скорость газа в области разряда составляет 180-230 м/с, а скорость газа в области химических реакций находится в интервале 100- 150 м/с.

Личный вклад автора

Все представленные в работе экспериментальные результаты получены автором самостоятельно и в составе научных коллективов Физического факультета МГУ и ОИВТ РАН. Разработаны уникальные методики, алгоритмы и программы обработки экспериментальных данных. Обработка экспериментальных данных и расчеты выполнены автором самостоятельно. Выполнено численное моделирование трехмерного течения в канале с внезапным расширением и в случае плазмоинициированного горения водорода на плоской стенке в сверхзвуковом потоке воздуха. Проведено сравнение экспериментальных данных, полученных различными методами. На основании полученных результатов автором сформулированы и обоснованы выводы диссертационной работы.

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

Полученные экспериментальные результаты являются важными для развития плазменной аэродинамики, газодинамики и физики плазмы. Результаты данной диссертационной работы могут быть использованы при исследовании различных плазменных образований в потоке газа, в экспериментах по плазмо-инициированному горению топливных смесей, а также в областях практической аэродинамики. Представленные результаты также могут быть применены для оптимизации лабораторных установок, предназначенных для исследования предварительно подогретых потоков газа, в том числе при наличии плазмы в потоке. Результаты работы могут быть использованы в ЦАГИ, ЦИАМ, Институт механики МГУ, ИТПМ СО РАН и ряде других организаций.

Достоверность полученных результатов

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

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

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

«Ломоносов», МГУ, Москва, 2007, 2008; Школа-семинар по магнито-плазменной аэродинамике, ОИВТ РАН, Москва, 2008; AIAA Plasmadynamics and Lasers Conference, San Antonio, Texas, USA, 2009; Международное Совещание по Магнитоплазменной Аэродинамике, ОИВТ РАН, Москва, 2009, 2010, 2011; AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Orlando, Florida, USA, 2010, 2011; IEEE International Conference on Plasma Science, Norfolk, Virginia, USA, 2010; Школа молодых ученых «Актуальные проблемы физики», ФИАН, Звенигород-Москва, 2010; AIAA International Space Planes and Hypersonic Systems and Technologies Conference, San Francisco, California, USA, 2011; Инженерные системы, РУДН, Москва, 2011;

Публикации

По материалам диссертации опубликовано 4 статьи в реферируемых журналах, входящих в перечень ВАК, и 16 статей в сборниках тезисов и трудов конференций. Основные результаты, представленные в диссертации, отражены в публикациях автора, список которых приведен в конце автореферата.

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

Диссертационная работа состоит из введения, четырех глав, заключения и списка литературы. Объем диссертации составляет 105 страниц, включая 73 рисунка и библиографию из 75 наименований.

СОДЕРЖАНИЕ ДИССЕРТАЦИИ

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

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

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

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

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

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

Рис. 1. Схема импульсной аэродинамической трубы (МГУ)

1 - компрессор; 2 - баллоны; 3 - электромагнитный клапан;

4 - термостатированный нагреватель; 5 - термопара; 6 - сверхзвуковое сопло;

7 - система подачи топлива; 8 - электроды; 9 - датчики давления;

10 - формирующая секция; 11 - тестовые секции № 1 -№4 (слева направо);

12 - барокамера; 13 - вакуумный насос

В разделе 2.1 приведено описание экспериментальной установки, а также используемого оборудования и методов диагностики. Исследования возможности применения импульсного поперечного разряда для определения скорости потока проводились на экспериментальном стенде кафедры Физической электроники МГУ имени М.В. Ломоносова. Аэродинамический

канал оснащен соплом, рассчитанным на число Маха М = 2, его выходной диаметр равен 25 мм. За секцией сопла размещена формирующая секция, представляющая собой канал с внутренним диаметром 25 мм и длиной 250 мм. За формирующей секцией следуют 4 тестовых секции прямоугольного сечения 25x40 мм и длиной 200 мм каждая. Аэродинамическая труба оснащена 16-ти канальным аналого-цифровым преобразователем, подключенным к компьютеру и предназначенным для сбора данных с термопары и датчиков статического и полного давления. Схема экспериментальной установки представлена на Рис. 1.

В разделе 2.2 представлено обоснование метода импульсного разряда и описана методика измерения скорости. В экспериментах использован разряд малой длительности импульса Д1 = 20 - 100 мкс при частоте повторения 3050 Гц, что позволяет утверждать, что в среднем мощность разряда мала по сравнению с мощностью набегающего потока. Электроды располагались так, что пробой происходил поперек потока в плоскости, равноудаленной от боковых (больших) стенок канала. Скорость потока определяется как отношение продольной протяженности разряда в потоке на интегральной фотографии (т.е. расстояния, на которое сносится поперечный потоку фрагмент разрядного канала) к длительности импульса разряда. За один пуск воздуха по каналу производится серия фотографий, по каждой из которых производится определение протяженности разряда и скорости потока, т.е. в результате получается зависимость скорости от времени для каждого пуска. Затем вычисляются средние скорость потока и протяженность разряда за пуск.

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

В разделе 2.4 представлено сравнение результатов, полученных методом импульсного разряда, с результатами численного моделирования в программном комплексе РЬ^Л^вюп и газодинамическими измерениями. Приведено описание

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

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

600

"о"

500

£ 400

§ 300

200

ОД 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1

X, (м)

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

В разделе 2.6 представлены результаты определения скорости предварительно подогретого потока газа. Для контроля измерений использовалась термопара, расположенная на стенке дозвуковой части канала. За характерное время пуска (1-2 с) термопара не успевала выходить на насыщение, что потребовало применения компьютерной обработки кривой прогрева. Однако полученные с помощью аппроксимации значения температуры оказались ниже ожидаемых теоретических, рассчитанных через скорость потока. Полученный результат наглядно демонстрирует сложности определения 1 температуры потока в установках кратковременного действия. В разделе показано, что метод импульсного разряда применим в условиях подогретого потока воздуха и позволяет оценить полную температуру газа. На представленном на Рис. 3 графике приведены значения скорости, определенные экспериментально методом импульсного разряда, и соответствующие им значения температуры: измеренные с помощью термопары и теоретические, т.е. вычисленные через скорость потока.

Скорость (по разряду) в 1й секции, (м/с) Рис. 3. Результаты измерений в потоке предварительно нагретого воздуха

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

В разделе 3.1 приведено описание экспериментальной установки и используемого оборудования. Исследования возможности применения лазерной системы для определения скорости потока проводились на экспериментальной

установке ИАДТ-50Н лаборатории Плазменной аэродинамики Объединенного института высоких температур РАН. Аэродинамический канал оснащен соплом, рассчитанным на число Маха М = 2, его выходное сечение составляет 60x72 мм и эквивалентно сечению тестовой секции. Конструкция тестовой секции позволяет измерять статическое давление на стенке канала, устанавливать различные вставки с электродами и трубками подачи топлива для проведения экспериментов по воспламенению топлива в сверхзвуковом потоке.

✓ —

Рис. 4. Схема импульсной аэродинамической трубы ИАДТ-50Н (ОИВТ РАН)

1 - цистерна высокого давления, оснащенная компрессором; 2 - электромагнитный клапан;

3 - форкамера; 4 - сопло; 5 - рабочая секция; 6 - электроды; 7 - трубки подачи топлива;

8 - баллон с топливом для эксперимента; 9 - барокамера, оснащенная вакуумным насосом;

10 - подогреватель воздуха; 11 - баллон с топливом для подогревателя;

12 - лазерная рефракционная система; 13 - осциллограф

В разделе 3.2 представлено описание взаимно-корреляционного метода определения скорости сверхзвукового потока по сигналам рефракционных датчиков. Метод основан на определении времени Аи которое требуется возмущениям плотности в потоке для преодоления известного расстояния Ь между двумя лазерными лучами. После прохождения сверхзвукового потока лучи попадают на фотодиоды, которые детектируют изменение интенсивности излучения, связанное с отклонением луча на градиентах плотности в потоке. Время А1 равно интервалу сдвига локального максимума взаимной корреляционной функции (ВКФ) для пар фрагментов и и2 регистрируемых рефракционной лазерной системой сигналов: В12(т) = и1(£)и2(р + т) йь. Фрагмент ВКФ показан на Рис. 5.

~ 0,5

0,0

-0,5

-0,4 -0,2 0,0 0,2 0,4

Сдвиг T (MC)

Рис. 5. Взаимная корреляционная функция (фрагмент)

В разделе анализируется влияние на результаты измерений интервала поиска локального максимума ВКФ и интервала разбиения исходного сигнала на фрагменты U, и U2, по которым рассчитывается взаимная корреляционная функция. Интервал поиска локального максимума задается через ожидаемые минимальную и максимальную скорости. Показано, что сужение интервала ожидаемых скоростей к реальному диапазону приводит к уменьшению количества ошибок автоматической обработки результатов эксперимента (т.е. значений, более чем на 25% отличающихся от средней величины). Увеличение 1 длительности интервала разбиения с 0,2 мс до I мс влияет аналогичным образом. Необходимо отметить, что дальнейшее увеличение интервала разбиения не уменьшает количество ошибок обработки, однако при этом , ухудшается временное разрешение корреляционного метода.

В разделе 3.3 представлены результаты измерений скорости сверхзвукового потока. Для потока без предварительного нагрева проведено сравнение результатов: для канала с М = 2 скорость потока, определенная предложенным методом, составила V = 495130 м/с, что хорошо согласуется I с результатом V = 487±7 м/с, полученным через газодинамические функции по ' измеренному статическому давлению.

При использовании в ИАДТ-50Н подогревателя воздуха измерена скорость потока предварительно нагретого газа. В нормальном режиме работы подогревателя нагрев газа приводил к увеличению скорости до 550 - 750 м/с. В данных опытах температура предварительно нагретого газа находилась

в интервале 360 - 670 К. В отдельных случаях использования подогревателя I нагрев приводил к ускорению газа до скорости 1200 м/с (нестационарный режим работы). Зависимость скорости от времени при работе подогревателя

Время, (с)

Рис. 6. Скорость потока при работе подогревателя

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

разряда. Показано, что электрический пробой приводит к повышению амплитуды спектра в широком диапазоне частот от 100 до 400 кГц. В случае

последовательности из трех электрических разрядов (с интервалами между импульсами 110-140мкс) повышение амплитуды на указанных частотах -происходит в более широком временном интервале и возмущения практически сливаются в одно. Графики зависимости среднего значения десятичного логарифма амплитуды в заданной полосе частот 2АГ от времени показаны

на Рис. 7.

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

Рис. 7. Повышение амплитуды колебаний при прохождении возмущений

а - Один разрядный импульс (2ДГ = 40 кГц); б - три импульса (2АГ = 25 кГц)

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

В разделе 4.1 приведено описание дополнительных элементов экспериментальной установки и используемого в данной главе оборудования. Получение изображений для дальнейшего корреляционного анализа I осуществлялось с помощью оптической шлирен-системы, регистрация | проводилась с помощью камеры линейного сканирования (КЛС), т.е. за один кадр получалась одна строчка пикселей. Под изображением с КЛС

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

В разделе подробно описана схема расположения трубок подачи топлива ' и электродов. В нижней стенке рабочей секции была установлена керамическая I вставка с пятью отверстиями 0 4 мм для подачи топлива - водорода или | этилена. Вставка также содержала 7 медных электродов 0 2 мм, торцы которых

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

В разделе 4.2 представлено описание корреляционного метода обработки изображений. Метод основан на определении смещения ёх неоднородности плотности в потоке за известное время Л, т.к. известна частота работы камеры линейного сканирования (50 или 100 кГц). Для определения смещения ёх была использована двумерная взаимная корреляционная функция В12(а,т) = М^™ и1(х,1)и2(х + а^ + т^йхМ, где и, и и2 - фрагменты исходного изображения, смещенные друг относительно друга по времени на одну строку. Типичный вид исходного изображения КЛС и двумерная взаимная корреляционная функция представлены на Рис. 8 и Рис. 9 соответственно.

X (вдоль потока)

Рис. 8. Типичный вид изображения КЛС

(по вертикали последовательные моменты времени)

Сдвиг между двумя соседними максимумами по а и есть искомое dx. Скорость определялась с помощью программы, реализованной в программном комплексе Lab VIEW. В общем случае результат представлял собой изменение профиля скорости во времени. Приведены результаты измерений скорости над i профилем крыла в дозвуковом потоке, а также в зоне плазмо-инициированного горения топлива (этилена или водорода) в сверхзвуковом потоке воздуха.

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

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

Координата со значением 0 мм на Рис. 10 соответствует центру трубки подачи I топлива. Полученные профили скорости качественно совпадают, однако количественные значения незначительно расходятся.

350

область присутствия топлива

50 -

-•— Численное моделирование ♦ Корреляционная обработка изображения

-20

-10

X, (мм)

Рис. 10. Продольный профиль скорости в области плазмо-инициированного горения водорода

В заключении сформулированы основные выводы диссертационной работы (положения, выносимые на защиту):

• Проведено исследование зависимости продольной длины импульсного электрического разряда в высокоскоростном воздушном потоке от параметров разряда. Показано, что при длительности импульса 20 - 30 мкс разряд может быть применен для измерения скорости потока. Разработан и успешно применен метод определения скорости по протяженности импульсного разряда для холодного и предварительно нагретого (до -950 К) воздушных потоков, а также в потоке реагирующей смеси газов. Значения скоростей, измеренных предложенным методом, находились в диапазоне от 250 до 850 м/с.

• Разработана новая модификация корреляционного метода определения скорости сверхзвукового потока, позволяющая проводить измерения по естественным маркерам (без внесения искусственных возмущений в поток). Метод успешно применен для холодного и предварительно нагретого воздушных потоков в диапазоне скоростей от 490 до 1200 м/с.

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

• Разработан новый метод определения скорости потока в присутствии плазмы электрического разряда и в области химических реакций в потоке посредством корреляционного анализа шлирен-изображений с камеры линейного сканирования. Впервые проведены прямые измерения скорости реагирующего газа при плазмо-инициированном горении газообразных топлив в сверхзвуковом потоке воздуха. Показано, что при скорости сверхзвукового потока 495 м/с скорость газа в области разряда составляет 180-230 м/с, а в области горения - 100-150 м/с.

Лицензия на программный комплекс вычислительной аэро-гидродинамики FlowVision предоставлена компанией ТЕСИС в рамках программы поддержки учащихся высших учебных заведений «Учись, студент!».

ОСНОВНЫЕ РАБОТЫ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ

Статьи в журналах, входящих в перечень ВАК:

1) Ершов А.П., Каменщиков С.А., Колесников Е.Б., Логунов А.А., Фирсов А.А., Черников В.А. — О возможности измерения скорости потока с помощью маломощного импульсно-периодического разряда // Изв. РАН. МЖГ. 2008. №4. С. 119-127.

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

3) Леонов С.Б., Фирсов А.А., Яранцев Д.А. - Определение скорости сверхзвукового потока по функции корреляции сигналов рефракционных датчиков // Письма в ЖТФ. 2009. том 35. В. 12. С. 25-31.

4) Леонов С.Б., Савелкин КВ., Фирсов А.А., Яранцев Д.А. - Зажигание топлива и стабилизация фронта пламени в сверхзвуковом потоке при помощи электрического разряда // ТВТ. 2010. Т.48. №6. С. 941-947.

Другие публикации:

1) Александров А.Ф., Леонов С.Б., Фирсов А.А. - Численное моделирование горения водорода в сверхзвуковом воздушном потоке // Инженерные системы -2011, РУДН, Москва, 2011.

2) Leonov S.B., Firsov А.А., Shurupov М.А. - Flow structure analysis by cross-correlation image processing // The 10th International Workshop on MagnetoPlasma Aerodynamics, JIHT RAS, Moscow, 2011.

3) Leonov S.B., Isaenkov Yu.I., Firsov A.A., Yarantsev D.A., Shneider M.N. - HighPower Filamentary Pulse Discharge in Supersonic Flow // 48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Orlando, Florida, AIAA 2010-259.

4) Leonov S.B., Isaenkov Yu.I., Firsov A.A., Yarantsev D.A. - Mixing acceleration in high-speed flow by long-spark electrical discharge // The 8th International Workshop on Magneto-Plasma Aerodynamics, Moscow, 2009, proceedings, p.61-67.

5) Leonov S.B., Isaenkov Yu.I., Firsov A.A.— Mixing Intensification in High-Speed Flow by Unstable Pulse Discharge // 40th AIAA Plasmadynamics and Lasers Conference, San Antonio, Texas, AIAA 2009-4074.

6) Яранцев Д.А., Леонов С.Б., Фирсов А.А. - Особенности применения диагностики в экспериментах по взаимодействию электрического разряда со сверхзвуковым потоком // Третья школа-семинар по магнито-плазменной аэродинамике, ОИВТ РАН, Москва, 2008 г., с. 171.

7) Фирсов А.А. - Применение взаимной корреляционной функции для определения скорости сверхзвукового потока // Материалы докладов XV Международной конференции студентов, аспирантов и молодых ученых «Ломоносов», МГУ имени М. В. Ломоносова, Москва, 2008 г., ч. 21-8, с. 19.

Подписано в печать 10.11.2011 г. Печать лазерная цифровая Тираж 100 экз.

Типография Aegis-Print 115230, Москва, Варшавское шоссе, д. 42 Тел.: 8 (495) 785-00-38, 8 (926) 850-53-16 www.autoref.ae-print.ru

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

Оглавление.

Введение.

Глава 1. Обзор литературы.

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

1.2. Методы измерения скорости.

Глава 2. Определение скорости потока с помощью импульсного разряда.

2.1. Описание экспериментальной установки.

2.2. Обоснование метода и методика измерения скорости.

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

2.4. Сравнение результатов метода импульсного разряда с результатами других методов.

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

2.6. Определение скорости потока в случае предварительного нагрева воздуха, поступающего в канал.

Глава 3. Рефракционный лазерный метод.

3.1. Описание экспериментальной установки.

3.2. Описание взаимно-корреляционного метода и влияние ряда параметров на результаты.

3.3. Применение корреляционного метода для случая холодного и предварительно нагретого потока.

3.4. Лазерный метод определения спектра возмущений.

Глава 4. Корреляционный метод обработки изображения.

4.1. Описание экспериментальной установки.

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

4.3. Сравнение профилей с результатами численного моделирования.

Введение диссертация по физике, на тему "Корреляционные методы измерения скорости и температуры высокоскоростного воздушного потока и реагирующих смесей"

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

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

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

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

• Исследование динамики температуры и скорости предварительно нагретого высокоскоростного потока газа.

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

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

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

Личный вклад автора

Все представленные в работе экспериментальные результаты получены автором самостоятельно и в составе научных коллективов Физического факультета МГУ и

ОИВТ РАН. Обработка экспериментальных данных и расчеты выполнены автором самостоятельно.

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

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

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

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

• Экспериментально определена скорость холодного и предварительно нагретого сверхзвуковых воздушных потоков с помощью созданного комплекса.

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

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

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

Полученные экспериментальные результаты являются важными для развития плазменной аэродинамики, газодинамики и физики плазмы. Результаты данной диссертационной работы могут быть использованы при исследовании различных плазменных образований в потоке газа, в экспериментах по плазмоинициированному горению топливных смесей, а также в областях практической аэродинамики. Представленные результаты таюке могут быть применены для оптимизации лабораторных установок, предназначенных для исследования предварительно подогретых потоков газа, в том числе при наличии плазмы в потоке. Результаты работы могут быть использованы в ЦАГИ, ЦИАМ, Институт’ механики МГУ, ИТПМ СО РАН и ряде других организаций.

Достоверность полученных результатов

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

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

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

Ломоносов», МГУ, Москва, 2007, 2008; Школа-семинар по магнито-плазменной аэродинамике, ОИВТ РАН, Москва, 2008; AIAA Plasmadynamics and Lasers Conference, San Antonio, Texas, USA, 2009; Международное Совещание по Магнитоплазменной Аэродинамике, ОИВТ РАН, Москва, 2009, 2010, 2011; А1АА Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Orlando, Florida, USA, 2010, 2011; IEEE International Conference on

Plasma Science, Norfolk, Virginia, USA, 2010; Школа молодых ученых «Актуальные проблемы физики», ФИАН, Звенигород-Москва, 2010; AIAA International Space Planes and Hypersonic Systems and Technologies Conference, San Francisco, California, USA, 2011; Инженерные системы, РУДН, Москва, 2011;

Публикации

Структура диссертации

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

Заключение диссертации по теме "Приборы и методы экспериментальной физики"

Заключение

• Проведено исследование зависимости продольной длины импульсного электрического разряда в высокоскоростном воздушном потоке от параметров разряда. Показано, что при длительности импульса 20-30 мкс разряд может быть применен для измерения скорости потока. Разработан и успешно применен метод определения: скорости по протяженности импульсного разряда для холодного и предварительно нагретого (до ~950 К) воздушных потоков, а также в потоке реагирующей смеси газов. Диапазон скоростей, измеренных предложенным методом, составил от 250 до 850 м/с.

• Разработана новая модификация корреляционного метода определения скорости сверхзвукового потока, позволяющая^ проводить, измерения по естественным маркерам (без внесения искусственных возмущений в поток)' Метод успешно применен для холодного и предварительно нагретого воздушных потоков в диапазоне скоростей от 490 до 1200 м/с. ;

• С , целыо количественного; описания* смешивания впервые проведено экспериментальное изучение влияния: искрового разряда на спектр возмущений плотности сверхзвукового потока. Показано, что разряд создает в потоке неоднородности с характерными размерами 1 - 5 мм.

• Разработан новый метод определения скорости потока в присутствии плазмы электрического разряда и t в области химических реакций в потоке посредством корреляционного анализа шлирен-изображений с камеры линейного сканирования. Впервые проведены прямые измерения скорости реагирующего газа при плазмо-инициированном горении газообразных топлив в сверхзвуковом потоке воздуха: Показано, что при скорости сверхзвукового потока 495 м/с скорость газа в области разряда составляет 180 - 230 м/с, скорость газа в области горения — 100 - 150 м/с.

Автор выражает благодарность научному руководителю д.ф.-м.н. Леонову С.Б. и научному консультанту профессору

Александрову А.Ф., коллегам с кафедры Физической электроники Физического факультета МГУ, в частности профессору Ершову А.П. и доценту Черникову В.А., а также коллегам по лаборатории Плазменной аэродинамики Объединенного института высоких температур РАН, в частности Савелкину К.В. и Яранцеву Д.А., за помощь в проведении экспериментов и обсуждение результатов.

Автор благодарит инжиниринговую компанию ТЕСИС за предоставленную в рамках программы поддержки учащихся высших учебных заведений «Учись, студент!» лицензию на программный комплекс вычислительной аэро-гидродинамики Пом> У1&чоп.

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

1. Финкелъбург В., Меккер Г. Электрические дуги и термическая плазма. / Под ред. Фабриканта В.А. М. : Изд-во иностр. лит. 1961. С. 370.

2. Райзер Ю.П. Лазерная искра и распространение разрядов. М. : Наука. 1974. С. 307.

3. Райзер Ю.П. -Основы современной физики газоразрядных процессов. М.: Наука. 1980. С. 416.

4. А. Райзер Ю.П. Физика газового разряда. М. : Наука. 1992 (le издание 1987 г.). С. 536.

5. Пащенко Н.Т., Райзер Ю.П. -Тлеющий разряд в продольном потоке газа // Физика плазмы. 1982. Т. 8. С. 1086-1092.

6. Райзер Ю.П., Шнейдер М.Н., Яценко H.A. Высокочастотный емкостный разряд: Физика. Техники эксперимента. Приложения. М. : Изд. МФТИ; Наука. 1995. С. 320.

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

8. Велихов Е.П., Голубев B.C., Пашкин С.В. -Тлеющий разряд в потоке газа // Успехи физических наук. 1982. Т. 137. № 1. С. 117-150.

9. Арделян Н.В., Космачевский КВ., Чувашев С.Н. Математическое моделирование энергетически эффективных режимов сверхзвукового движения летательного аппарата с вдувом плазмы (газа) из иглы на носу // Прикладная физика. 1999. №6.

10. Fomin V.M., Tretyakov P.К., Taran J.-P. -Flow control using various plasma and aerodynamic approaches (short review) I I Aerospace Science and Technology. 2004. V. 8. №5. P. 411-421.

11. Bletzinger P., Ganguly B.N., Van Wie D., Garscadden A. Plasmas in high speed aerodynamics // J. Phys. D: Appl. Phys. 2005. V. 38. № 4. P. R33-R57.

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

13. Chernikov VA., Dvinin S.A., Ershov A.P., Timofeev I.B., Shibkov V.M. The Theory of Direct Current Discharge in Transversal Gas Flow // 41st Aerospace Sciences Meeting and Exhibit. 2003. 6-9 Jan. Reno. Nevada. USA. AIAA-2003-1193.

14. Герман В.О., Козлов П.В., Любимов Г.А. Экспериментальное взаимодействие электрической дуги с потоком газа // Изв.РАН. МЖГ. 2005. № 5. С. 169-179.

15. Алферов В.И. Исследование структуры электрического разряда большой мощности в высокоскоростном потоке воздуха // Изв.РАН. МЖГ. 2004. № 6. С. 136-175.

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

17. Абрамович Г.Н. Прикладная газовая динамика. М. : Наука. 1991. С. 600.

18. George A., Wilson W W, Cook R.L. Laser velocimetry measurements in an MHD aerodynamic duct //Plasmadynamics and Lasers Conference. 1992. 6-8 July. Nashville, m AIAA-1992-2986.

19. Beck M.S. Correlation in instruments: cross correlation flowmeters // J. Phys. E: Sci. Instrum. 1981. T. 14. C. 7-19.

20. Kumar R.A., Cook R.L., Norton O.P. -An Extinction Based cross-Correlation Technique for Determination of Flow Velocities // 25th AIAA Plasmadynamics and Lasers Conference. 1994. 20-23 June. Colorado Springs. CO. USA. AIAA 94-2444.

21. Smrekar J., Oman J., SirokB., Mori M. and Hocevar M. A cross-correlation optical velocimeter//Meas. Sci. Technol. 2007. V. 18. № 3. P. 555-560.

22. Norton O., Kumar R. Cross Correlation Velocity Measurements in a Simulated MHD Environment //28th Aerospace Siences Meeting . 1990. 8-11 January. Reno. Nevada. USA. AIAA 90-0350.

23. Norton O.P., Kumar R.A., Srikantaiah D.V. A Direct Comparison of CrossCorrelation Velocity Measurement with Laser Doppler Measurements // 29th Aerospace Sciences Meeting. 1991. 7-10 January. Reno. Nevada. USA. AIAA 910213.

24. Adrian R.J. Twenty years of particle image velocimetry // Experiments in Fluids. 2005. №39. C. 159-169.

25. Prasad, Ajay K. Particle image velocimetry // Current Science. 2000. T. 79. № 1. C. 51-60.

26. Raffel М., Willert C., Wereley S., Kompenhans J. Particle image velocimetry. Berlin : Springer. 2007. C. 448.

27. Петров А.П., Наливайченко Д.Г. Метод искрового трассирования для измерения скорости сверхзвукового потока // X International Conference on the Methods of Aerophisical research. 2000. 9-16 July. Novosibirsk. 158-161 (Proceedings).

28. Запрягаев В.И., Петров А.П., Солотчин A.B. Исследование неравномерности распределения скорости в слое смешения недорасширенной струи методом электроразрядного трассирования потока // Прикладная механика и техническая физика. 2004. Т. 45. № 6. С. 58-64.

29. Каменщиков С.А. Взаимодействие плазмы продольно-поперечного и плазмодинамического разрядов со сверхзвуковым воздушно-пропановым потоком. М.: МГУ. Физический факультет. 2011. С. 125. диссертация на соискание ученой степени к.ф.-м.н.

30. Горлин С.М., Слезингер И.И. Аэромеханические измерения. М. : Наука. 1964. С. 370.

31. Громов В.Г., Ершов А.П., Левин В.А., Шибкое В.М. Поперечные разряды в сверхзвуковых потоках воздуха. Моделирование эффектов, влияющих на нагрев газа в разрядном канале // ТВТ. 2006. Т. 44. № 2. С. 185-194.

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

33. Митин И.В., Русаков B.C. Анализ и обработка экспериментальных данных. М. : Физический факультет МГУ. 2002. С. 44.

34. FlowVision НРС. Руководство пользователя. М.: ТЕСИС. 2011. С. 314. www.flowvision.ru.

35. Колесников Е.Б. Взаимодействие плазмы импульсных разрядов со сверхзвуковыми потоками воздуха. М. : МГУ. Физический факультет. 2009. С. 156. диссертация на соискание ученой степени к.ф.-м.н.

36. Ершов А.П., Каменщиков С.А., Колесников Е.Б., Логунов A.A., Фирсов A.A., Черников В.А. О возможности измерения скорости потока с помощью маломощного импульсно-периодического разряда // Изв.РАН. МЖГ. 2008. № 4. С. 119-127.

37. Ершов А.П., Каменщиков С.А., Колесников Е.Б., Логнов A.A., Фирсов A.A., Черников В.А. Измерение скорости потока с помощью поперечного разряда // Вестник МГУ. серия 3. 2008. № 3. С. 70-72.

38. Leonov S.B., Savelkin K.V., Firsov A.A. -Contactless determination of surfacc temperature within the supersonic channel // The 8th International Workshop on Magneto-Plasma Aerodynamics. 2009. 31 March 2 April. Moscow. P. 127-129.

39. Васильев JI.A. Теневые методы. М. : Наука. 1968. С. 400.

40. Черников В.А. Диагностика импульсной плазмы. М. : Физический факультет МГУ. 2004. С. 180.

41. Зайдель А.Н., Островская Г.В. Лазерные методы исследования плазмы. Л : Наука. 1977. С. 220.

42. Алексеенко М.Д., Бараночников М.Л. Приемники * оптического излучения. Справочник. М.: Радио и связь. 1987. С. 296.

43. Леонов С.Б., Фирсов А.А., Яранцев Д.А. Определение скорости сверхзвукового потока по функции корреляции сигналов рефракционных датчиков // Письма в ЖТФ. 2009. Т. 35. № 12. С. 25-31.

44. Leonov S.B., Isaenkov Y.L., Firsov A.A.,Nothnagel S.L., Gimelshein S.F., Shneider M.N. Jet regime of the afterspark channel decay //Phys. Plasmas. 2010. T. 17. №053505. C. 1-8.

45. Leonov S.B., Isaenkov Y.I., Firsov A.A. Mixing intensification in High-Speed Flow by Unstable Pulse Discharge // 40thiAIAA Plasmadynamics and Lasers Conference. 2009. 22-25 June. San Antonio. Texas. USA. AIAA 2009-4074.

46. Leonov S.B., Firsov A.A., Isaenkov Y.I., Shurupov M.A., Yarantsev D.A., Shneider

47. M.N. Long-spark discharge in mixing layer of two molecular co-flown gases // The 37th IEEE International Conference on Plasma Science. 2010. 20-24 June. Norfolk. Virginia. USA. P. 437. ’

48. Leonov S.B., Isaenkov Y.I., Firsov A.A., Yarantsev D.A. Mixing acceleration in highspeed flow by long-spark electrical discharge // The 8th International Workshop on Magneto-Plasma Aerodynamics . 2009. 31 March - 2 April. JIHT RAS. Moscow. P. 61-67.

49. Кулманаков С.П., Яковлев С.В. Моделирование рабочего процесса сгорания в цилиндре дизеля // Вестник алтайской науки. 2010. № 2 (9). С. 35-38.

50. Firsov A.A , Leonov S.B., Miller A.B. Numerical Simulations of Sypersonic Flow in Covergent-DivergentDuct with Discharge Supply //The 9th International Workshop on Magneto-Plasma Aerodynamics. 2010. 13 - 15 April. Moscow.

51. Фирсов А.А., Леонов С.Б. Управление положением ударной волны с помощью электрического разряда // XIII Школа молодых ученых «Актуальные проблемы физики». ФИАН. 2010. 14-19 ноября. Звенигород-Москва. С.225-227.

52. Leonov S.B., Firsov А.А., Shurupov М.А. -Flow structure analysis by crosscorrelation image processing // The 10th international Workshop on Magneto-Plasma Aerodynamics. 2011. 22-24 March. JIHT RAS. Moscow.

53. Александров А.Ф., Леонов С.Б., Фирсов А.А. -Численное моделирование горения водорода в сверхзвуковом воздушном потоке // Инженерные системы -2011. 5-8 апреля. РУДН. Москва.

54. Леонов С.Б., Савелкин К.В., Фирсов А.А., Яранцев Д.А. Зажигание топлива и стабилизация фронта пламени в сверхзвуковом потоке при помощи электрического разряда // ТВТ. 2010. Т. 48. № 6. С. 941-947.