Исследование воспламенения и горения водорода и метана в газовых потоках при возбуждении электронных степеней свободы молекул кислорода тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

Титова, Наталия Сергеевна АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Москва МЕСТО ЗАЩИТЫ
2012 ГОД ЗАЩИТЫ
   
01.02.05 КОД ВАК РФ
Диссертация по механике на тему «Исследование воспламенения и горения водорода и метана в газовых потоках при возбуждении электронных степеней свободы молекул кислорода»
 
Автореферат диссертации на тему "Исследование воспламенения и горения водорода и метана в газовых потоках при возбуждении электронных степеней свободы молекул кислорода"

МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

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

Специальность 01.02.05 — Механика жидкости, газа и плазмы

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

им. М.В. Ломоносова

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

005009232

Титова Наталия Сергеевна

2 ОЕЗ Ш

Москва-2012

005009232

Работа выполнена в ФГУП "Центральный институт авиационного моторостроения им. П.И. Баранова"

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

A.M. Старик

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

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

Ведущая организация: Учреждение РАН Институт общей физики

им. A.M. Прохорова РАН

Защита состоится 2 марта 2012 г. в 15 часов на заседании диссертационного совета Д.501.001.89 при Московском государственном университете им. М.В. Ломоносова по адресу: 119991, г. Москва, Ленинские горы, д. 1, Главное здание МГУ, аудитория 16-10.

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

Автореферат разослан " ¿Г- января 2012 г. Ученый секретарь

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

А.Н. Осипцов

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

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

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

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

В данной работе рассматривается метод интенсификации горения, основанный на возбуждении метастабильных электронных состояний а!Д„ и 1>'Е„+ молекул Ог, которые с высокой эффективностью можно получить в специально организованном электрическом разряде [6, 7] или при воздействии резонансного лазерного излучения [8]. Поскольку молекулярный кислород (или воздух) является окислителем большинства используемых в промышленности топлив, то данный метод может оказаться эффективным во многих практических приложениях. Всё это определяет актуальность исследований по кинетике процессов в реагирующих средах при наличии в них возбужденных молекул 02(а'л„) и 02(Ь1Т.„+), созданию новых детальных реакционных механизмов в случае возбуждения электронных состояний молекул 02 либо элек-

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

Цель и задачи исследования

Целью работы является численное исследование возможности ускорения воспламенения (в том числе в сверхзвуковом потоке) и увеличения скорости ламинарного пламени в смесях Н2-воздух и СН4-воздух при возбуждении молекул 02 в метастабилыше электронные состояния ¡з'А8 и ¿>'28+, а также определение энергетической эффективности этого метода инициирования горения. Для достижения поставленной цели в работе ставились следующие задачи:

1. разработать детальные кинетические модели окисления водорода и метана в воздухе при наличии в рассматриваемых смесях электронно-возбужденных молекул 02(а'Д8) и 02(б'18+);

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

3. исследовать влияние возбуждения молекул 02 в состояния а'Д8 и

при импульсном воздействии лазерного излучения с длиной волны Л/=1.268 мкм и 762 нм, соответственно, на время задержки воспламенения, а также на температуру воспламенения в смесях Н2-воздух и СН4-воздух;

4. исследовать неравновесные процессы за фронтом наклонной ударной волны в смесях Н2-воздух и СН4-воздух при возбуждении электронных состояний молекул 02 резонансным лазерным излучением и провести анализ возможности существенного сокращения длины зоны индукции в сверхзвуковом потоке за фронтом ударной волны;

5. проанализировать возможность использования метода возбуждения молекул 02 в состояния а'Дц и для ускорения конверсии обогащенных топливом метано-воздушных смесей в синтез-газ (Н2+СО) или для осуществления этого процесса при более низкой температуре;

