О влиянии моделей обменных химических реакций на параметры течения воздуха за сильными ударными волнами Текст научной статьи по специальности «Физика»

УДК 533.6.011

Вестник СПбГУ. Сер. 1. Т. 1 (59). 2014. Вып. 1

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

О. В. Кунова, Е. А. Нагнибеда

Санкт-Петербургский государственный университет,

Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7/9

В работе численно исследованы неравновесные течения пятикомпонентной воздушной смеси N2/02 /МО/К/О за сильными ударными волнами на основе наиболее строгого поуровнево-го приближения. Уравнения для заселенностей колебательных уровней молекул и числовых плотностей атомов совместно с уравнениями сохранения численно решались для одномерного квазистационарного течения. Учитывались различные переходы колебательной энергии, диссоциация, рекомбинация и обменные реакции. На основе полученного решения показано влияние разных реакций на изменение колебательных распределений и макроскопических параметров в релаксационной зоне за фронтом ударной волны. В вычислениях использовались разные модели для зависящих от колебательных уровней коэффициентов скоростей обменных реакций. Изучено влияние выбора модели обменных реакций на заселенности колебательных уровней, состав смеси и температуру газа. Библиогр. 14 назв. Ил. 6.

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

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

В рамках поуровневого приближения рассматривается течение пятикомпонент-ной воздушной смеси (N2, О2, N0, N О) за ударной волной с учетом возбуждения колебательных уровней молекул; принимаются во внимание следующие химические реакции диссоциации, рекомбинации и обмена:

