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

• Анохин, Евгений Михайлович
• кандидат физико-математических наук
• 2005

Оглавление диссертации кандидат физико-математических наук Анохин, Евгений Михайлович

Оглавление

1 ВВЕДЕНИЕ

2 ЭКСПЕРИМЕНТАЛЬНОЕ ОБОРУДОВАНИЕ И МЕТОДИКА ИЗМЕРЕНИЙ

2.1 Ударная труба.

2.2 Экспериментальная установка.

2.3 Измерение скорости ударной волны.

2.4 Определение параметров газа за ударной волной.

2.5 Регистрация собственного излучения газа.

2.5.1 Регистрация собственного излучения с помощью фото-электронного умножителя

2.5.2 Регистрация собственного излучения с помощью стрик-камеры.

2.5.3 Регистрация излучения ICCD камерой LaVision PicoStar HR12.

2.6 Абсолютная калибровка оптических систем.

3 КИНЕТИЧЕСКАЯ СХЕМА

3.1 Состав кинетической схемы.

3.2 Термодинамика веществ.

3.3 Уравнения химических реакций.

3.3.1 Реакции диссоциации.

3.3.2 Реакции ионизации.

3.3.3 Реакции переноса заряда.

3.3.4 Обменные реакции

3.3.5 Реакции электронного возбуждения.

3.3.6 Излучение.

3.3.7 Реакции диссоциативного тушения.

3.4 Колебательный энергообмен при столкновении нейтральных частиц

3.4.1 Процессы колебательного энергообмена в системе C-N-0.

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

3.5.1 Поуровиевая кинетика реакций.

3.5.2 Изменение порога реакции и вероятности перехода при возбуждении реагентов

3.5.3 Оценка величины энергетического барьера реакции.

3.5.4 Оценка вероятности перехода для выделенных уровней.

3.6 Расчет уровневых констант диссоциации.

ОГЛАВЛЕНИЕ

4 ЭКСПЕРИМЕНТАЛЬНЫЕ ДАННЫЕ И СРАВНЕНИЕ С РАСЧЕТОМ

4.1 Переход к лабораторной системе коордииат.

4.2 Смесь CO-N2 70

4.3 Смесь C02-N2-Ar 9.7-0.3-90.

4.4 Смесь C02-N2-Ar 48.5-1.5

4.5 Смесь C02-N2 97

Введение диссертации (часть автореферата) на тему «Физико-химические процессы за фронтом сильных ударных волн в смесях CO2-N2»

