Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Проблемы кинетики элементарных химических реакций к семидесятилетию акад. В. Н. Кондратьева [сборник]

.pdf
Скачиваний:
2
Добавлен:
24.10.2023
Размер:
9.51 Mб
Скачать

198

Н. И. Ющевкова

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

Скорости перечисленных реакций могут отличаться на поряд­ ки величин в зависимости от условий истечения, главным образом от степени ионизации газа в области торможения. Для полно­ стью ионизированной плазмы, содержащей двухзарядные и одно­ зарядные ионы, процесс рекомбинации приводит к образованию метастабильных ионов и атомов, что существенно влияет на кине­ тическую картину течения, особенно, при введении молекуляр­ ного газа (SFß, СС12, С02) в струю плазмы [8, 21].

Несмотря на то, что время релаксации возбуждения тв элект­ ронного состояния обычно меньше характерного газодинамиче­ ского времени течения, распределение концентраций возбужден­ ных частиц по уровням существенно отличается от больцмановского вследствие высоких неравновесных концентраций двухзаряд­ ных, однозарядных ионов и электронов в сверхзвуковой части струи плазмы. При этом крайне существенным фактором является релаксация электронной температуры: чем резче падает элект­ ронная температура в струе (при аналогичных параметрах тормо­ жения), тем выше параметр неравновесности х = пе/пе . Однако

большая величина отношения тг/т* [4] обеспечивает квазистацио-

нарное распределение заселенностей, что позволяет отделить систему релаксационных уравнений, описывающих поле Те, пе от системы уравнений для заселенностей. Указанный подход был применен для оценки заселенности электронных состояний азота при течении плазмы азота в сопле [4] и в настоящей работе для расчета течения ионизированного аргона. Из кинетической картины течения в струях газа и плазмы следует, что наибольший практический интерес представляют струи молекулярных газов с колебательной неравновесностью в области I жс химической неравновесностью в остальных областях и неравновесностью элект­ ронных процессов в струе разреженной плазмы. Ниже приводятся результаты исследований указанных течений.

ХИМИЧЕСКАЯ И КОЛЕБАТЕЛЬНАЯ РЕЛАКСАЦИЯ

ВСТРУЯХ ГАЗОВЫХ СМЕСЕЙ

Всверхзвуковых расширяющихся потоках двуокиси углерода

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

Релаксационные процессы в струях низкотемпературной плазмы 199

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

импульса) совместно с

уравнением состояния может быть сведена

к дифференциальному

уравнению (8), описывающему изменение

поступательной температуры вдоль расширяющейся трубки

тока [3]:

dT

2 RT - Щ (Т)

R T ^

И

dl

V2

dl

x ß cA 1 —

/{TZ h

V* - * 2 ^

dS (

 

 

S(l) dl J

X

(8)

-1

S(l)

— расширение

Aef}I

CK0P0CTb течения;

где

у = I /

2 |я„

2$іні (T) +

2

Hi(T) — энтальпия

i-й компоненты; ß*

= p;/p — массовая кон­

центрация

г-й компоненты, отнесенная к молекулярному весу;

 

 

 

трубки тока;

сг — теплоемкость г-й компонен­

ты; R — газовая постоянная; Ае$г) — отклонение от равновесного значения энергий /-й колебательной степени свободы г-й моле­ кулы.

Распределение состава смеси описывается системой уравнений неразрывности, сформулированных с учетом кинетики химических реакций:

dßi

 

 

(9)

Tu

 

 

 

 

 

где lg] — скорости соответствующих

реакций;

 

п = ?iKp^Kp2ßi^‘S'^ ( 2 ß i)

р~~ число

частиц в

единице объема;

икр, F KP — полное число

частиц в единице объема и скорость в

критическом сечении соответственно.

 

 

Рассматривая совокупность колебательных

степеней свободы

многокомпонентной смеси (С02, СО, 0 2, N2) как ансамбль гармони­ ческих осцилляторов, можно получить 18] систему релаксацион­ ных уравнений для средних энергий колебательных степеней

200 Н. И. Юшенкова

свободы Ej в потоке многокомпонентной смеси:

леЧ»)

,

 

 

ч

E f

ѵ ~ і г = ^ 2

2і г Л Ж ' ' Ч {

И

1) + в М 1)] п

ß

ц ? > X