6. изучить влияние наличия возбужденных молекул 02(<7'Д„) в смесях Н2-воздух и СЫрВоздух на скорость и концентрационные пределы распространения ламинарного пламени;

7. провести сравнение энергетической эффективности метода интенсификации горения, основанного на возбуждении молекул 02 в состояния а и ¿>'Х8+, и метода, основанного на нагреве газа;

Научная новизна

1. Предложена полуэмпирическая методика расчета константы скорости реакции с участием электронно-возбужденных молекул 02(а'Д„) и 02(б'Е„+), в случае, когда константа скорости аналогичной реакции с участием молекулы 02, находящейся в основном электронном состоянии Х3!^", известна.

2. Созданы детальные кинетические модели окисления водородно-воздушных и метано-воздушных смесей при наличии в них электронно-возбужденных молекул 02(аАе) и О2(£'28+). Проведено тестирование разработанных кинетических моделей на имеющихся экспериментальных данных, в том числе и при наличии в смеси возбужденных молекул 02, и дано объяснение наблюдаемых в экспериментах эффектов.

3. Проведено планирование эксперимента по влиянию возбужденных молекул 02 на процесс горения водородно-кислородной смеси, а именно, определены параметры течения смеси Н2-02 в проточном реакторе, при которых можно визуально наблюдать влияние наличия молекул синглет-ного кислорода 02(а'Д8), генерируемого в электрическом разряде, на длину задержки воспламенения.

4. С использованием разработанных моделей проведено исследование влияние возбуждения молекул 02 в состояния а] Ае и на время задержки воспламенения в замкнутом реакторе и на длину задержки воспламенения в сверхзвуковом потоке за наклонной ударной волной в смесях Н2-воздух и СН4-воздух. Показано, что наличие в горючей смеси электронно-возбужденных молекул 02(а'Д„) и 02(й1Е8+) приводит к сокращению задержки воспламенения и к уменьшению температуры воспламенения рассматриваемых смесей. Также продемонстрирована возможность интенсификации процесса конверсии метана в синтез-газ (Н2+СО) при возбуждении молекул 02.

5. Изучено влияние присутствия молекул 02(а'Д8) на скорость и концентрационные пределы распространения ламинарного пламени в смесях Н2-воздух и СН4-воздух различного состава и показана возможность существенного (до 20-70%) увеличения скорости ламинарного пламени даже при небольшом содержании молекул 02(а'Д„) в 02.

6. Исследованы механизмы протекания цепных реакций при горении смесей Н2-воздух и СН4-воздух, содержащих в небольшом количестве электронно-возбужденные молекулы 02(а'Д8) и 02(¿>'Е8+). Показано, что определяющую роль в интенсификации воспламенения и горения в таких смесях играют реакции инициирования и разветвления цепи с участием молекул 02(а'Д3).

7. Установлено, что для рассмотренных в работе горючих смесей метод интенсификации горения, основанный на возбуждении молекул 02 в состояния а'д„ и б'Х„+, намного эффективнее метода, основанного на нагреве смеси, при одинаковой подведенной к газу энергии.

Практическая значимость

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

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

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

На защиту выносятся следующие основные результаты и положения:

1. детальные кинетические модели для описания процессов, протекающих при инициировании воспламенения и горения смесей Н2-воздух и СН4-воздух путем возбуждения молекул 02 в метастабильные состояния а'Д8 и Ъ\+\

1. результаты численного исследования влияния возбуждения молекул 02 в состояния а'Д„ и при импульсном воздействии лазерного излучения на время задержки воспламенения, а также на температуру воспламенения в смесях Н2-воздух и СРЦ-воздух;

3. обоснование возможности существенного сокращения длины зоны индукции в сверхзвуковом потоке за фронтом ударной волны при инициировании горения смесей Н2-воздух и СН4-воздух путем возбуждения электронных состояний молекул 02 резонансным лазерньм излучением;

4. анализ возможности использования метода интенсификации горения, основанного на возбуждении молекул 02 в состояния а1 Д8 и для ускорения конверсии обогащенных топливом метано-воздушных смесей в синтез-газ (Н2+СО) или для осуществления этого процесса при более низкой температуре;

5. результаты численного исследования влияния наличия молекул 02(а'Д8) в смесях Н2-воздух и СН4-воздух на скорость и концентрационные пределы распространения ламинарного пламени;

6. анализ механизмов развития цепного процесса при инициировании горения смесей Н2-воздух и СН4-воздух путем возбуждения молекул 02 в состояния я'Д„ и

7. сравнительный анализ эффективности метода интенсификации горения, основанного на возбуждении молекул 02 в электронные состояния а'Д8 и

и метода, основанного на нагреве смеси.

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

Основные результаты докладывались на 21 российской и международной

конференциях:

1. International Colloquium on Control of Detonation Processes, 4-7 July 2000, Moscow, Russia;

2. VIII Всероссийский съезд по теоретической и прикладной механике, 2329 августа 2001, Пермь, Россия;

3. International Colloquium on Advances in Confined Detonations, 2-5 July 2002, Moscow, Russia;

4. VI International School-Seminar "Nonequilibrium Processes and their Applications", 31 August - 5 September 2002, Minsk, Belarus;

5. 25th International Congress on High-Speed Photography and Photonics, 29 September - 4 October, 2002, Beaune, France;

6. International Symposium on Combustion and Atmospheric Pollution, 8-11 July 2003, St. Petersburg, Russia;

7. 7th International Conference on Laser Ablation, October 5-10, 2003, Crete, Greece;

8. International Colloquium on Application of Detonation for Propulsion, 6-9 July 2004, St. Petersburg, Russia;

9. International Workshop "Nonequilibrium processes in combustion and plasma based technologies", 21-26 August 2004, Minsk, Belarus;

10. 17th International Symposium on Air Breathing Engines, 4-9 September 2005, Munich, Germany;

11. 2nd International Symposium on Nonequilibrium Processes, Combustion and Atmospheric Phenomena, 3-7 October 2005, Sochi, Russia;

12. ICONO/LAT, 28 May -1 June 2007, Minsk, Belarus;

13. 3rd International Symposium on Non-Equilibrium Processes, Plasma, Combustion and Atmospheric Phenomena, 25-29 June 2007, Sochi, Russia;

14. IEEE International Conference on Plasma Science, 15-19 June 2008, Karlsruhe, Germany;

15. 7th International Symposium on Hazards, Prevention, and Mitigation of Industrial Explosions, 7-11 July 2008, St.-Petersburg, Russia;

16. 6th International Colloquium on Pulse and Continuous Detonations, 10-12 November 2008, Moscow, Russia;

17. 6th International Seminar on Flame Structure, Brussels, Belgium, 14-17 September 2008;

18. 19th International Symposium on Plasma Chemistry, Bochum, Germany, 26 -31 July 2009;

19. 22nd International Colloquium on the Dynamics of Explosions and Reactive Systems, 27-31 July 2009, Minsk, Belarus;

20. 4th International Symposium on Non-Equilibrium Processes, Plasma, Combustion, and Atmospheric Phenomena, Sochi, Russia, 5-9 October 2009;

21. III Международная научно-техническая конференция "АВИАДВИГАТЕЛИ XXI ВЕКА", 30 ноября-3 декабря 2010, ЦИАМ, Москва, Россия.

Разработанные в диссертационной работе кинетические модели окисления смесей Н2-воздух и СН4-воздух при наличии в них молекул синглетного кислорода используются различными исследовательскими группами Украины, Беларуси, Франции, США. Реакции с участием молекул 02(a'Ag) и соответствующие им константы скорости включены в кинетическую базу данных

[9].

Результаты работы обсуждались на семинаре по механике сплошных сред под руководством академика РАН А.Г. Куликовского, проф. В.П. Кар-ликова и чл.-корр. РАН О.Э. Мельника в НИИ механики МГУ им. М.В. Ломоносова в 2011 г.

Публикации и личный вклад автора Основное содержание и результаты диссертационного исследования изложены в 29 работах [А1-А29], в том числе в 23 статьях [А1-А23] в рекомендованных ВАК журналах. Во всех работах соискателю принадлежит участие в постановке задачи, численном моделировании и анализе результатов. Все положения, выносимые на защиту, получены лично соискателем.

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

Диссертация состоит из введения, четырех глав и заключения. Работа изложена на 176 листах, содержит 65 рисунков, 12 таблиц и список библиографических ссылок из 188 источников.

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

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

Первая глава содержит обзор работ, посвященных изучению элементарных химических процессов с участием молекул синглетного кислорода 02(a'Ag) и а также экспериментальным и теоретическим исследова-

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

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

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

Впервые влияние предварительной активации кислорода тлеющим электрическим разрядом на скорость распространения пламени в бедной водо-родно-воздушной смеси при низком давлении было обнаружено в [10]. Авторы связывали этот эффект с наличием в продуктах разряда молекул 02(а'Д„). И лишь в последние годы были проведены эксперименты [11, А21], в которых непосредственно изучалось влияние ОгСа'До) на процесс горения смесей Н2-02-Аг(Не) и С2Н4-С>2-Аг(Не). Во всех работах наблюдалось пром.о-тирующее действие молекул 02(а'Дв) на процесс воспламенения и горения.

Первая кинетическая модель, призванная объяснить экспериментально наблюдаемое в [10] увеличение скорости пламени в водородно-воздушной смеси, была предложена в [2]. Расчеты относительного увеличения скорости пламени при добавке молекул 02(а'Д8) качественно соответствовали наблюдаемой в эксперименте тенденции. Обнаружено, что наибольший вклад в увеличение скорости пламени дает реакция Н+02(а'Д8)=0Н+0. Из-за недостатка данных на тот период эта модель не включала некоторые очень важные процессы: (1) электронно-поступательную релаксацию молекул 02(а1Д6); (2) диссоциацию 02(а'Д„); (3) образование молекул 0:(/>'Х„~) в так называемом пулинг-процессе; а также (4) химические реакции с молекулами 0:(й'Е„+).

В 2001 году совместно с научным руководителем была предложена идея интенсификации цепных реакций в горючих смесях путем возбуждения молекул 02 в метастабильные состояния а1 А„ и резонансным лазерным излучением [АЗ]. Для теоретического обоснования этой идеи требовались соответствующие кинетические модели. Однако, модель [2] для водородно-воздушной смеси требовала пересмотра на основе современных данных по кинетике элементарных процессов. Необходимым было тестирование модели на экспериментальных данных. Кинетические механизмы, описывающие горение даже самых простых углеводородов, например метана, при наличии в смеси синглетного кислорода, в литературе вообще отсутствовали.

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

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

Возбуждение колебательных или электронных степеней свободы реагирующих молекул приводит к снижению барьера эндотермических реакций [12]. Пусть константа скорости обменной реакции АВ+С=А+ВС имеет Ар-рениусовский вид к{Т) = Л-Т" ■ ехр(- Еа¡Т), где Т- температура газа, Еа -энергия активации этой реакции. Согласно [12] поверхность потенциальной энергии для прямой реакции можно представить в виде (рис. 1)

их=АН + Еа ехр(г/гх),

где АН - тепловой эффект реакции. Потенциальная поверхность для обратной реакции при этом будет иметь следующий вид

и 2 = [аН + Еа )ехр(- г! г2). Здесь г\ и г2 - радиусы действия обменных сил для реагентов и продуктов соответственно. При протекании реакции АВ(е)+С=А+ВС с участием молекулы АВ(е), возбужденной в электронное состояние 'е' с энергией возбуждения Ее, поверхность потенциальной энергии сместится вверх на величину Ее по отношению к соответствующей поверхности для реакции с участием невозбужденной молекулы

= АН + Ее+ Еа ехр(г/7-,). (1)

Будем считать, что форма потенциальной поверхности и[ реакции с участием электронно-возбужденной молекулы аналогична форме потенциальной поверхности 1/\ для молекулы в основном электронном состоянии (в общем случае это предположение может не выполняться). Для многих реакций [12] Г]-Г2. В этом случае в точке пересечения потенциальных поверхностей и' =и2, и справедливо соотношение

Еа12+(АН + Ее)1-(АН + Еа) = 0, (2)

Рис. 1. Энергетическая диаграмма экзотермической обменной реакции АВ+С= А+ВС при участии невозбужденной АВ и возбужденной АВ(е) молекул.

где t = ехр(/-/г,). Подставляя решение уравнения (2) в (1), для энергии активации Ееа реакции АВ(е)+С=А+ВС получаем следующее выражение

К =\ (У И + Ее)2+4Еа(АН + Еа)-(АН + Ее)}.

Константу скорости реакции с участием молекулы АВ(е) записываем в виде ке(Т) = А -Т" ехр(- El/Т), где Ееа<Еа.

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

Например, константа скорости реакции инициирования цепи в метано-воздушной смеси СЬЦ+С^а'Дп^СНз+НОз рассчитывалась с помощью "bond-energy-bond-order" [13] и ab initio [14] методов. В частности, в [14] сделан вывод, что энергия активации реакции с участием синглетного кислорода меньше, чем при взаимодействии метана с кислородом, находящемся в основном электронном состоянии. Расчет константы скорости реакции по предложенной в диссертационной работе методике дает выражение, которое неплохо согласуется с полученными в [13, 14] значениями (рис. 2).

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

Третья глава содержит описание разработанной в данной работе кинетической модели воспламенения и горения смеси Н2-воздух, содержащей электронно-возбужденные молекулы 02, а также сравнение результатов расчета по этой модели с экспериментальными данными как для обычной смеси Н2-воздух, так и для водородно-воздушной смеси, содержащей молекулы 02(а'Д8) и 02(i1Sg+). Здесь же проводится исследование влияния возбуждения

к, см /(мольх) 1.Е+12

СН4+02(й Дд)=СН3+Н02 [14]

Рис. 2. Константы скорости реакций СН4+02=СН3+Н02 и СН4+02(а1Д8)= СН3+Н02 (сплошная и штриховые линии), рассчитанные различными методами.

молекул 02 в состояния <з'Лц и Ь]1.Ё+ на время индукции и температуру воспламенения смеси в адиабатическом реакторе, а также на длину задержки воспламенения за ударной волной и на скорость распространения ламинарного пламени. Проводится анализ механизмов, ответственных за интенсификацию процесса горения, и сравнительный анализ эффективности метода интенсификации горения, основанного на возбуждении молекул 02 в состояния а]А„ и и обычного теплового метода инициирования горения.

Представлена детальная кинетическая модель окисления смеси Н2-воздух при неравновесном возбуждении молекул 02 в состояния а'Д„ и ¿'хД включающая химические реакции, в том числе с электронно-возбужденными молекулами 02(ахА„) и 02(Ь^1,„+) и атомами О('О), а также процессы электронно-электронного обмена и электронно-поступательной релаксации. При построении кинетической модели для констант скоростей реакций, где это возможно, использовались экспериментальные данные. В отсутствии таких данных, константы скоростей реакций с участием возбужденных компонентов определялись по методике, предложенной в главе 2. Константы скоростей обратных реакций вычислялись через константу равновесия на основе принципа детального равновесия. Для определения константы равновесия были рассчитаны статистические суммы для электронно-возбужденных компонентов ОгСа'Лз), 02(б'Х„т) и О('О) и на их основе определены приведенные энергии Гиббса, входящие в выражение для константы равновесия.

Тестирование модели проводилось на доступных экспериментальных данных по времени задержки воспламенения, скорости распространения ламинарного пламени и эволюции концентрации различных компонентов за фронтом отраженной ударной волны и в проточном реакторе для обычной смеси Н2-02(воздух), а также на имеющихся экспериментальных данных в смеси Н2-02, содержащей молекулы 02(а'Л„) [10, 15, А21]. Тестирование показало, что модель правильно описывает все наблюдаемые в экспериментах особенности процессов окисления, воспламенения и горения как в обычных смесях Н2-02(воздух), так и в смесях, содержащих молекулы 02(а'Д„).

Так в [10] исследовалось влияние образующихся в тлеющем разряде молекул 02(а'Д„) на скорость распространения ламинарного пламени (/„ в очень бедных смесях Н2-02 при низком давлении. Расчет, также как и эксперимент, показал увеличение скорости распространения пламени (таблица 1). Причем расчетное значения ип увеличивается сильнее для смеси с меньшей концентрацией Н2, что наблюдалось и в эксперименте. Отличие между рассчитанны-

Таблица 1. Скорость ламинарного пламени для различных смесей Н2-02.

[/„, см/с £//*, см/с и„"'Ж

мольная доля Н2 расчет расчет расчет эксперимент [10]

10.5% 17.07 36.51 2.14 1.66

8.5% 6.07 16.91 2.78 1.92

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

Анализ результатов работы [15], в которой изучалось тушение полученных в разряде молекул 02(a'As) различными компонентами, в том числе и молекулами Н2, показал, что авторы [15] неправильно интерпретировали результаты эксперимента, предположив, что уменьшение концентрации молекул 02(a'Ag) за фронтом ударной волны обусловлено лишь процессом тушения 02(а'Л8). Это привело их к нефизичному значению константы скорости этого процесса, поскольку не были учтены химические реакции с участием возбужденных молекул 02(л'Д8), которые могут протекать в смеси Н2-02 при 7>600 К. Расчет по представленной в диссертационной работе кинетической модели с общепринятым значением константы тушения 02(а'Ag) молекулами Н2 позволил описать динамику уменьшения концентрации 02(а'Д8), полученную в эксперименте при всех значениях Г за фронтом ударной волны.

Изменение длины задержки воспламенения при введении в водородно-воздушную смесь молекул 02(а'Д„) впервые было зарегистрировано в [А21]. Для экспериментального исследования влияния синглетного кислорода на длину зоны индукции был выбран метод, основанный на организации стационарного горения в проточном реакторе низкого давления, на вход которого раздельно подаются холодные газообразные водород и кислород. При этом

молекулярный кислород мог подвергаться действию разряда, который генерирует возбужденные молекулы 02(a'Ag). Предварительно было проведено планирование эксперимента и показано, что для обеспечения вос-730 к пламенения обычной водородно-

кислородной смеси на длине реактора 7Ю к 675 к 70 см при скорости потока 10-20 м/с, g.......о стенки реактора надо нагревать до

9 775 К 50'

"40

30 20 -\

О

10 I 700-1000 К. При наличии же

о 0.01 0.02 0.03 0.04 0.05 у02(а1) 5°/о 02(а1Д8) в потоке кислорода мож-

. „ . но ожидать пятикратного уменыпе-Рис. 3. Длина зоны индукции ь.„ в смеси

Н2/О2=5/2сР0=ЮТорр,и0=17м/свза- ния длины зоны воспламенения,

висимости от мольной доли молекул Проведенные эксперименты полно-

02(а'ле) в потоке молекулярного кисло- стыо подтвердили результаты расче-

рода. »-эксперимент [А21], линия- та. На рис. 3 показана длина зоны

расчет. Над расчетными точками даны индукции смеси Н2-02 как функция

значения температуры газа перед сече- мольной доли молекул 02(аЛЕ) в мо-

нием воспламенения. лекулярном кислороде, полученная в

эксперименте [А21] и рассчитанная по модели данной работы. Следует отметить, что вследствие неоднородности изменения температуры газа по продольной координате реактора, воспламенение при наличии молекул 02(сг'Д„) в смеси происходило не только на более коротком расстоянии от начала реактора, но и при более низкой температуре.

Проведен анализ процесса воспламенения смеси Н2-воздух при воздействии импульса резонансного лазерного излучения на покоящуюся среду. Для определения параметров газа в центральной области лазерного пучка, на основе полной системы уравнений, включающих процессы макро- и микропереноса, действующих в поле импульса лазерного излучения, выделены характерные времена определяющих процессов и проведена их оценка при исследуемых в работе параметрах среды Р0=104Па, Г0=500-Ю00 К, интенсивности /=1-10 кВт/см2 и длительности т^Ю"4 с импульса воздействующего излучения в области облучения Ла=1 см. Будем рассматривать режимы, когда характерное время воспламенения (период индукции) т¡„ удовлетворяет условию

т0«Тр~Т/ < Т/СА~Т/Я«(ТС/, Тту, Ч, Гр) <т„ (3)

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

сЬ

индуцированных переходов, т, - характерное время химическои реакции, приводящей к образованию 1-ого компонента, тд, и - времена многокомпонентной и термодиффузии, т> - время теплопроводности, хр - время изменения состояния среды при воздействии стрикционной силы, т5 - время спонтанной эмиссии. Из неравенства (3) следует, что на интервале [0, т,„] процессами макропереноса и спонтанным излучением можно пренебречь и считать, что заметного движения газа поперек пучка не происходит. В этом случае определяющая система уравнений имеет следующий вид

+ , к Ь ^ я сА +«лМУ/0 ,

ч Т,-

У? = 0,

пдЕ-1гТ\ь Р у К « ТЛ Р ^ К тР ^ VI ^ Я I сут10 ) п , СутУ0 т,-

где N1 - концентрация /-ого компонента смеси, /г0, - энтальпия его образования, N = Nj , у¡=А',/ЛГ, Бц - скорость образования /-ого компонента в д-ой реакции, ^ , ку - коэффициент поглощения лазерного излучения,

1п~ число квантов, теряемых (приобретаемых) молекулой /-ого сорта при индуцированных переходах.

Проведенный численный анализ показал, что возбуждение молекул 02 в состояния а'Л„ и Ь1Т,„+ лазерным излучением приводит к значительному

уменьшению времени индукции т,„ и температуры воспламенения смеси Н2-воздух (рис. 4). При одном и том же значении вложенной энергии Е„ сокращение времени индукции при воздействии обоих типов излучения практически одинаково. Однако с практической точки зрения легче организовать возбуждение молекул 02 в состояние т.к. при одинаковых параметрах среды скорость индуцированных переходов X —существенно (в 75 раз) больше, чем скорость переходов Х%~->а1 Д8.

Сокращение времени задержки воспламенения при возбуждении молекул 02 связано с появлением новых каналов образования активных атомов О, Н и радикалов ОН: Н2+0= Н+Н02, Н+02(я'Дг) = ОН+О. Константы скорости этих процессов существенно больше, чем константы скорости соответствующих реакций с участием невозбужденных молекул 02. Поэтому атомы Н, О и радикалы ОН образуются раньше и их концентрации больше, по сравнению со случаем отсутствия облучения.

Аналогичный эффект в смеси Н2-воздух наблюдается при инициировании процесса горения в сверхзвуковом потоке за фронтом наклонной ударной волны. Это иллюстрирует рис. 5. При числе Маха набегающего потока М0=8 и отсутствии излучения (за фронтом ударной волны температура 7'1=924 К и давление Рх=1.3x104 Па) длина зоны индукции составляет 7 м, т.е. можно считать, что смесь не воспламеняется. При возбуждении молекул 02 перед фронтом ударной волны в состояние излучением с к¡=762 нм

воспламенение реализуется на длине 1 м при подведенной энергии £¡=0.1 эВ/(молекулу 02). Если эта энергия лазерного излучения идет на нагрев газа, то температура перед фронтом повышается до 369 К, а за фронтом ударной волны достигает 71=980 К (рис. 6). И хотя длина зоны воспламенения сокращается и составляет уже 4.2 м, это значение все же существенно больше, чем при возбуждении молекул 02 в состояние 02(б'2е+).

Длина зоны индукции зависит также от давления газа. Для смеси Н2-воздух для каждого £, существует две области изменения Р0. В одной Ьы уменьшается с ростом Ро, а в другой - увеличивается. Граничное значение

Рис. 4. Зависимость т,„ от начальной температуры Т0 стехиометрической смеси II;-воздух при отсутствии излучения (1) и при воздействии излучения с Я/=1.268 мкм и 762 нм (сплошные и пунктирные линии) при £,=0.01; 0.03; 0.06; 0.12 эВ/(молекулу 02) (2-5), х= 10"4 с. Штриховая линия -т,„(7о), когда энергия £,=0.12 эВ/(молекулу 02) идет на нагрев газа. Рь= 10 кПа.

Роь, разделяющее эти области, растет с увеличением Существование Роь и увеличение Ь,„ при Р<у>Роъ обусловлено тем, что при увеличении Р0 происходит ускорение трехчастичной реакции рекомбинации Н+02+М=Н02+М, а последующая реакция 2Н02=Н202+02 приводит к образованию химически инертных молекул Н202 и уменьшению концентрации атомов О, Н и радикалов ОН. Чем выше Е„ а, следовательно, и концентрация 02(я'Д8) в смеси, тем при больших значениях Р0 эти процессы начинают доминировать. Таким образом, возбуждение молекул Ог в синглетные состояния и ¿>'28+ позволяет существенно расширить диапазон давлений, в котором воспламенение реализуется на приемлемых расстояниях (¿,„<1 м).

Важной характеристикой процесса дефлаграционного горения является скорость распространения ламинарного пламени и„, которая, по сути, определяет скорость выгорания смеси. Поэтому в работе был проведен специальный анализ возможности увеличения скорости распространения пламени в смеси Н2-воздух путем возбуждения молекул 02 в синглетное состояние а'Д„. Проведенный численный анализ показал, что присутствие молекул 02(д1Д8) в смеси Н2-воздух приводит к заметному увеличению 11„ даже при умеренном содержании молекул 02(а'Дв) (рис.7). Наибольший относительный рост £/„ (до 55%) наблюдается в обедненных топливом смесях.

Рис. 5. Зависимость ¿,„(М0) при горении смеси Н2-воздух (ф=1) за фронтом УВ при 7"0=300 К, ДНО3 Па, (5=25°, при отсутствии излучения (1) и при воздействии излучения с £,=0.01; 0.03; 0.05 и 0.1 эВ/(молекулу 02) (2-5). Пунктирная линия соответствует 1,>,(Мо), когда энергия £,=0.1 эВ/(молекулу 02) идет на нагрев газа.

Рис. 6. Изменение температуры за фронтом ударной волны в смеси Н2-воздух (ф=1) при £5=0 (1), когда энергия £ =0.1 эВ/(молекулу 02) идет на нагрев (2) или возбуждение молекул 02 (3). Го=300 К, Р0=Ю3 Па, М0=8, р=25°.

Скорость пламени определяется двумя процессами: образованием активных атомов и радикалов, определяющих скорость цепного процесса, и диффузией активных атомов и радикалов из горячей зоны пламени в холодную. Скорость реакции разветвления цепи Н+02(<з' Д8)=0+0Н с участием молекул ОгСа'Дц) значительно выше, чем скорость реакции, в которой 02 находится в основном электронном состоянии. Поэтому при возбуждении молекул 02 увеличивается концентрация активных атомов О и Н и радикалов ОН - носителей цепного механизма в горячей зоне пламени. При этом увеличивается и скорость диффузии этих атомов и радикалов в холодную зону, где они инициируют протекание цепных реакций. Это и приводит к увеличению скорости распространения пламени.

В четвертой главе приводится описание кинетической модели окисления метана в воздухе при наличии в смеси молекул 02(а'Д3) и 02(&'Е8+) и атомов О('Б), проводится тестирование модели на экспериментальных данных и анализ влияния присутствия в смеси возбужденных молекул 02 на динамику воспламенения метано-воздущных смесей в замкнутом реакторе, в сверхзвуковом потоке за фронтом ударной волны, а также на скорость распространения ламинарного пламени, определяются кинетические механизмы, ответственные за интенсификацию процессов воспламенения и горения.

Кинетическая модель окисления метана в воздухе основана на реакционном механизме [16]. Для определения минимального набора процессов, позволяющего правильно описать динамику воспламенения и тепловыделение, была проведена редукция полной кинетической схемы и показано, что использование модели, включающей 29 компонентов, позволяет получить правильные значения периода индукции т,„ и конечной температуры газа Те при воспламенении смеси СН4-02. Эта схема была дополнена реакциями с участием молекул 02(а1Д£) и О2(Ь1И&+) и атомов О('В). Из-за отсутствия экспериментальных данных по константам скоростей реакций с возбужденными компонентами, эти константы были вычислены по приведенной в главе 2 методике. И только для небольшого числа процессов электронно-поступательной релаксации были найдены экспериментальные данные. Разработанная модель окисления метана включает как составную часть механизм окисления водорода.

О -I-1-1-1-1-

О 1 2 3 4 ф

Рис. 7. Зависимость £/„(ф) для смеси Н2-

воздух с мольной долей у® , =0;

02(а Л),)

5 и 10% у°02 при 7о=300 К, /-о=105 Па.

Метан воспламеняется намного хуже водорода и для его воспламенения требуется более высокая температура. Однако возбуждение молекул 02 в состояния 02(а,1А8) и О2(^>1Е8+) позволяет существенно снизить температуру воспламенения В первом приближении Т,^ определяется из уравнения т,„(7о, Р0, Е^=Хо,(Т0, Р0), где тс, - время диффузии активного компонента, определяющего развитие цепного механизма (для смеси СН4-воздух ¡=0). Численный анализ показал (рис. 8), что при Р0= 104 Па Г/^,уменьшается с -1350 К при Е5=0 до ~850К при £<=0.12 эВ/(молскулу 02). Если же воспламенение реализуется, то возбуждение позволяет его заметно ускорить. Такой метод интенсификации горения более выгоден по сравнению с нагревом газа, так как позволяет сильнее сократить т/и при одинаковой вложенной энергии. В случае же обеспечения одинакового времени индукции метод, основанный на возбуждении молекул 02 в состояние Ь'1.,/, требует в 7 раз меньше энергии, чем метод, основанный на нагреве газа. При этом инициирование горения путем возбуждения молекул 02 с /./ =762 нм обеспечивает более высокую (~2 раза) долю химической энергии реагентов АНс/„ выделяющейся в тепло

Рис. 8. Время индукции т,„ и диффузии х0 атомов О (сплошные и пунктирная линии) в смеси СН4-воздух (Ро=\04 Па, ф=1) при отсутствии излучения (1) и при воздействии излучения с /,/=762 нм и £,=0.01; 0.03; 0.06; 0.12 эВ/(молекулу 02) (2-5). Штриховая линия - зависимость т,„(7,0), когда вся энергия £,=0.12 эВ/(молекулу 02) идет на нагрев газа.

и Т/^г - значения температуры воспламенения при возбуждении молекул 02 излучением с £,=0.12 эВ/(молекулу 02) и при нагреве.

Таблица 2. Энергия, подводимая к газу Е„ температура продуктов сгорания Те и химическая энергия реагентов, выделяющаяся в тепло, ДД^ при различных способах инициирования горения в смеси СН4-воздух при Г0=Ю00 К, Ро=Ю4 Па, ф=1 для дос^ тижения т,„=0.013 с._

Процесс эВ/(молек. 02) Е„ Дж/см Т/я« С Те, К Д#сь Дж/см

Нагрев газа 0.834 0.0184 0.013 2663 0.0294

Возбуждение молекул 02 0.120 0.0026 0.013 2574 0.0565

(см. таблицу 2).

В настоящее время рассматривается возможность использования в различных устройствах синтетических топлив, в том числе синтез-газа, который можно получить при частичном окислении углеводородов. Проведенные расчеты показали, что возбуждение молекул 02 в состояния а'Д„ или Ьх1,8+ позволяет, во-первых, ускорить процесс частичного окисления богатой метано-воздушной смеси, а, во-вторых, проводить этот процесс при более низкой температуре. При энергии £,=0.12 эВ/(молекулу 02) процесс конверсии смеси СН4-воздух с ф=3 может быть осуществлен при 7о=900 К, что на 350 К меньше, чем при отсутствии излучения. Подвод энергии в поступательные степени свободы, т.е. нагрев смеси, позволяет понизить температуру только на -60 К. Воздействие лазерного излучения не приводит к изменению концентрации Н2 и СО в продуктах конверсии, что видно из рис. 9, иллюстрирующем влияние возбуждения молекул 02 на динамику окисления обогащенной топливом метано-воздушной смеси в адиабатическом реакторе. С другой стороны, можно говорить о том, что для заданного времени пребывания газа в реакторе тт (или длины проточного реактора) возбуждение молекул позволяет существенно увеличить выход Н2 и СО при конверсии метана. Например, для тге5=0.1 с мольная доля Н2 и СО при частичном окислении смеси СН4-воздух с £,==0.06 эВ/(молекулу 02) на выходе из реактора составляет соответственно 24% и 12%, а при £5=0.12 эВ/(молекулу 02) их концентрации равны 29 и 16%. Без возбуждения процесс конверсии метана при данных начальных параметрах для такого тге5 вообще не реализуется.

Проведенный численный анализ показал, что возбуждение молекул 02 в метано-воздушной смеси также эффективно с точки зрения сокращения длины зоны индукции в сверхзвуковом потоке за фронтом ударной волной. Уменьшение длины зоны индукции составляет 20 и более раз при вложенной энергии Е$=0.1 эВ/(молекула 02) по сравнению со случаем отсутствия облучения. Тепловой метод воздействия лазерного излучения дает выигрыш только в 2 - 2.5 раза, несмотря на более высокую температуру газа за фронтом ударной волны.

Присутствие синглетного кислорода увеличивает скорость ламинарного пламени в смеси СН4-воздух даже больше, чем в водородно-воздушной смеси. Так, для бедной смеси (ф=0.45) это увеличение максимально и составляет 67% (против 55% в смеси Н2-

1.Е-05 1.Е-03 1.Е-01 Рис. 9. Эволюция мольных долей СН4, н2 и СО при объемной реакции в смеси СН4-воздух с ф=3, Г0=Ю00 К, Ро=Ю5 Па при воздействии лазерного излучения с Я,/ =762.346 нм и = 0, 0.3, 0.6, 0.12 эВ/(молекулу 02) (1-4).

воздух). Аналогичная тенденция наблюдается и в богатой смеси (см. табл. 3). Также при добавлении синглетного кислорода увеличивается конечная температура Тс и химическая энергия, выделившаяся в тепло при горении, АНС},.

Анализ показал, что при наличии молекул 02(£г'Да) в смеси СН4-воздух за интенсификацию горения ответственны реакции инициирования СН4+02(а'Д8)=СН3+Н02 и разветвления цепи СН3+02(а1Д8)=СН20+0Н и СНз+02(а'Л8)= СН30+0, константы скорости которых существенно больше, чем в соответствующих реакциях с 02 в основном состоянии. Эти реакции приводят к образованию активных радикалов СН3, ОН и атомов О и, в итоге, к ускорению горения.

Таблица 3. Параметры пламени при Г0=450 К и Рг,-1 атм.

СИ,- воздух Н2- воздух

0% 02(а'Д„) 10%02(а'Дв) 0% 02(а'Д„) Ю%02(я'Дв)

<¡>=0.45 1/„, см/с 10.9 18.1 115.1 178.9

ДС//С/„о, % - 67 - . 55

Те, К 1496 1553 1617 1679

ДНа, кДж/г 1.27 1.34 1.46 1.55

ф=1 1/„, см/с 70.7 84.9 409.2 486.6

Д Ши„0, % - 20 - 19

Те, К 2276 2304 2430 2453

АНС1,, кДж/г 2.52 2.55 2.97 3.0

ф=1.9 £/„, см/с 10.9 15.2 515.7 591.0

АШим, % - 38 - 15

Те, К 1770 1807 2224 2253

ДНс11, кДж/г 2.03 2.09 3.24 3.29

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

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

2. Созданы детальные кинетические модели окисления смесей Н2-воздух и СН4-воздух при наличии в них электронно-возбужденных молекул 02(а'Д„) и 02(б'Е8+), а также атомов О('Б). Проведено тестирование этих кинетических моделей на известных экспериментальных данных по времени задержки воспламенения, скорости распространения ламинарного

пламени и изменению концентрации компонентов в проточном реакторе и за фронтом ударной волны, в том числе и для случаев, когда в смеси Н2-02 присутствуют электронно-возбужденные молекулы 02(о'Д8), и продемонстрирована хорошая предсказательная способность разработанных моделей. Дано объяснение наблюдаемых в экспериментах эффектов по влиянию возбуждения молекул кислорода на динамику воспламенения и горения.

3. На основе разработанной кинетической модели проведено планирование эксперимента, а именно, определены параметры потока в проточном реакторе, при которых можно четко наблюдать влияние присутствия молекул 02(а'Д8) в потоке на воспламенение смеси Н2-02. Выполненный по разработанной кинетической модели расчет длины зоны индукции для смеси Н2-02 в зависимости от концентрации молекул 02(а'А„) показал хорошее соответствие с экспериментальными данными.

4. Проведен анализ иерархии характерных времен процессов микро- и макропереноса и получена система уравнений, описывающая изменение параметров в зоне облучения смесей Н2-воздух и СН4-воздух при импульсном воздействии лазерного излучения с длиной волны Я,/=1.268 мкм или 762 нм, приводящем к возбуждению молекул 02 в электронные состояния а1 Д8 или На основе численного моделирования показано, что возбуждение молекул 02 в состояния аА„ и позволяет существенно (до 10 и более раз, в зависимости от начальной температуры и вложенной энергии) уменьшить время задержки воспламенения, а также понизить температуру воспламенения горючих смесей. При вложенной энергии £¡=0.12 эВ/(молекулу 02) температура воспламенения смеси Н2-воздух может быть уменьшена с -800 до -500 К при давлении 1 атм, а для смеси СН4-воздух - с -1350 до -850 К при давлении 0.1 атм.

5. Исследован процесс воспламенения и горения гомогенных смесей Н2-воздух и СН4-воздух в сверхзвуковом потоке за фронтом наклонной ударной волны при возбуждении электронных состояний молекулярного кислорода лазерным излучением. Показано, что даже при малой подведенной к газу энергии £5=0.1 эВ/(молекула 02) лазерное излучения позволяет инициировать горение на расстоянии менее 1 м от фронта относительно слабых ударных волн, когда температура газа не превышает 900 К для смеси Н2-воздух и 1200 К для смеси СН4-воздух.

6. Изучена возможность интенсификации процесса конверсии обогащенной топливом смеси СНгвоздух в синтез-газ (Н2+СО) при возбуждении молекул 02 в состояния а1 Д8 или Показано, что в этом случае конверсия может быть осуществлена при достаточно низкой начальной температуре (-900 К) и атмосферном давлении. При заданной же температуре присутствие даже незначительного количества молекул 02(а'Дв)

(~3% от концентрации 02) в исходной смеси приводит к существенному сокращению времени окисления метана и тем самым, позволяет увеличить выход Н2 и СО в продуктах конверсии в реакторе с заданным временем пребывания по сравнению с обычной метано-воздушной смесью.

7. Исследовано влияние присутствия в небольшом количестве (5-10% от концентрации 02) молекул синглетного кислорода на скорость распространения ламинарного пламени в смесях Н2-воздух и CHt-воздух. На основе численного моделирования установлено, что относительное увеличение скорости пламени максимально для обедненной топливом смеси: для смесей CHt-воздух и Н2-воздух при ф=0.45, Г0=450 К, 7° , =0.1 Yn оно составляет-67% и 55%, соответственно, для обога-

02(а Ag) Ю2

щенных топливом смесей относительный рост U„ несколько меньше - 38% и 15%, соответственно. Установлено, что основная причина роста скорости пламени - увеличение концентрации активных атомов и радикалов в горячей зоне пламени в результате ускорения реакций разветвления цепи с участием молекул 02(a'Ag) и, как следствие, увеличение скорости их диффузии из горячей области пламени в холодную предпламенную зону. Присутствие молекул 02(a'Ag) в горючей смеси расширяет бедный и богатый концентрационные пределы распространения пламени.

8. Проведен анализ кинетических механизмов, ответственных за интенсификацию воспламенения смесей Н2-воздух и СН4-воздух при наличии в них молекул 02(a1Ag) и О2(6'Sg+). Показано, что сокращение периода индукции и уменьшение температуры воспламенения обусловлены увеличением скорости образования активных атомов О, Н и радикалов ОН, СН3, СН20, носителей цепного механизма горения водородно-воздушных и ме-тано-воздушных смесей, в реакциях инициирования и разветвления цепи с участием молекул 02(а'Д„) и 02(6'Zg+).

9. Выполнено сравнение эффективности метода инициирования горения, основанного на возбуждении молекул 02 в состояния а]А., или bxZg+, и обычного метода, основанного на нагреве смеси. Показано, что с точки зрения уменьшения температуры воспламенения и сокращения периода индукции, лазерно-индуцированное возбуждение молекул 02 намного эффективнее простого нагрева среды. Так, при поглощенной энергии излучения £/=0.12 эВ/(молекулу 02) температуры воспламенения смеси Н2-воздух уменьшается на 300 К при возбуждении молекул 02 и лишь на 100 К при нагреве смеси. Для смеси СН4-воздух соответствующее уменьшение равно 500 и 100 К. Время индукции при возбуждении молекул 02 сокращается до 10 и более раз (в зависимости от начальной температуры смеси и вложенной энергии) по сравнению с нагревом газа.

ЦИТИРОВАННАЯ ЛИТЕРАТУРА

1. R.C. Brown. A theoretical study of vibrationally induced reactions in combustion processes // Combusion and Flame. - 1985. - Vol. 62, no. 1. - P. 1-12.

2. В.Я. Басевнч, A.A. Беляев. Расчет увеличения скорости водородно-кислородного пламени при добавках синглетного кислорода // Химическая физика. - 1989. -Т. 8, № 8. - С. 1124-1127.

3. A.M. Старик, Н.Г. Даутов. О возможности ускорения горения смеси Н2+О2 при возбуждении колебательных степеней свободы молекул Н2 или 02IIДоклады Академии наук. - 1994. - Т. 336, № 5. - С. 617-622.

4. G.C. Light. The effect of vibrational excitation on the reaction of О(ЗР) with H2 and the distribution of vibrational energy in the product OH // J. Chem. Phys. - 1978. - Vol. 68, no. 6. -P. 2831-2843.

5. В.Я. Басевич, В.И. Веденеев. Константа скорости реакции Н+0:('Дя)= ОН+О // Химическая физика. - 1985. - Т. 4,№8.-С. 1102-1106.

6. K.F. Pliavaka, S. V. Gorbatov, S. V. Shushkov, et al. Singlet oxygen production in electrical non-self-sustained HV pulsed + DC cross discharge at atmospheric pressure with application to plasma assisted combustion technologies // Nonequilibrium processes in combustion and plasma based technologies / Ed. By S.A. Zhdanok. - Minsk: HMTI Publ., 2006. - P. 186-191.

7. A. Hicks, S. Norberg, P. Shawcross, et al. Singlet oxygen generation in a high pressure non-self-sustained electric discharge // J. Phys. D: Appl. Phys. -2005. - Vol. 38. - P. 3812-3824.

8. IV. C. Eisenberg, A. Snelson, R. Butler, et al. Gas phase generation of 02(+Ag) at atmospheric pressure by direct laser excitation II J. Photochemistry. - 1984. - Vol. 25, no. 2-4. - P. 439-448.

9. L.B. Jbraguimova, G.D. Smekhov, O.P. Shatalov. Recommended rate constants of chemical reactions in an H2-02 gas mixture with electronically excited species 02('Д), O('D), OH(2£) involved / Institute of Mechanics of Lomonosov Moscow State University, M., 2003. http://www.chemphys.edu.ru/me-dia/files/-2003-01-20-001 .pdf.

10. В.Я. Басевич, C.M. Когарко. О механизме влияния продуктов тлеющего разряда на скорость водородно-кислородных пламен в условиях полуострова воспламенения // Кинетика и катализ. - 1966. - Т. 7, № 3. - С. 393-398.

11. Т. Ombrello, S.H. Won, Y. Ju, S. Williams. Flame propagation enhancement by plasma excitation of oxygen. Part II: Effects of 02(a'Ag) // Combustion and Flame. -2010.-Vol. 157,no. 10.-P. 1916-1928.

12. В.Д. Русанов, A.A. Фридман. Физика химически активной плазмы. - М.: Наука, 1984.-416 с.

13. S.W. Mayer, L. Schieler. Activation energies and rate constants computed for reactions of oxygen with hydrocarbons // J. Phys. Chem. - 1968. - Vol. 72, no. 7.-P. 2628-2631.

14. A. Starik, A. Sharipov. Theoretical analysis of reaction kinetics with singlet oxygen molecules // Phys. Chem. Chem. Phys. - 2011. - Vol. 13. - P. 1642416436.

15. P. Borrell, D.S. Richards. Quenching of singlet molecular oxygen, 02(й'а„) and 02(/)'S„+), by H2, D2, HC1 and HBr //,/. Chem. Soc., Faraday Trans. 2. -1989. -Vof. 85, no. 9. - P. 1401-1411.

16. Н.Г. Даутов, A.M. Старик. К вопросу о выборе кинетической схемы при описании объемной реакции метана с воздухом // Кинетика и катализ. -1997. -Т. 38, № 2. - С. 207-230.

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

Al. A.M. Старик, Н.С. Титова. О возможности уменьшения порога воспламенения горючих смесей путем селективного возбуждения молекулярных колебаний исходных реагентов. // Доклады Академии наук. - 2000. - Т. 370, № 1. - С. 38-42.

А2. A.M. Старик, Н.С.Титова. О кинетических механизмах инициирования горения водородо-воздушных смесей в сверхзвуковом потоке за ударной волной при возбуждении молекулярных колебаний исходных реагентов II Журнал технической физики. - 2001. - Т. 71, вып. 8. - С. 1-12.

A3. A.M. Старик, Н.С. Титова. Инициирование горения и детонации в горючих смесях при возбуждении молекулярного кислорода в состояние 02(д'А8) И Химическая физика. - 2001. - Т. 20, № 5. - С. 1725.

А4. A.M. Старик, Н.С. Титова. Низкотемпературное инициирование детонационного горения газовых смесей в сверхзвуковом потоке при возбуждении молекулярного кислорода в состояние 02(r/'Ag) И Доклады Академии наук. - 2001. - Т. 380, № 3. - С. 332-337.

А5. А.М. Starik, N.S. Titova. Initiation of combustion and detonation by laser induced electronical excitation of 02 molecules to я'Д8 and states // Proceedings ofSPIE. - 2002. - Vol. 4760. - P. 609-620.

A6. A.M. Старик, Н.С. Титова. О кинетике инициирования детонации в сверхзвуковом потоке смеси Н2+О2 (воздух) при возбуждении молекул 02 резонансным лазерным излучением II Кинетика и катализ. - 2003. - Т. 44, № 1. - С. 35-46.

А7. А.М. Старик, Н.С. Титова. О кинетических механизмах инициирования горения водородно-кислородных смесей при возбуждении электронных степеней свободы молекулярного кислорода лазерным излучением // Журнал технической физики. 2003. - Т. 73, № 3. - С. 59-68.

А8. A.M. Старик, Н.С. Титова. О возможности интенсификации цепных реакций в горючих смесях при возбуждении электронных со-

стояний молекул 02 лазерным излучением II Доклады Академии наук. - 2003. - Т. 391, № 4. - С. 471-477.

А9. A.M. Starik, B.I. Loukhovitski, N.S. Titova. Control of combustion and detonation by means of resonance laser radiation: analysis and potentialities // Proceedings of SPIE. - 2003. - Vol. 4948. - P. 411-417.

A10. A.M. Старик, H.C. Титова. О возможности инициирования горения смесей СН4-02 (воздух) при возбуждении молекул 02 лазерным излучением // Физика горения и взрыва. - 2004. - Т. 40, № 5. - С. 315.

All. A.M. Старик, H.C. Титова. Инициирование горения метано-воздушной смеси в сверхзвуковом потоке за ударной волной при возбуждении молекул 02 лазерным излучением Н Журнал технической физики. - 2004. - Т. 74, вып. 9. - С. 15-22.

А12. A.M. Starik, B.I. Loukhovitski, N.S. Titova. Laser-induced excitation of target molecules as an efficient approach to control the combustion and technological chemical processes // Proceedings of SPIE. - 2005. - Vol. 6053.-P. 1-12.

A13. A.M. Старик, Б.М. Луховицкий, H.C. Титова. Об инициировании горения в сверхзвуковом потоке водородовоздушной смеси излучением СОг-лазера // Механика жидкости и газа. — 2005. - № 2. - С. 157-167.

А14. A.M. Старик, Н.С. Титова. Об интенсификации окисления богатых метановоздушных смесей при возбуждении молекул 02 в состояние a1 Ag // Кинетика и катализ. - 2006. - Т. 47, № 4. - С. 504-513.

А15. A.M. Starik, N.S. Titova, L.V. Bezgin, V.I. Kopchenov, V.V. Naumov. Control of combustion by generation of singlet oxygen molecules in electrical discharge // Czechoslovak Journal of Physics. - 2006. - Vol. 56, Supplement 2. - B1357-B1363.

A16. B.E. Козлов, A.M. Старик, H.C. Титова. Об ускорении горения во-дородно-воздушной смеси при возбуждении молекул 02 в состояние a1 Ag // Физика горения и взрыва. - 2008. - Т. 44, Л» 4. - С. 3-12.

А17. A.M. Starik, N.S. Titova, L.V. Bezgin, V.I. Kopchenov. The promotion of ignition in a supersonic H2-air mixing layer by laser-induced excitation of 02 molecules: Numerical study // Combustion and Fflame. - 2009. -Vol. 156, no. 8. - P. 1641-1652.

A18. A.M. Starik, V.E. Kozlov, N.S. Titova. On mechanisms of a flame velocity increase upon activation of 02 molecules in electrical discharge // Journal of Physics D: Applied Physics. - 2008. - Vol. 41, no. 12. - 125206.

A19. A.M. Starik, N.S. Titova, L.V. Bezgin, V.I. Kopchenov. Initiation of diffusion combustion in a supersonic flow of H2-air mixture by electrical-discharge-excited oxygen molecules // Journal of Physics D: Applied Physics. - 2008. - Vol. 41, no. 12. -125210.

А20. A.M. Старик, П.С. Кулешов, Н.С. Титова. К вопросу об инициировании горения водородно-воздушных смесей лазерным излучением // Журнал технической физики. - 2009. - Т. 79, вып. 3. - С. 28-38.

А21. В.В. Смирнов, A.M. Старик, О.М. Стельмах, Н.С. Титова, Д.Н. Козлов, В.И. Фабелинский. Интенсификация горения водородно-кислородных смесей в дозвуковом потоке при возбуждении молекул 02 в состояние а1 А% электрическим разрядом II Доклады Академии наук. - 2009. - Т. 424, № 6. - С. 760-764.

А22. A.M. Starik, P.S. Kuleshov, N.S. Titova. Comprehensive analysis of combustion initiation in methane-air mixture by resonance laser radiation // Journal of Physics D: Applied Physics. - 2009. - Vol. 42, no. 17. -175503.

A23. A.M. Starik, V.E. Kozlov, N.S. Titova. On the influence of singlet oxygen molecules on the speed of flame propagation in methane-air mixture // Combustion and Flame. - 2010. - Vol. 157, no. 2. - P. 313-327.

A24. A.M. Starik, N.S. Titova. Initiation of combustion and detonation in H2+02 mixtures by excitation of electronic states of oxygen molecules // HighSpeed Deflagration and Detonation. Fundamentals and Control / Ed. by G. Roy, S. Frolov, D. Netzer, A. Borisov. - M.: ELEX-KM Publishers, 2001. -P. 63-78.

A25. A.M. Starik, N.S. Titova. Kinetics of detonation initiation in supersonic flow of combustible mixtures under excitation of vibrational and electronic states of reacting molecules // Nonequilibrium processes and their applications / VI International School-seminar, Minsk, Belarus, 2002. - P. 99-103.

A26. A.M. Starik, N.S. Titova. Kinetics of ignition of gaseous mixtures by resonant laser radiation // Confined detonations and pulse detonation engines / Ed. by G.D. Roy, S.M. Frolov, R. J. Santoro, S.A. Tsyganov. - M.: Torus press, 2003.-P. 87-104.

A27. V.V. Naumov, A.P. Chernukho, A.M. Starik, N.S. Titova. Modeling of plasma-chemical initiation of detonation in a supersonic flow of combustible mixtures // Combustion and Atmospheric Pollution / Ed. by G.D. Roy, S.M. Frolov, A.M. Starik. - M.: Torus press Ltd., 2003. - P. 312-317.

A28. A.M. Starik, N.S. Titova, B.I. Loukhovitski. Novel approach to control the combustion processes by laser induced excitation of target molecules // Combustion and Atmospheric Pollution / Ed. by G.D. Roy, S.M. Frolov, A.M. Starik. - M.: Torus press Ltd., 2003. - P. 46-50.

A29. A.M. Starik, V.E. Kozlov, N.S. Titova. Effect of oxygen molecules activation in electric discharge on flame propagation in H2-air and CRt-air mixtures / Explosion dynamics and hazards / Ed. by S.M. Frolov, F. Zhang, P. Wolan-ski. - M.: Torus press, 2010. - P. 79-104.

Титова Наталия Сергеевна

ИССЛЕДОВАНИЕ ВОСПЛАМЕНЕНИЯ И ГОРЕНИЯ ВОДОРОДА И МЕТАНА В ГАЗОВЫХ ПОТОКАХ ПРИ ВОЗБУЖДЕНИИ ЭЛЕКТРОННЫХ СТЕПЕНЕЙ СВОБОДЫ МОЛЕКУЛ КИСЛОРОДА

01.02.05 - Механика жидкости, газа и плазмы

АВТОРЕФЕРАТ

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

Усл.п.л. -1.5 Заказ №07560 Тираж: ЮОэкз.

Копицентр «ЧЕРТЕЖ.ру» ИНН 7701723201 107023, Москва, ул.Б.Семеновская 11, стр.12 (495) 542-7389 www.chertez.ru

 
Текст научной работы диссертации и автореферата по механике, кандидата физико-математических наук, Титова, Наталия Сергеевна, Москва

московским государственный университет

им. М.В. Ломоносова

61 12-1/536

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

Титова Наталия Сергеевна

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

Специальность 01.02.05 - Механика жидкости, газа и плазмы

ДИССЕРТАЦИЯ на соискание ученой степени кандидата физико-математических наук

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

Старик А. М.

Москва-2012

Содержание

Введение...................................................................................................................4

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

1.1. Методы получени синглетного кислорода..........................................15

1.2. Эксперименты по интенсификации горения синглетным кислородом............................................................................................20

1.3. Экспериментальные исследования элементарных процессов с участием синглетного кислорода........................................................23

1.3.1. Реакция 202(а1Аё)=02(Ь11:ё+)+02(Х31:ё").................................23

1.3.2. Реакция 02(Ь11ё+)+М=02(а1Аё)+М.........................................23

1.3.3. Реакция 02(а1Аё)+М=02(Х32:ё")+М........................................ 24

1.3.4. Реакции Н2+02(а1Аё)=Н+Н02, Н2+02(Ь1^ё+)=Н+Н02(А')... 24

1.3.5. Реакция Н+02(а1Аё)=0Н+0................................................... 25

1.3.6. Реакция 02(а1Аё) + Н02...........................................................25

1.3.7. Реакции с 0(1Б).......................................................................26

1.4. Теоретические исследования элементарных процессов с участием молекул синглетного кислорода..........................................................27

1.5. Модели процессов в реагирующих средах при наличии синглетного кислорода.........................................................................28

Глава 2. Теоретический анализ интенсификации цепных реакций и

процесса воспламенения при возбуждении внутренних степеней

свободы молекул...........................................................................................31

2.1. Теоретическое обоснование возможности интенсификации цепных реакций и процесса воспламенения при возбуждении внутренних степеней свободы молекул.............................................31

2.2. Методика расчета констант скорости реакций с участием электронно-возбужденных молекул...................................................39

Глава 3. Интенсификация процессов воспламенения и горения смеси Н2-

воздух при возбуждении электронных степеней свободы молекул 02... 49

3.1. Кинетическая модель процессов в смеси Н2-02(воздух)...................49

3.2. Кинетическая модель для описания воспламенения и горения смеси Н2-воздух при наличии электронно-возбужденных молекул 02.............................................................................................................62

3.3. Термодинамических свойств электронно-возбужденных компонентов..........................................................................................63

3.4. Тестирование модели.............................................................................68

3.5. Воспламенение смеси при воздействии импульса излучения..........78

3.6. Интенсификация процесса воспламенения при инициировании горения за ударной волной..................................................................93

3.7. Распространение ламинарного фронта пламени..............................102

Глава 4. Влияние возбуждения электронных степеней свободы молекул 02

на характеристики воспламенения и горения смеси СЩ-воздух..........111

4.1. Кинетическая модель процессов в смеси СН4-воздух.....................111

4.2. Кинетическая модель для описания воспламенения и горения смеси СН4-воздух при наличии электронно-возбужденных молекул 02...........................................................................................121

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

4.4. Конверсия метана в синтез-газ...........................................................130

4.5. Инициирование воспламенения за ударной волной......................... 137

4.6. Распространение ламинарного фронта пламени..............................142

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

Литература...........................................................................................................157

Введение

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

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

Однако в последние годы проявляется значительный интерес к использованию других физических методов воздействия на реагирующую смесь с тем чтобы ускорить процесс воспламенения. Так, например, обсуждается возможность интенсификации цепных реакций в молекулярных системах, а, следовательно, и процессов воспламенения и горения, путем возбуждения внутренних степеней свободы реагирующих молекул [1—3]. Интерес к этой проблеме обусловлен тем, что химические реакции с участием возбужденных даже в нижние колебательные или электронные состояния молекул протекают на несколько порядков величины быстрее, чем с участием невозбужденных [4, 5]. При этом возникает большое число взаимосвязанных процессов, в которых колебательное или электронное возбуждение может тем или иным способом передаваться от одного компонента к другому и, таким образом, влиять как на локальные, так и на интегральные характеристики процесса горения. При этом реализуются так называемые энергетически разветвленные цепные реакции. В последнее время на этом направлении исследований был достигнут значительный прогресс. Связано это, в первую

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

В данной работе рассматривается метод интенсификации горения, основанный на возбуждении метастабильных электронных состояний <я'д§ и молекул 02 которые с высокой эффективностью можно получить в специально организованном электрическом разряде [6, 7] или при воздействии резонансного лазерного излучения [8]. Поскольку молекулярный кислород (или воздух) является окислителем большинства используемых в промышленности топлив, то данный метод может оказаться эффективным во многих практических приложениях. Всё это определяет актуальность исследований по кинетике процессов в реагирующих средах при наличии в них возбужденных молекул О2(«1Аё) и О^б1!^), созданию новых детальных реакционных механизмов в случае возбуждения электронных состояний молекул 02 либо электрическим разрядом, либо резонансным лазерным излучением, а также оценке эффективности этого метода интенсификации горения.

Цель и задачи исследования

Целью работы является численное исследование возможности ускорения воспламенения (в том числе в сверхзвуковом потоке) и увеличения скорости ламинарного пламени в смесях Н2-воздух и СН4-воздух при возбуждении молекул 02 в метастабильные электронные состояния аА& и а также

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

1. разработать детальные кинетические модели окисления водорода и метана в воздухе при наличии в рассматриваемых смесях электронно-возбужденных молекул 02(а!Аё) и 02(^1Е§+);

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

3. исследовать влияние возбуждения молекул О2 в состояния и при импульсном воздействии лазерного излучения с длиной волны X]=1.268 мкм и 762 нм, соответственно, на время задержки воспламенения, а также на температуру воспламенения в смесях Н2-воздух и СН4-воздух;

4. исследовать неравновесные процессы за фронтом наклонной ударной волны в смесях Н2-воздух и СН4-воздух при возбуждении электронных состояний молекул 02 резонансным лазерным излучением и провести анализ возможности существенного сокращения длины зоны индукции в сверхзвуковом потоке за фронтом ударной волны;

5. проанализировать возможность использования метода возбуждения молекул 02 в состояния аАо и для ускорения конверсии обогащенных топливом метано-воздушных смесей в синтез-газ (Н2+СО) или для осуществления этого процесса при более низкой температуре;

6. изучить влияние наличия возбужденных молекул Ог(а в смесях Н2-воздух и СН4-воздух на скорость и концентрационные пределы распространения ламинарного пламени;

7. провести сравнение энергетической эффективности метода интенсификации горения, основанного на возбуждении молекул 02 в состояния и Ь1Её+, и метода, основанного на нагреве газа;

Научная новизна

1. Предложена полуэмпирическая методика расчета константы скорости реакции с участием электронно-возбужденных молекул СЬСя'Лг) и

в случае, когда константа скорости аналогичной реакции с участием

•2

молекулы 02, находящейся в основном электронном состоянии X известна.

2. Созданы детальные кинетические модели окисления водородно-воздушных и метано-воздушных смесей при наличии в них электронно-возбужденных молекул 02(а1Аё) и 02(&1Её+). Проведено тестирование разработанных

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

3. Проведено планирование эксперимента по влиянию возбужденных молекул 02 на процесс горения водородно-кислородной смеси, а именно, определены параметры течения смеси Н2-02 в проточном реакторе, при которых можно визуально наблюдать влияние наличия молекул синглетного кислорода 02(а1Аё), генерируемого в электрическом разряде, на длину задержки воспламенения.

4. С использованием разработанных моделей проведено исследование влияние возбуждения молекул 02 в состояния а'Дё и Ьхна время задержки воспламенения в замкнутом реакторе и на длину задержки воспламенения в сверхзвуковом потоке за наклонной ударной волной в смесях Н2-воздух и СН4-воздух. Показано, что наличие в горючей смеси электронно-возбужденных молекул 02(а1Аё) и 02(6!Е§+) приводит к сокращению задержки воспламенения и к уменьшению температуры воспламенения рассматриваемых смесей. Также продемонстрирована возможность интенсификации процесса конверсии метана в синтез-газ (Н2+СО) при возбуждении молекул 02.

5. Изучено влияние присутствия молекул 02(а!Аё) на скорость и концентрационные пределы распространения ламинарного пламени в смесях Н2-воздух и СН4-воздух различного состава и показана возможность существенного (до 20-70%) увеличения скорости ламинарного пламени даже при небольшом содержании молекул 02(а1Аё) в 02.

6. Исследованы механизмы протекания цепных реакций при горении смесей Н2-воздух и СН4-воздух, содержащих в небольшом количестве электронно-возбужденные молекулы 02(а!Аё) и 02(й1Хё+). Показано, что определяющую роль в интенсификации воспламенения и горения в таких смесях играют реакции инициирования и разветвления цепи с участием молекул 02(а1Аё).

7. Установлено, что для рассмотренных в работе горючих смесей метод интенсификации горения, основанный на возбуждении молекул 02 в состояния а Д„ и , намного эффективнее метода, основанного на нагреве смеси, при одинаковой подведенной к газу энергии.

Практическая значимость

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

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

На защиту выносятся следующие основные результаты и положения:

1. детальные кинетические модели для описания процессов, протекающих при инициировании воспламенения и горения смесей Н2-воздух и СН4-воздух путем возбуждения молекул 02 в метастабильные состояния а А„ и Ь1

2. результаты численного исследования влияния возбуждения молекул 02 в состояния а Аё и при импульсном воздействии лазерного излучения на время задержки воспламенения, а также на температуру воспламенения в смесях Н2-воздух и СН4-воздух;

3. обоснование возможности существенного сокращения длины зоны индукции в сверхзвуковом потоке за фронтом ударной волны при инициировании горения смесей Н2-воздух и СН4-воздух путем возбуждения электронных состояний молекул 02 резонансным лазерным излучением;

4. анализ возможности использования метода интенсификации горения, основанного на возбуждении молекул 02 в состояния a]Ag и ¿1Dg+, для ускорения конверсии обогащенных топливом метано-воздушных смесей в синтез-газ (Н2+СО) или для осуществления этого процесса при более низкой температуре;

5. результаты численного исследования влияния наличия молекул 02(<2!Ag) в смесях Н2-воздух и СН4-воздух на скорость и концентрационные пределы распространения ламинарного пламени;

6. анализ механизмов развития цепного процесса при инициировании горения смесей Н2-воздух и СН4-воздух путем возбуждения молекул 02 в состояния a1 Ag и b1Sg+;

7. сравнительный анализ эффективности метода интенсификации горения, основанного на возбуждении молекул 02 в электронные состояния а А„ и

g+, и метода, основанного на нагреве смеси.

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

Основные результаты докладывались на 21 российской и международной

конференциях:

1. International Colloquium on Control of Detonation Processes, 4-7 July 2000, Moscow, Russia;

2. VIII Всероссийский съезд по теоретической и прикладной механике, 23-29 августа 2001, Пермь, Россия;

3. International Colloquium on Advances in Confined Detonations, 2-5 July 2002, Moscow, Russia;

4. VI International School-Seminar "Nonequilibrium Processes and their Applications", 31 August - 5 September 2002, Minsk, Belarus;

5. 25th International Congress on High-Speed Photography and Photonics, 29 September - 4 October, 2002, Beaune, France;

6. International Symposium on Combustion and Atmospheric Pollution, 8-11 July

2003, St. Petersburg, Russia;

7. 7th International Conference on Laser Ablation, October 5-10, 2003, Crete, Greece;

8. International Colloquium on Application of Detonation for Propulsion, 6-9 July

2004, St. Petersburg, Russia;

9. International Workshop "Nonequilibrium processes in combustion and plasma based technologies", 21-26 August 2004, Minsk, Belarus;

10.17th International Symposium on Air Breathing Engines, 4-9 September 2005, Munich, Germany;

11.2nd International Symposium on Nonequilibrium Processes, Combustion and

Atmospheric Phenomena, 3-7 October 2005, Sochi, Russia; 12.ICONO/LAT, 28 May - 1 June 2007, Minsk, Belarus;

rd

13.3 International Symposium on Non-Equilibrium Processes, Plasma, Combustion

and Atmospheric Phenomena, 25-29 June 2007, Sochi, Russia; 14.IEEE International Conference on Plasma Science, 15-19 June 2008, Karlsruhe, Germany;

15.7th International Symposium on Hazards, Prevention, and Mitigation of Industrial

Explosions, 7-11 July 2008, St.-Petersburg, Russia; 16.6th International Colloquium on Pulse and Continuous Detonations, 10-12

November 2008, Moscow, Russia; 17.6th International Seminar on Flame Structure, Brussels, Belgium, 14-17 September 2008;

18.19th International Symposium on Plasma Chemistry, Bochum, Germany, 26-31 July 2009;

19.22nd International Colloquium on the Dynamics of Explosions and Reactive Systems, 27-31 July 2009, Minsk, Belarus;

20.4th International Symposium on Non-Equilibrium Processes, Plasma, Combustion,

and Atmospheric Phenomena, Sochi, Russia, 5-9 October 2009; 21.Ill Международная научно-техническая конференция "АВИАДВИГАТЕЛИ XXI ВЕКА", 30 ноября-3 декабря 2010, ЦИАМ, Москва, Россия.

Разработанные в диссертационной работе кинетические модели окисления смесей Н2-воздух и СН4-воздух при наличии в них молекул синглетного кислорода используются различными исследовательскими группами Украины, Беларуси, Франции, США. Реакции с участием молекул 02(ó1Ag) и соответствующие им константы скоростей включены в кинетическую базу данных [9].

Результаты работы обсуждались на семинаре по механике сплошных сред под руководством академика РАН А.Г. Куликовского, проф. В.П. Карликова и чл.-корр. РАН О.Э. Мельника в НИИ механики МГУ им. М.В. Ломоносова в 2011 г.

Публикации и личный вклад автора

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

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

Диссертация состоит из введения, четырех глав и заключения. Работа изложена на 176 листах, содержит 65 рисунков, 12 таблиц и список библиографических ссылок из 188 источников.

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

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