В последнее время снова наблюдается рост интереса мировой общественности к развитию космических программ. Американское космическое агентство (NASA) в рамках своей программы по исследованию Марса (Mars Exploration Program) запустило два марсохода 10 июня и 7 июля 2003 г., совершивших успешную посадку на Марс 3 и 24 января 2004 г. соответственно. Европейское космическое агентство (ESA) разрабатывает собственную программу исследований солнечной системы "Аврора"(Aurora program). В целях продолжения изучения Марса в нее входят запуск марсохода на поверхность планеты и запуск возвращаемого модуля, который попытается доставить на Землю образцы марсианской породы. В предварительные испытания входят выведение небольшого летательного аппарата на эллиптическую земную орбиту для эмуляции возвращения космического аппарата с Марса на больших скоростях, а также запуск на Марс небольшого космического корабля, который должен проверить возможность аэродинамического торможения в атмосфере Марса. Пилотируемый полет на Марс в рамках программы "Аврора" запланирован на 2025-2030 г.г. (при условии, что все предыдущие миссии окажутся успешными). История создания космических аппаратов показала, что одной из важных проблем проектирования космического аппарата является разработка и тестирование его тепловой защиты при спуске аппарата на поверхность Земли. Такие исследования широко велись в XX веке при проектировании космических кораблей многоразового использования (Space Shuttle в США, 11 Буран "в СССР) и разработке баллистических ракет. Множество фундаментальных и прикладных работ было посвящено решению данной задачи. Успешно летающие космические корабли многоразового использования позволяют утверждать, что проблема тепловой защиты для космического аппарата, входящего в атмосферу Земли с орбитальной скоростью (6 - 7 км/с), успешно решена. Однако для исследования планет солнечной системы необходимо иметь данные о процессах, протекающих в атмосферах других планет, так как параметры процессов в ударном слое, образующемся перед входящим в верхние слои атмосферы космическим аппаратом, в химических смесях различного состава могут сильно отличаться от тех, которые характерны для земных условий. Более того, при возвращении таких космических аппаратов на Землю их скорость составит уже не б - 7 км/с, как при орбитальных полетах, а 10 - 15 км/с, выводя нас в новую область параметров гиперзвуковой аэродинамики. Подобные обстоятельства ставят перед нами с новой актуальностью задачу создания эффективной тепловой защиты космического корабля в отличных от земной атмосферы условиях и новом диапазоне скоростей. В ударных слоях, возникающих при движении гиперзвуковых космических аппаратов в атмосфере, происходит новый для классической газовой динамики процесс переноса энергии излучением, который может существенно повлиять на нагрев аппарата и на параметры газа в ударном слое. Поглощение и излучение квантов света в газе происходит при изменении энергетического состояния частиц газа. Вид переходов, а следовательно, энергия и частота квантов зависят от природы газа, его температуры и давления. Поэтому знаний, накопленных по переносу тепла излучением в воздушной атмосфере, не хватает для адекватных оценок радиационного потока в атмосферах других планет.

Что касается Марса, то вторая космическая скорость для него составляет 5 км/с и полет к нему от Земли по траектории наименьшей энергии приводит ко входу в атмосферу Марса со скоростью примерно 6 км/с. Такой полет занимает примерно один год в одну сторону, что не очень удобно для пилотируемого полета. Возможен более быстрый полет при использовании более энергоемкой траектории, однако это приведет к более скоростному входу в верхние слои атмосферы. Поиск компромисса между скоростью входа в атмосферу и длительностью путешествия приводит на сегодняшний день к 5-8 месяцам пути в одну сторону и скорости входа в атмосферу 7-9 км/с. Использование аэродинамического торможения позволяет существенно снизить начальную массу корабля для полетов такого типа.

Марсианская атмосфера состоит из 1-1.5% Аг, 2.5-3% N2, оставшаяся часть атмосферы приходится на СО2. Значит основными компонентами в ударном слое будут С, СО, С2 и CN, которые, как известно, обладают сильно излучающими состояниями. Так же известно, что в земной атмосфере при скоростях около 10 км/с большую роль играет излучение из неравновесной зоны. В случае марсианской атмосферы неизвестно соотношение излучения из равновесной и неравновесных зон при различных скоростях полета, поэтому было бы интересно узнать размер неравновесной зоны и ее роль в переносе тепла излучением.

На сегодняшний день ряд отечественных и зарубежных исследовательских коллективов вовлечен в решение данной проблемы. Как уже было сказано, спектральный состав излучения слоя, нагретого ударной волной, напрямую зависит от состава газа, в котором движется летательный аппарат, а также от тех физико-химических процессов, которые протекают в ударном слое, непосредственно перед термозащитной поверхностью аппарата. Поэтому наряду с небольшим количеством экспериментальных работ существует множество расчетов, которые, опираясь на имеющиеся экспериментальные данные, пытаются предсказывать спектральный состав и интенсивность излучения газа за фронтом ударной волны. Из работ, непосредственно посвященных исследованию излучения атмосферы Марса, стоит отметить работу С.А. Лосева и др. [1], в которой рассматривается излучение за фронтом ударной волны в смесях C02-N2-Ar в диапазоне скоростей 3-4 км/с. Эксперименты проводились в двух типах смесей: 1 - C02-N2-Ar = 0.096:0.003:0.901, 2 - C02-N2-Ar = 0.048:0.0015:0.9595 на однодиафрагменной ударной трубе. Были проведены три серии экспериментов: первая - в смеси 1 для начального давления Р1=1 торр и средней скорости ударной волны Vs=3.45 км/с; вторая - в смеси 1 для начального давления Р1=5 торр и средней скорости ударной волны Vs=2.94 км/с; третья - в смеси 2 для начального давления Р1=0.5 торр и средней скорости ударной волны Vs=3.75 км/с. В работе были получены экспериментальные спектры излучения за фронтом ударной волны в диапазоне 350 - 850 нм. Диагностика велась с помощью видиконовой камеры в широком диапазоне длин волн в режиме интегрирования сигнала и проводилось измерение интенсивности излучения на длинах волн 387, 650 и 470 нм с помощью ФЭУ. Перечисленные линии соответствуют фиолетовой системе CN, красной системе CN и полосе Свана молекулы С2 соответственно. Также измерялась интенсивность излучения атомов кислорода (777.19, 777.42 и 777.54 нм). В работе отмечается, что основными излучающими молекулярными полосами являлись фиолетовая и красная системы CN и полоса Свана молекулы С2. На всех экспериментальных спектрах присутствовали атомарные линии кислорода и углерода. Работа [1] наиболее близка к теме данных исследований. Однако вышеописанные эксперименты, в силу конструктивных особенностей установки, проводились в достаточно узком диапазоне скоростей ударной волны. Используемый интегральный тип диагностики не позволяет проследить динамику изменения интенсивности излучения, за исключением нескольких длин волн, для диагностики которых использовались ФЭУ. Также интересен факт наличия или отсутствия излучения четвертой положительной системы СО, излучение которой лежит в области вакуумного ультрафиолета, на этот вопрос работа С.А.Лосева и его коллег также не дает ответа.

Под руководством Г.Г. Черного и С.А. Лосева была выполнена работа [2], в которой один из разделов был посвящен экспериментальному исследованию излучения, однако экспериментальные результаты, приводимые в отчете, в основном те же, что и в работе [1]. Отдельная часть работы посвящена экспериментальным данным, которые получены на электроразрядной ударной трубе в ЦАГИ, к сожалению отсутствует детальное описание экспериментальной установки и используемой диагностики, что не позволяет сделать адекватные выводы о результатах по этому разделу. Стоит, однако, отметить, что в работе обращается внимание на важность излучения четвертой положительной системы СО как системы, которая вносит основной вклад в неравновесное излучение за фронтом ударной волны. Такой вывод был сделан на основе того, что интенсивность излучения начинает расти по направлению к области вакуумного ультрафиолета.

В ряде работ для исследования спектрального состава излучения используются различные типы разрядов. Например, в [3] для сравнения с численным кодом и получения эмиссионного спектра смеси C02-N2=97:3 использовались электродуговой и СВЧ разряд. Из работы следует, что спектральный состав излучения газовой смеси зависит как от типа возбуждения газа, так и от вкладываемой в разряд мощности. Для дугового разряда измерения проводились в диапазоне длин волн 275 - 950 нм, основной вклад в излучение вносили фиолетовая система CN и система Свана молекулы С2, также значительную интенсивность имели атомарные линии кислорода. В случае СВЧ разряда, как уже упоминалось, значительную роль играли начальное давление и мощность разряда. При давлении 10 торр и мощности 800 Вт, основной вклад в излучение вносили третья положительная система и система Ангстрема молекулы СО, при давлении 50 торр и мощности 1300 Вт, основной вклад в излучение вносит фиолетовая система CN.

В работе [4] для исследования спектрального состава излучения также использовался СВЧ разряд. Проведено сравнение спектрального состава атмосферы Земли и Марса в диапазоне 150 - 400 нм, показано, что в отличие от земной атмосферы, значительная часть излучения атмосферы Марса сосредоточена в области вакуумного ультрафиолета, за что ответственна четвертая положительная система СО.

Также проводились эксперименты по исследованию неравновесного излучения в СО2 [5]. Однако несмотря на все эти исследования до сих пор неясна роль радиационного теплопереноса и доля излучения неравновесной зоны в Марсианской атмосфере.

Ударные трубы широко используются для исследования высокотемпературных физико-химических процессов, с их помощью были измерены константы скоростей многих важных в горении реакций. Наряду с "горельны-ми" реакциями также можно измерять и другие не менее важные реакции, например реакции протекающие за фронтом ударной волны при входе космического корабля в планетарную атмосферу. При наличии констант скоростей необходимых химических реакций, не возникло бы необходимости проводить экспериментальные исследования и вся задача свелась бы к решению системы дифференциальных уравнений. Поэтому ряд исследований направлен на получение констант скоростей различных реакций, в том числе с участием молекул CN, СО, СО2, N2 и т.д. Достаточно широкий спектр экспериментальных данных существует по реакциям диссоциации СО2, СО и CN, а также по реакциям обмена, содержащим указанные компоненты. Обзор по таким реакциям можно найти в [6] и [7]. Однако большинство экспериментальных данных охватывают диапазон температур до 15000 К, поступательная же температура за ударной волной в Марсианской атмосфере при скорости в 6 км/с составляет 25000 К, а при скорости 9 км/с - 59000 К, поэтому для описания процессов за фронтом ударной волны приходится пользоваться экстраполяцией более низкотемпературных данных по константам скоростей химических реакций.

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

В случае неравенства поступательной, вращательной, колебательной и электронной температур между собой, но при сохранении внутри подсистем больцмановского распределения, удобно пользоваться двух- или многотемпературными моделями. В таких случаях рассматриваются константы скоростей реакций k(T,Tv), зависящие от различных температур, в данном случае от поступательной и колебательной. Такие константы получаются преобразованием равновесных констант скоростей химических реакций: k(T,Tv) = Z(T,Tv)-k°(T) (1.1) где Z{T,TV)— фактор неравновесности, к°(Т)— равновесная константа скорости реакции.

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

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

Наиболее приближен к действительности микроскопический уровень описания, когда среда рассматривается как множество взаимодействующих частиц, движение которых описывается классическими или квантовыми уравнениями. Примером может служить численное решение задачи о возбуждении колебаний и диссоциации двухатомных молекул при столкновении с атомами инертного газа за фронтом сильной ударной волны [9]. Однако такой подход требует огромных затрат машинного времени, поэтому таким методом обсчитывают единичные реакции. Например в [10] сделан расчет методом квазиклассичсских траекторий обменной реакции CO+N —» CN+0.

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

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

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

Глава 2

ЭКСПЕРИМЕНТАЛЬНОЕ ОБОРУДОВАНИЕ И МЕТОДИКА ИЗМЕРЕНИЙ

2.1 Ударная труба

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

Ударная труба представляет собой длинную трубу, обычно круглого или прямоугольного сечения, которая разделяется тонкой диафрагмой на две части [12] . Одна из них, камера низкого давления (КНД), наполняется исследуемым газом. Во вторую часть, камеру высокого давления (КВД), напускается рабочий газ. Обычная длина ударной трубы составляет несколько метров, а внутренний диаметр - порядка нескольких сантиметров. Давление исследуемого газа бывает различным: от долей торра до нескольких сотен торр. В камере высокого давления необходимо сделать как можно более высокое давление. В нужный момент диафрагма, разделяющая камеры низкого и высокого давлений разрывается и по исследуемому газу в камере низкого давления распространяется ударная волна, которая практически мгновенно нагревает исследуемый газ до высокой температуры. Способы создания высоких температур и давлений в камере высокого давления ударной трубы могут быть различными: от простого нагнетания давления до электрического разряда или взрыва горючей смеси. Ударная волна распространяется в виде скачка уплотнения в исследуемом газе, длина такого скачка определяется интенсивностью ударной волны и свойствами исследуемого газа. За скачком уплотнения следует область перемешивания рабочего и исследуемого газов, так называемая контактная поверхность. Из-за образования пограничного слоя и как следствие сужения эффективного сечения трубы, расстояние между фронтом ударной волны и контактной поверхностью в режиме установившегося течения практически не изменяется. Данное обстоятельство позволяет оценить время диагностики за падающей ударной волной, для сильных ударных волн оно составляет от десятков до сотен микросекунд.

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

Р± Pi

27! М2 7i-l 7i + l 7I +1

2.1) где р - давление, 7 - показатель адиабаты, с - скорость звука в газе. Индексы 1 относятся к исследуемому газу, 4 - к рабочему газу в камере высокого давления.

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