x-4-p* г

 

 

 

 

Je=Si

3

*

ГП(Д

 

 

 

.ку) (т1.. .т^)}

хпm=mi

 

 

 

А£{(/) (*і..

ÄT

П

( 10)]-

 

т=т,

3 т

 

 

 

 

столкновений молекулы і-го сорта с молекулами

где Zir — числоb S H “p[-

 

 

 

 

 

 

г-го сорта в единицу времени; ( Q ^ ^ ' ^ ^ ^ ' " mx))ir— вероятность дезактивации первого уровня ;-й колебательной степени свободы і-й молекулы при одновременном возбуждении К осцилляторов и

дезактивации ц осцилляторов, q — вырождение

осциллятора,

Vjl)32— характеристическая частота ;'-го

колебания і-й молекулы;

ДЕ — дефект энергии,

зависящий от механизма перераспределе­

ния энергий:

 

"‘и.

 

 

 

 

А£ {(/') (*) ("»)}= h

 

 

 

 

2

Ѵт}-

 

 

 

к—kt

7П=ТПі

 

 

 

 

При расчетах учитывалась кинетика следующих процессов:

химических

 

 

 

 

 

 

1. N + CO2= CO + N O

4 . 0 + С 0 2 = С 0 +

0 2

7 .

0 +

N 0 = 0 2 +

N

2 . N + 0 a = N 0 + 0

5 . О “ t~ С 0 + М = С 0 2 - | - М

8 .

0 +

N 2 = N 0 +

N

3. N + N + M = N Ü+ M 6. о + о + м = о 2 + н

9 . 2 N + M = N 2 +

М

колебательных

I . СО2(01Ю) + Х=С02(00°0) + X

II.С02(10°0) + Х =С02(02«0) + X III. С02(00°1) + Х=С02(І1Ю) + X IV. С02 (00<Ч) + Х=СО2(30Ю) + X

V. С0(0) + С02(00°1)=с0(1) + С02(00°0)

VI.

СО(1) +

С02(00°0)=СО(0) +

С02(11Ю)

V II.

N2 (0) +

CO2(0001)=NÜ(1) +

С02 (0000)

V III.

N2(l) +

CO(0)=N2(0) + СО(1)

IX.

N0(0) + CO(l)=NO(l) + С0(0)

X.

N2(l) +

C02(00«0)=N2 (0) +

С02(11Ю),

г д е X = N 2 , C 0 2 , C 0 .

Релаксационные процессы в. струях низкотемпературной плазмы 201

Рис. 5. Колебательная неравновесность ширеиии в области I ( М а — 1)

а — распределение колебательных темпе­ ратур It = I см, Ткр = 3000° К , РКр -

=40 атплс; пунктирные линии соответствуют расчетам при вероятностях, полученных из теории Шина и SSH; сплошные линии — при расчете вероятностей процессов II—IV с учетом ангармоничности и вращения

б — распределение инверсной заселенно­ сти уровней СОі AN = N (001) — N (010>

вдоль потока при расширении смеси в ко­ ническом сопле. Іо = 0,2 см, Ркр = 40 атм, ТКр - 1500 °К

Константы скоростей химических реакций были взяты из [13]. Вероятности колебательных процессов II и III рассчитывались в [8] с учетом ангармоничности и Кориолисова взаимодействия [15]

 

3000° К

2000° К

800° К

600° к

Q J J

2.10-а

1,3-10-2

8.10-3

6-10-3

Q n l

2,4-Ю -з

1 ,0 -Ю-з

4 . 1 0 -4

1 , 2 -1 0 -«

Вероятности остальных колебательных процессов рассчитыва­ лись по теории S S H и Шина [16]. Уравнения (8) — (10) интегриро­ вались методом Рунге — Кутта. Параметры в выходном сечении звукового сопла варьировались в диапазоне Р а = 50 1 атмг Т а = 2000 -г- 3000° К, га = 0,1 -г-10 см. Результаты расчета по­ казали (рис. 5), что замораживание температуры антисимметрич­ ного валентного колебания происходит на меньших расстояниях от среза, чем деформационного и валентного симметричного ко­ лебаний.

Увеличение концентрации азота приводит к повышению замо­ роженных колебательных температур. Использование вероятно­ сти процесса III из работ [15, 16] приводит к более высоким значе­ ниям замороженных колебательных температур антисимметрично­

2 0 2 Н. И. Ющенкова

го валентного колебания, что указывает на существенное влияние Кориолисова взаимодействия на релаксационные параметры смеси С02—N2. При превышении ТѵС02 (ѵ3) более чем в 1,7 раза Тѵ С02

(ѵ2) в С02 образуется инверсная

заселенность уровней Дп =

= п (001) — п (010). Зависимость

Ап от концентрации N2 носит

немонотонный характер (рис. 5,6). Расчеты Ап с учетом Кориоли­ сова взаимодействия для процесса III показывают снижение Ап но сравнению с расчетами, при которых вероятность процесса III рассчитывалась по методике [16].

В смесях, содержащих атомы и радикалы, образующиеся при введении Н20 в дозвуковую часть течения, кинетика химических реакций приводит к образованию высоких (неравновесных) кон­ центраций Н, ОН [17].

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

Поскольку

процесс диффузии начинает существенно влиять

на состав

газа только в областях IV —V [1], то предполагалось,

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

Параметры течения в изобарическом участке струи были опре­ делены в [19] путем численного решения уравнений пограничного слоя, описывающих усредненные по времени свойства течения в изобарическом участке струи. Эти уравнения записываются в ко­ ординатах Мизеса х, ф, где ф — функция тока, связанная с декар­ товыми цилиндрическими координатами формулами

ЭТ

= Р и«,

 

5Y

рѴу.

 

(14)

^

 

 

 

При этом уравнение неразрывности удовлетворяется автомати

чески, а остальные уравнения имеют вид:

 

 

 

уравнение движения

 

 

/>ТРиУг

dU \

 

 

(15)

dU

 

1

д

 

дх

=

Y

dY

\

Y

3Y /

 

уравнение энергии

 

 

 

 

 

 

 

 

 

дТ

1

д

[7 V

ХН-тР^2

дт ]

 

>/ дх =

ЧГ ЭТ ІА Рг /

^

 

+

 

Р ТР V у 2 ( dU V

 

Le

 

 

 

 

 

1

( і б )

\ а т /

+

Рг

\ 4 * сРі дчг )

dir J

 

pи 4*

 

 

Релаксационные процессы в струях низкотемпературной плазмы 203

уравнение неразрывности z-й компоненты смеси

(17)

Здесь U — составляющая скорости потока, параллельная оси струи; у — координата, отсчитываемая поперек струи; рт — коэф­ фициент турбулентной вязкости; cvf — «замороженная» тепло­ емкость смеси; Le — число Льюиса; Рг — число Прандтля; срі — теплоемкость і-й химической компоненты смеси; а* — массовая скорость образования г'-й компоненты.

Для учета процесса переноса была принята модель турбулент­ ной вязкости, полученная из эксперимента [17]:

И'т^^Рь W b- Uo I,

где г — полуширина струи; pb, Ub — плотность и

скорость на

границе струи; U0 — скорость на оси струи; 0 ,0 1

4 ^ x ^ 0 ,0 3 .

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

1)

ОН + СО =г= Н + С02

7)

о + Н2 =^= 20Н

2) ОН + Н2 =г= Н + Н20

8)