(М = N2, О2, N0, N,0; г, г' —колебательный уровень молекулы).

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

* Работа выполнена при финансовой поддержке РФФИ (гранты №11-01-00408, 12-08-00826) и СПбГУ (НИР 6.38.73.2012).

N^2) + 0 ^ N0(2')

02(г) + N ^ N0(2') + О, N (г) + М ^ N + N + М, 02(г) + М ^ 0 + 0 +М, N0(2) + М ^ N + 0 + М,

(1) (2)

(3)

(4)

(5)

вычислительных ресурсов. В настоящее время детальное поуровневое описание неравновесных процессов вызывает большой интерес. Однако применению данного подхода при изучении различных течений многокомпонентных смесей посвящено небольшое число работ (см., например, [1, 2, 3]). В настоящей работе изучены особенности колебательных распределений и макропараметров воздуха в релаксационной зоне за фронтом ударной волны при использовании разных моделей для определения коэффициентов скоростей реакций (1)—(2), зависящих от колебательных уровней молекул.

Основные уравнения. Система уравнений, описывающих течение воздушной смеси газов в поуровневом приближении, состоит из уравнений кинетики для засе-ленностей колебательных уровней молекул (п^2«о2», «N0») и числовых плотностей атомарных компонентов (п^, «о), а также уравнений сохранения массы, импульса и энергии смеси [4].

В работе при изучении течения за фронтом ударной волны поток считается одномерным и квазистационарным, эффекты вязкости и теплопроводности не учитываются. Используется часто принимаемое предположение (см., например, [2, 3]), что в результате реакций (1)—(2) молекулы N0 образуются в основном колебательном состоянии и остаются в нем в процессе релаксации. В этом случае система уравнений принимает вид

<1{уп_ тЬг

^ - "Г" "N2® Т ^^

Й^Г + + * = 0 1 • • • ^N2 - 1, (6)

(1{уПо2г) „тЬг , т?2^2 , • _ п 1 т (7\

- К02г + К02г + К02г > « — 0, 1, . . . Ьо2 - 1, (7)

3.x

Ь N2 — 1 ^02 — 1

- - - + ^N0 >

3х

»=0 »=0 ЬN2 — 1 Ь^ — 1 Ьо2 — 1

^ - N0 - + > I9)

»=0 »=0 »=0 Ьо2 — 1 Ьо2 — 1 ЬN2 — 1

3х

-2 £ до:? - дм03 - Е Д07»2 + Е (10)

»=0 »=0 »=0

Р0^0 = ри, (11)

+ Р0 = ри2 + р, (12)

22

+ ^ = (13)

Здесь индекс «0» обозначает параметры газа в набегающем потоке, LN2 и Ьо2 —число колебательных уровней молекул, х — расстояние от фронта ударной волны, и — скорость, р — плотность, р — давление, Н — энтальпия в расчете на единицу массы:

н = Е НСУС'

УС = рС/р — массовая доля компонента смеси сорта с,

5 £

кс = —КСТ -\--для с = О,

2 тС

1 V- ,

hc = -RCT + - > £lnci + — для с = N2, 02, NO.

2 Pc^ mc

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

Запишем выражения для релаксационных членов уравнений (6)-(10). Для реакции обмена (Д2^2), диссоциации и рекомбинации (Д2^3) имеем

с>2^2 > N , о , о ^

ДN2i — nNOnNkNO, N2» - nN2¿nOkN2i,NO,

с>2^2 > O^ , N , O

^2» — nNOnOkNO,O2i - "«2inNkо2i,NO,

— "М (^^^^^^гсс,^» - "N2^33,N2i), М

— "М(nOnOkMCJO2i - nO2ikMMss,O2i), М

ДNO гес, NO - nNOkMis fi.NO ),

М

(М — 02, N0, N 0).

Колебательный энергообмен (Д"Ьг — + ) включает УУ- и УУ'-обмены

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

ДУУ \ ^ \ ^ ( 7 -^ ^^ 7 ^^-^т^Л^ \

N2i — ("N2^"N2™'— nN2inN2mkN2^ ) ,

т=0 i' =^1

(т'=т^ 1)

¿02 — 1

ДУУ ^ ^ ( 7 т^П^-^ ^^ 7 ^^-\

O2i — "^т' ^Ю^' —^ — nO2i"^т^^—^ ) ,

т=0 i' ^±1

(т' =т^1)

¿02 — 1

пУУ' /* ^г.т'— т ,O2,m—т'\

ДN2i — (^N2^nO2m'kN2,i'—i - nN2inO2т%2^—^ ^ ,

т=0 ^ =i±1

(т'=т^ 1)

—1

пУУ' /* > N2 ,т'—т , N2,™—тЛ

ДO2i — "N2™'kO2,i'—i - nO2inN2mkо2,i—i' ^ ,

т=0 ^ =i±1

(т'=т^ 1)

TjVT V^ / iM , ,M /7 M >M \\

Rci = nM lnc,i-1kc,i-1^i + nc,i+1kc,i+1^i - nc,i + fcc,i^i+JJ ,

M

c =N2, O2, M = N2, O2, NO, N, O, i = 0,1 ...Lc - 1.

ЗДесь ^NO' kN;°NO> kM2i,diSS. kM2 i,diSS. kMO,diSS. и ^¿^¿±1 -коэффИЦИ-

енты скоростей обменных реакций, реакций диссоциации и энергообмена с размерностью м3/сек; ^Mc N2i, ^Mc o2i, no —коэффициенты скоростей рекомбинации с

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

Для описания коэффициентов скоростей энергообмена использовались формулы теории Шварца, Славского, Герцфельда, приведенные в [5]. Коэффициенты скорости диссоциации определялись с помощью модели Тринора—Маррона [6], обобщенной в [4] для поуровнего приближения. Коэффициенты скоростей обменных реакций рассмотрены ниже.

Модели химической кинетики. Поуровневые модели термически неравновесных химических реакций предложены в ряде работ на основе аппроксимации экспериментальных данных, результатов траекторных расчетов, обобщения модели Три-нора—Маррона для обменных реакций и других соображений. Однако вопрос о выборе адекватной модели для решения конкретных газодинамических задач изучен недостаточно. В настоящей работе изучается влияние моделей обменных реакций на параметры потока воздуха за ударной волной. Для определения коэффициентов скоростей обменных химических реакций (1)-(2) в работе рассматриваются три модели, приведенные в [7, 8, 9].

1. а-модель [7] учитывает эффективность вклада колебательной энергии в энергию активации:

СТо^о (Т ) = А(Т )ехр

Е„ — а£

N2 (о2)

кТ

-в(Еа-ае^)

(14)

где а — коэффициент эффективности колебательной энергии, © — функция Хэвисай-да, Еа — энергия активации реакции. Предэкспоненциальный множитель А(Т) имеет значение, близкое к частоте столкновений при единичной концентрации Z = Д0 — газокинетический радиус столкновений, р — приведенная масса сталкивающихся частиц, а = 0.51 для реакции (1), а = 0.24 для реакции (2).

2. В описанной Полаком и др. модели (см. [8]) предлагается формула, основанная на аппроксимации результатов траекторных расчетов:

(Т) = А %2(о2)*,№(Т ) = А еХР

7—I N

Еа - 1£г

N2 (о2 )

вкТ

©(Еа - 7£^(о2))

(15)

Здесь А « пД2(8кТ/пр)1/2, параметры в и 7 характеризуют участие поступательной и колебательной энергий в реакции: в = 0.90, 7 = 0.52 для реакции (1), в = 0.46, 7 = 0.12 для реакции (2).

3. Варнатцем [9] предложено выражение для коэффициента скорости обменных реакций, основанное на теории химически активных столкновений:

кок^о (Т ) = С (г +1)Тв ехр

где

N2 (о2 )»,№

С = 4.17 • 1012' С = 1.15 - 109'

Еа £

N2 (о2)

кТ

-©(Еа - £

N2^2)

(16)

в = 0 для реакции ^(г) + О ^ N0 + N в =1.0 для реакции 02(г) + N ^ N0 + О.

Все представленные выше коэффициенты скоростей реакций имеют размерность см3/(моль • сек). Некоторые другие уровневые модели приведены в [10, 11], однако их

)

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

Сравнение коэффициентов скоростей реакций, рассчитанных с помощью представленных выше моделей, показано на рис. 1, 2. На рис. 1 показано сравнение тем-

1 1 , О^ ,N,0 ,

пературной зависимости коэффициентов А^'г N0 и «О г N0 при фиксированном колебательном уровне реагента. Видно, что соответствующие коэффициенты а-модели и модели [8] имеют близкие значения и совпадают на верхних уровнях. В то же время коэффициенты, рассчитанные по модели Варнатца, отличаются от найденных с учетом [7], [8] до двух порядков величины, однако с ростом номера уровня эта разница уменьшается. Видно, что во всех трех случаях с увеличением уровня г растут и коэффициенты скоростей реакций.

Далее на рис. 2 приведено сравнение осредненных по больцмановскому распределению коэффициентов скоростей прямых и обратных реакций в зависимости от температуры. Заметим, что в работе [12] показано удовлетворительное согласие коэффициентов прямых реакций (1), (2), найденных на основе осредненной модели Варнатца и по формуле Аррениуса [10]. Отметим также, что в термически равновесном газе обратная реакция (1) с молекулами азота протекает быстрее, чем прямая, а вероятность прямой реакции (2) с молекулами кислорода гораздо выше обратной. Расчеты также показали, что обменные реакции протекают более активно, чем реакции диссоциации и рекомбинации (3)—(5).

Результаты численных расчетов. В ходе численного решения задача об ударной волне разбивалась на два этапа. На первом этапе определялись газодинамические параметры непосредственно за фронтом ударной волны. При этом обмен колебательной энергией и химические реакции внутри фронта считались замороженными с неизменными составом смеси и заселенностями колебательных уровней молекул. На следующем этапе полученные значения использовались в качестве начальных данных при решении системы уравнений (6)—(13), состоящей из ¿N2 + ¿о2 + 3 уравнений ки-

«-(1) —в—(2) -«-(3)

Рис. 1. Коэффициенты скорости прямых обменных реакций в зависимости от температуры. Сплошные линии — г = 0, пунктирные — г = 20, штриховые — г = 30. Маркеры соответствуют трем рассматриваемым моделям. Здесь и далее на графиках: (1)— а-модель, (2) — модель [8], (3) — модель Варнатца.

Рис. 2. Осредненные по больцмановскому распределению коэффициенты скорости прямых (сплошные линии с маркерами) и обратных (штриховые линии с маркерами) обменных реакций в зависимости от температуры.

Рис. 3. Относительная числовая плотность компонент смеси в зависимости от расстояния от фронта ударной волны. а) — молекул N2, О2 и N0; Ь) —атомов N и О.

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

На рис. 3-5 представлены результаты расчетов при следующих условиях в невозмущенном потоке: Т0 = 293 К, р0 = 100 Па, М0 = 15. Предполагалось, что в набегающем потоке смесь состоит только из молекул азота п^ = 0.79«0 и кислорода «о2 = 0.21«0 («0 —числовая плотность смеси в набегающем потоке), а распределения считались равновесными с температурой газа Т> Спектр колебательной энергии моделировался осциллятором Морзе, при этом LN2 =48 и Lо2 = 37.

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

На рис. 3, 4 показано изменение состава газа и его температуры в зависимости от расстояния от фронта ударной волны. При удалении от фронта волны уменьша-

Рис. 4- Температура Т в зависимости от расстояния х.

ется температура газа и числовые плотности молекул азота и кислорода, и, соответственно, увеличиваются числовые плотности атомов. Можно отметить немонотонное изменение числовой плотности молекул окиси азота. Этот факт объясняется тем, что на ранней стадии релаксации быстрая диссоциация молекул кислорода (4) и более медленная диссоциация азота (3) приводят к появлению свободных атомов и увеличению молекул N0 за счет прямых реакций (1)—(2). Дальнейшее накопление атомов азота вызывает активную обратную реакцию (1), приводящую к уменьшению числовой плотности окиси азота. Немонотонный характер изменения числовой плотности молекул N0 за фронтом ударной волны отмечался и при рассмотрении химической релаксации в термически равновесном газе (см., например, [13]). Сравнение макропараметров и колебательных распределений в пятикомпонентной воздушной смеси за фронтом ударной волны, найденных в поуровневом и однотемпературном (термически равновесном) приближениях, показано в [14].

На рис. 5 представлено изменение заселенностей фиксированных колебательных уровней молекул N2 и О2 в зависимости от расстояния от фронта ударной волны. Наблюдается немонотонность изменения заселенностей колебательных уровней и по2г при возрастании х: уменьшение и последующее возрастание заселенности 0-го уровня молекул азота (заселенность нулевого уровня молекул кислорода постоянно уменьшается) и, наоборот, увеличение, а затем уменьшение заселенностей уровней г > 0 молекул N2 и О2. Это объясняется переходом поступательной энергии ударно нагретого газа в колебательную, что приводит к резкому возбуждению колебательных уровней молекул N2 и О2. Затем при уменьшении температуры газа происходит уменьшение числовой плотности возбужденных молекул вследствие дезактивации и химических реакций.

На рис. 6 показаны относительные заселенности колебательных уровней молекул азота и кислорода при различных значениях х. Отметим разный характер распределений п^2г и по2г. Связано это с более высокой скоростью колебательного энергообмена молекул О2, в результате чего стадия возбуждения молекул кислорода за фронтом волны протекает в более узкой зоне, нежели молекул азота.

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

х, ш X, т

(1) ---(2) - — (3)

Рис. 5. Относительные заселенности отдельных уровней молекул в зависимости от расстояния х: а— N2; Ь — О2. Сплошные линии — модель (1), штриховые — модель (2), штрих-пунктирные — модель (3).

а) Ь)

1 - х=10 4 т; 2 - х=10 3 т; 3 - х=0.01 т; 4 - х=0.02 т

Рис. 6. Относительные заселенности молекул в зависимости от уровня г для различных значений х: а — N2; Ь — О2. Сплошные линии — модель (1), штриховые — модель (2), штрих-пунктирные — модель (3).

[8] дают очень близкие значения макропараметров течения (рис. 3, 4) и заселенностей колебательных уровней молекул (рис. 5, 6). Связано это с тем, что коэффициенты скоростей реакций со средних и высших уровней практически не отличаются (верхние кривые на рис. 1, а, Ь), а невозбужденные молекулы гораздо менее активно участвуют в реакциях (1)-(2).

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

плотности молекул кислорода (25% и 22% соответственно). Наименьшие различия — в распределении температур (1.5%), числовых плотностей молекул и атомов азота (1% и 8% соответственно) и атомов кислорода (3%).

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

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

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

Литература

1. Capitelli M., Colonna G., Kustova E., Nagnibeda E. State-to-state kinetics and transport properties in supersonic air nozzle flows // Proceedings of the 4th European Symposium on Aerothermodynamics for Space Vehicles. ESTEC, Noordwijk, The Netherlands, ESA Publication Division, ESA SP-487. P. 137143.

2. Bazilevich S. S., Sinitsyn K. A., Nagnibeda E. A. Non-Equilibrium Flows of Reacting Air Components in Nozzles // AIP Conference Proceedings, 2009. Vol. 1084, N1. P. 843-848.

3. Capitelli M., Armenise I., Gorse C. State-to-state approach in the kinetics of air components under re-entry conditions //J. Thermophys. Heat Transfer. 1997. Vol. 11. N4. P. 570-578.

4. Нагнибеда Е. А., Кустова Е. В. Кинетическая теория процессов переноса и релаксации в потоках неравновесных реагирующих газов. СПб.: Изд-во С.-Петерб. ун-та, 2003. 270 c.

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

6. Marrone P., Treanor C. Chemical relaxation with preferential dissociation from excited vibrational levels // Phys. Fluids. 1963. Vol. 6. N9. P. 1215-1221.

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

8. Полак Л. С., Гольденберг М.Я., Левицкий А. А. Вычислительные методы в химической кинетике. М.: Наука, 1984.

9. Warnatz J., Riedel U., Schmidt R. // Advanced in Hypersonic. Vol. 2: Modeling Hypersonic Flows. Birhauser, Boston, 1992. P. 67-103.

10. Физико-химические процессы в газовой динамике / под ред. Г. Г. Черного и С. А. Лосева. М.: Науч. мир, 2007. Т. 1. 399 с.

11. Ковач Э. А., Лосев С. А., Сергиевская А. Л., Храпак Н. А. Каталог моделей физико-химических процессов. 4. Химические реакции в уровневом представлении // Физико-химическая кинетика в газовой динамике. 2010. №2. С. 413-443.

12. Кунова О. В., Нагнибеда Е. А. Поуровневое описание колебательной и химической релаксации воздухе // Вестн. С.-Петерб. ун-та. Сер. 1: математика, механика, астрономия. 2013. №3. С. 103112.

13. Неравновесные физико-химические процессы в аэродинамике / под. ред. Г. И. Майкапара. М.: Машиностроение, 1972. 344 с.

14. Кунова О. В., Нагнибеда Е. А., Шарафутдинов И. З. Неравновесная колебательная и химическая кинетика в потоках воздуха за сильными ударными волнами // Всероссийская конференция с участием иностранных ученых «Современные проблемы динамики разреженных газов». Новосибирск, 2013. С. 127-129.

Статья поступила в редакцию 24 октября 2013 г.

Сведения об авторах

Кунова Ольга Владимировна — аспирант; kunova.olga@gmail.com

Нагнибеда Екатерина Алексеевна — доктор физико-математических наук, профессор; e_nagnibeda@mail.ru

THE INFLUENCE OF EXCHANGE CHEMICAL REACTIONS MODELS ON THE PARAMETERS OF AIR FLOWS BEHIND STRONG SHOCK WAVES

Olga V. Kunova, Ekaterina A. Nagnibeda

St.Petersburg State University, Universitetskaya nab., 7/9, St. Petersburg, 199034, Russian Federation; kunova.olga@gmail.com, e_nagnibeda@mail.ru

In the paper non-equilibrium flows of the 5-component air mixture N2/O2/NO/N/O behind strong shock waves are numerically studied using the most accurate state-to-state approach. This approach describes strongly non-equilibrium vibrational and chemical kinetics in a flow under the assumption of equilibrium distributions over translation and rotation degrees of freedom. The characteristic times of vibrational and chemical relaxation are supposed to be of the same order as macroscopic parameters changing. The equations for vibrational level populations of molecular species and atomic number densities are coupled to the conservation equations for impulse and total energy and solved numerically for a 1-D quasi-stationary flow. Various vibrational energy transitions as well as dissociation, recombination and exchange reactions are taken into account. Four different models for state-depending rate coefficients of exchange reactions for NO formation are used in numerical calculations and the comparison of the results obtained using various models is presented. The influence of a choice of a model for exchange reactions on level populations, mixture composition and gas temperature in the shock heated air mixture flow behind a shock wave is studied. Refs 14. Figs 6.

Keywords: non-equilibrium kinetics in the air, vibrational distributions, chemical composition, dissociation, recombination, exchange reactions.