Заключение диссертации по теме «Физика плазмы», Анохин, Евгений Михайлович

ЗАКЛЮЧЕНИЕ

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

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

Проведенные эксперименты по исследованию спектрального состава излучения за фронтом ударных волн в смесях C0-N2(70:30), C02-N2(97:3), С02-N2-Ar(9.7:0.3:90), C02-N2-Ar(48.5:1.5:50) при скоростях ударных волн от 4 до 8 км/с и начальных давлениях смеси 0.5-10 торр показали, что основным источником излучения является радиационный переход CN(B2E+) —» CN(X2E+).

Излучение четвертой положительной системы СО растет относительно излучения фиолетовой системы CN с ростом температуры и при скоростях 8-9 км/с может стать определяющим в процессах радиационного тсплоперсноса.

Как показывают экспериментальные данные, интенсивность излучения системы Свана молекулы Сг является существенным в смесях CO-N2(70:30),

C02-N2-Ar(9.7:0.3:90), C02-N2-Ar(48.5:1.5:50) в исследованном диапазоне скоростей ударных волн.

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

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

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

1. С.А. Лосев, П.В. Козлов, J1.A. Кузнецова, В.Н. Макаров, Ю.В. Романен-ко, С.Т. Суржиков, Г.Н. Залогин. Излучение смеси C02-N2-Ar в ударных волнах: эксперимент и теория. Препринт 40-98, Институт Механики МГУ. Москва 1998 г.