2N + M=N2 + М

3) Н -f 02 =^= ОН + О

9)

О -г О + М =02 + М

4)

Н + О Н ^ Ш + О

10) Н + ОН + М=Н20 + М

5) О + СОз ^ СО -f О2

И) Н + Н + М=Н2 + М

6) о + со+м^со2+м

 

 

Выбор реакций основан на приведенных в [12] механизмах окисления СО и Н2 и возможных реакций обмена с участием ради­ калов. Предполагалось, что в турбулентном потоке можно исполь­ зовать те же температурные зависимости констант скоростей реак­ ций, что и в ламинарном. Константы скоростей прямых реакций были взяты из справочника [13].

Параметры в начале изобарического участка принимались из работы [18]: Т0 = 1600° К, Ѵ0 = 3000 м/сек, радиус струи г0 варь­

ировался в диапазоне от 2- ІО2 до 103 см, 1 ^ Р 0 ^

Ю2 mop.

Пара­

метры внешней среды: Р ь = Р 0, Тъ = 220° К и

Ѵь = 800

м/сек,

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

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

сивно.

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

204

п т ,см -з

Н. И. Юшенкова

Рис. 6. Влияние догорания на распределение температур и состава продуктов сгора­ ния углеводорода в изобари­ ческом участке недорасшиРенной струи