2. Проблемы аэротермобаллистики, радиационной газодинамики и тепломассообмена в экспедициях возвращения с планетными образцами. Проект МНТЦ №1549, под рук. Г.Г. Черного и С.А. Лосева, Москва 2000-2003 г.г.

3. М. Lino da Silva. Martian Plasma: Arc-Jet Experiments and the SESAM Simulation Tool. Euromech colloquia on Aerodynamics and Thermochemistry of High Speed Flows, Marseille, France, 2002.

4. E.M. Anokhin, S.M.Starikovskaia, A.Yu.Starikovskii. Energy transfer in hypersonic plasma flow and flow structure control by low temperature nonequilibrium plasma. 42nd AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, USA, 2004.

5. Nealy J.E. An Experimental Study of Ultraviolet Radiation Behind Incident Normal Shock Waves in C02 at Venusian Entry Speeds. AIAA Paper 75-1150, Sept. 1975

6. Baulch D.L., Drysdale D.D., Duxbury J. and Grant S.J. Evaluated Kinetic Data for High Temperature Reactions. Vol. 3, Homogeneous Gas Phase

7. Reactions of the 02-03 System, the C0-02-H2 System, and of Sulphur-Containing Species. Butterworth, London, 1987.

8. Физико-химические процессы в газовой динамике. Под ред. Г.Г. Черного и С.А. Лосева. T.l М: Издательство Московского Университета, 1995.

9. Б.В. Куксенко, С.А. Лосев.//Докл. АН СССР. 1969. Т.185. М. С.69

10. М.Ю. Погосбекян, С.А. Лосев. Исследование реакции CO+N —> CN+0 методом квазиклассичсских траекторий с использованием вычислительного комплекса "MD Trajectory ".//Химическая физика. 2003. Т. 22. С. 38-46.

11. Bhaskaran R.F., Roth P. The shock tube as wave reactor for kinetic studies and material systems.// Progress in energy and combustion science. 2002.

12. Я.Б.Зельдович, Ю.П.Райзер. Физика ударных волн и высокотемпературных гидродинамических явлений. Издание второе, дополненное. М.: Наука, 1966.

13. Е.В. Ступоченко, С.А. Лосев, А.И. Осипов. Релаксационные процессы в ударных волнах. М.: Наука, 1965.

14. J. Kiefer. Laser-schlieren technique in shock-tube kinetics.// Shock waves in chemistry. Ed. by A. Lifshits, NY, 1981. P. 219.

15. Физические величины: Справочник. Под ред. И.С. Григорьева, Е.З. Мсй-лихова. М.:Энергоатомиздат, 1991.

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

17. Термодинамические свойства индивидуальных веществ. Справочное издание в четырех томах. Издание третье, переработанное и расширенное. Под редакцией В.П. Глушко. М.:Наука, 1978-1982 гг.

18. JI.A. Кузнецова, Н.Е. Кузьменко, Ю.Я. Кузяков, Ю.А. Пластинин. Вероятности оптических переходов двухатомных молекул. Под ред. Р.В. Хох-лова. М.: Наука, 1980