1 Рі = 100 тор,

Го = 200 см, х = 0,104, Т ьіТ0 = 0,2,

ѵ,/ѵв= 2;

2 Р1 — 100 mop,

r0 =

200

см,

X =

0,03;

 

3

— Рі =

10 mop,

r0 =

10 см,

X

=

0,014;

d

Рі =

10 mop,

r0 =

200 см,

X

=

0014;

 

5

— расчет Т в струе без

 

 

учета химических реак­

 

 

ций

 

параметры струи, то при расчетах варьировались: характерный размер — начальный радиус струи, отношение скорости внешнего потока к скорости на оси струи, значение температуры внешнего потока и величина турбулентной вязкости. Расчеты показали, что основное влияние на поле температур оказывает соотношение вре­ мени релаксации и газодинамического для процессов 1) — 4) и 10) (тю^>Ті,4) Во всей области изобарического течения концент­ рации радикалов ОН’ и Н’ превышают равновесные значения, а кон­ центрации СО и Н2 резко падают в зоне догорания. Основными реакциями, влияющими на состав газа, являются обменные реакции 2)—4), а увеличение температуры в струе связано с образованием Н20 в результате реакции 10). При низких давлениях, когда вре­ мя релаксации процесса 10) больше газодинамического, концент­ рация радикалов ОН’ продолжает возрастать в области изобари­ ческого течения, что позволяет применить сверхзвуковые недорас­ ширенные струи для изучения кинетики процессов с участием радикалов при введении примесей в изобарический участок струи.

Релаксационные процессы в струях низкотемпературной плазмы 205

НЕРАВНОВЕСНЫЕ ПАРАМЕТРЫ СТРУЙ ПЛАЗМЫ ИНЕРТНЫХ ГАЗОВ

Определение неравновесных электрофизических параметров (концентраций и температур электронов) в струях плазмы инертных газов с примесями молекулярных газов и инверсной населен­ ностью возбужденных уровней атомов и ионов основывалось на экспериментальных исследованиях струй в плазменных газодина­ мических установках [5, 6, 8, 20]. Теоретические значения населен­ ностей возбужденных уровней оценивались из анализа кинетиче­ ских уравнений с учетом механизмов возбуждения атомов, одно­ кратно заряженных ионов, рекомбинации однократно и двухкрат­ но заряженных ионов с электронами, а также процессов высвечи­ вания. В качестве рабочего газа рассматривался химически чистый аргон. Основные процессы возбуждения и рекомбинации в аргоне связаны с взаимодействием атомов и ионов с электронами.

Для области свободного расширения, в которой распределение плотности и концентрации электронов (при а е = const) описы­ вается формулой (1), резкое падение Те вдоль оси струи и высокие относительные концентрации электронов и ионов (двухзарядных и однозарядных) создают условия для образования неравновесных концентраций возбужденных атомов и ионов. Для оценки концент­ рации возбужденных частиц применен метод работы [4]. Прибли­ женно распределение концентраций возбужденных частиц вдоль оси плазменной струи при заданных Те (I) и пе (Г) описывается за­ висимостью:

Ех Е ь — разность потенциалов ионизации и энергии возбуж­ дения; Те, Тео — температура электронов в струе и на срезе сопла; пге — масса электрона; z — статистическая сумма.

Выполненные оценки распределения концентраций возбужден­ ных частиц в рекомбинирующей плазме аргона (рис- 7) показыва­ ют, что концентрации N* (АІІ4р) крайне чувствительны к рас­ пределению Те и очень сильно отличаются от равновесных. Рас­ пределение пе и релаксация Те в струях исследовались экспери­ ментально. Эксперименты проводились в аэродинамической уста­ новке низкой плотности, описанной в [6, 8, 10]. В качестве гене­ ратора плазмы использовался плазмотрон постоянного тока со стабилизацией дуги газовым слоем [19]. Струя истекала через

206

Н. II. Ющенкова

Рис. 7. Влияние режима истечения на параметры струи