19. Laher R.R., Gilmore F.R. Improved fits for the vibrational and rotational constants of many states in nitrogen and oxygen. //J. Phys. and Chemical Reference Data, 20, 1991

20. Baulch D.L., Drysdale D.D., and Home D.G. Evaluated Kinetic Data for High Temperature Reactions. Vol. 2, Homogeneous Gas Phase Reactions of the H2-N2-O2 System. Buttcrworth, London, 1973.

21. Физико-химическая кинетика в газовой динамикс./Под ред. С.А. Лосева, О.П. Шаталова. М.: Изд-во Моск. Ун-та, 1986 г.

22. Park С., Howe J.T., Jaffe R.L. and Candler G.V. Review of Chemical-Kinetic Problems of Future NASA Missions, II: Mars Entries.//Journal of Thermophysics and Heat Transfer. V. 8, N. 1, 1994.

23. Park C., JafTe R.L. and Partridge H. Chemical-Kinetic Parameters of Hyperbolic Earth Entry.// Journal of Thermophysics and Heat Transfer. V. 15, N. 1, 2001.

24. Н.Е. Афонина, В.Г. Громов. Исследование на основе модели вызкого ударного слоя течения в области торможения при входе космического аппарата в марсианскую атмосферу.//Препринт №31-97 Ин-та Механики МГУ, М.1997.

25. В.Н. Макаров. Определение механизма физико-химических процессов в высокотемпературном воздухе.// Прикладная механика и техническая физика. Т. 37, № 2, 1996.

26. Л.Б. Ибрагимова. Константы скорости химических реакций в высокотемпературном газе С02.//Математическое моделирование. Т. 12, № 4, 2000.

27. Gupta R.N., Lee К.Р. An aerothcrmal study of MESUR Pathfinder aeroshell. AIAA Paper 94-2025,1994.

28. Неравновесная колебательная кинетика. / Под. ред. М.Капителли. М.Мир. 1989.

29. Никитин Е.Е., Осипов А.И., Уманский С.Я. // Химия плазмы. М. Энер-гоатомиздат. 1989. Вып. 15.

30. Никитин Е.Е., Уманский С.Я., Шалашилин Д.В. // Химическая физика. 1989. Т. 8. N 8. С. 1011.

31. Moor С.В. // J.Chem.Phys. 1965. V. 43. N 9. Р. 2979.

32. Sharma R.D., Brau С.А. // J.Chem.Phys. 1960. V. 50. N 2. P. 924.

33. Лозовский В.А., Саркисов О.М., Уманский С.Я., Цанава Р.А. // Химическая физика. 1983. Т. 2. N 2. С. 201.

34. Никитин Е.Е. Теория элементарных атомно-молекулярных процессов в газах. М. Химия. 1970.

35. Корниенко И., Уманский С.Я. //Химическая физика. 1982. Т. 1. N 7. С. 901.

36. В.А.Салышков, А.М.Старик // ТВТ. 1995. Т.ЗЗ. N 1. С.121.

37. О.В.Ачасов, Д.С.Рагозин // Константы колебательного энергообмена в лазерно-активных средах с добавками О2, Н2, Н2О, СО. Препринт N 16. Минск: ИТМО АН БССР. 1986. 53 с.

38. А.П.Зуев, Б.К.Ткаченко // Хим.физика. 1988. Т.7. N 11. С.1451.

39. А.П.Зуев, С.А.Лосев, А.И.Осипов, А.М.Старик // Хим.физика. 1992. Т.П. N 1. С.4.

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

41. Диссертационная работа на соискание степени доктора физико-математических наук "Нелинейные волны и энергообмен в химически активных системах", Долгопрудный, 2000.

42. Landau L. // Phys. Z. Sow.,1942. V.2. P.46.

43. Zener С. // Proc. Roy. Soc., 1932. V.A137. P.696.

44. JI.A. Кузнецова, С.Т. Суржиков. Атлас спектральных сечений поглощения электронных и колебательных систем полос двухатомных молекул. Препринт № 603. Институт проблем механики РАН. Москва, 1997.