а— распределение концентраций и температур электронов и заселенностей уровней АгІІ4р вдоль струи плазмы аргона, истекающей из звукового сопла в неподвижную среду аргона. Режим] 1 Р0 = 40 тор, I = 500 а, V — 15 в (чистая струя аргоновой

плазмы),ссе ж 1; режим 2 — Р„ = 150 тор, I = 200 а, V = 18 в, Рь = 6-10—®тор; че 10-2 (чистая струя слабоионизирующего аргона)

б — интенсивность излучения плазмы аргона АгІІ на переходах 4р — 4s в зависимости от Р0при I = 500 а, V = 15 в для области струи, отстоящей от сопла на 1 см

сопло в вакуумную камеру с давлением Рь = ІО-2 — 10-1 тор. Давление Р 0 в плазмотроне (на входе в сопло) было существенно

выше Рь. Параметры истечения изменялись в диапазоне Р0 =

ІО2 -f- 10 тор, М а = 8 - і- 1 .

Мощность, вкладываемая в разряд,

составляла 5 —20 кет при I

= 100 -ь 1000 а.

В качестве рабочих газов использовались спектрально чистые

аргон и ксенон, в вакуумную камеру вводились молекулярные

газы SF6, CC12F2, СС14, С02.

 

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

Эксперименты показали, что режим истечения существенно влияет на физико-химические параметры, излучение и структуру

плазменной струи. При режиме 1 соответствующему /

200 а,

Р о )> 100 тор, струя плазмы аргона имеет розовый цвет и ясно различимую волновую газодинамическую структуру (см. рис. 3, б). Спектр плазмы (рис. 8, а) содержит в основном линии возбужден­ ных атомов. При уменьшении давления до 60—30 тор и увеличе-

Релаксационные процессы в струях низкотемпературной плазмы 207

нии тока дуги / 500 а степень ионизации повышается. При этом плазма приобретает голубой цвет (рис. 8, б), в спектре исче­ зают линии атомов и резко увеличивается интенсивность линий возбужденных ионов, соответствующих переходам 4р — 4s (см. рис. 8, а, режим 2).

При переходе к бесстолкновительному режиму истечения (в ре­ зультате снижения P L) качественно изменяется структура струи плазмы за счет влияния магнитогазодинамических эффек­ тов. Наблюдается (рис. 8, в) образование плазменного «шнура», сжимающегося под воздействием собственного магнитного поля. Параметры шнура характеризуются концентрациями nt = ІО13 -г-

1014 смГ3и температурой ионов порядка 10 эв.

Для режимов истечения, при которых в спектре струи наблю­ дались интенсивные линии излучения однозарядных ионов, соот­ ветствующих переходам пр ns, в вакуумную камеру вводились молекулярные газы (SF8, N2, CC12F2, С02), которые вследствие раз­

реженности плазменного потока при (^кр/^кр) Ра/^ь ІО-2 диф­ фундировали во все области струи. При этом происходили изме­ нения в электрофизических параметрах и излучении струй. Ви­ зуально наблюдалась резкая граница светящегося плазменного образования с формой, близкой к эллипсоиду и ярко-синей окрас­ кой для струй аргона (рис. 8, г) и зеленой —для струй ксенона (рис. 8, д).

Зондовые

и спектроскопические

исследования параметров

струй

показали

уменьшение

Те в

областях

II — IV, сни­

жение

пе в

области за границей образования и резкое увеличе­

ние

интенсивности

излучения

ArII,

ХеІІ из

данной области

[А, =

4376, 4880, 4965 Â — для аргона (рис. 8, е) и 4603, 4683, 4983

4919, 5185 Â — для

ксенона]. Увеличение интенсивности излуче­

ния может быть связано как с изменением Те, так и с процессом нейтрализации ионов инертных газов ионами электроотрицатель­ ных газов или механизмом резонансной перезарядки. Следова­ тельно, если плазма содержит высокие концентрации двухзаряд­ ных ионов (что реализуется при электронных температурах на срезе сопла 20000—30000° К), то в области сверхзвукового рас­ ширения образуются высокие концентрации возбужденных метастабильных однократно заряженных ионов. Последнее позволяет применять струи плазмы для проведения кинетических исследо­ ваний взаимодействия возбужденных заряженных частиц с моле­ кулярными газами, вводимыми в струю плазмы.

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

станием температур проводилось за системой

ударных

волн

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

показали,

что в

Соседние файлы в папке книги из ГПНТБ