Например, Бобцов

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

КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ПРИ ИЗУЧЕНИИ ФИЗИЧЕСКИХ...

УДК 537.5 КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ПРИ ИЗУЧЕНИИ ФИЗИЧЕСКИХ ПРОЦЕССОВ В ТЛЕЮЩЕМ РАЗРЯДЕ В ВОЗДУШНОЙ СМЕСИ ПРИ НИЗКИХ ДАВЛЕНИЯХ М.В. Чернышеваa, В.П. Марекb, А.С. Чирцовc, Д.А. Швагерc
a Институт фундаментальной электроники, Париж, Франция, Maria.Tchernycheva@ief.u-psud.fr b Эколь - Политехник, Париж, Франция, nika.marek@mail.ru c Университет ИТМО, Санкт-Петербург, Россия, alex_chirtsov@mail.ru
В соответствии с запланированным циклом исследований по развитию методов численного моделирования многокомпонентных газоразрядных сред, основанных на современной концепции нелокальной плазмы, построены простые тестовые численные модели газового разряда в воздухе при пониженных давлениях. В рамках полуэмпирического подхода рассчитана зависимость от давления падения напряжения на разрядном промежутке как суммы падений в прикатодном слое и на положительном столбе тлеющего разряда, вычисляемых в соответствии с их моделями, учитывающими ионизацию электронным ударом из основного состояния, рекомбинацию носителей на стенках разрядной трубки, ионно-электронную эмиссию с катода и потерю энергии электронами при столкновениях с невозбужденными молекулами. Рассчитанные зависимости падения напряжения от давления сопоставлены с результатами численного моделирования разряда, выполненного на основе решения системы гидродинамических уравнений для концентраций в двумерном случае при учете перечисленных процессов. В рамках вышесказанных приближений обе модели привели к удовлетворительно согласующимся между собой результатам, которые оказались существенно заниженными по сравнению с экспериментальными результатами в соответствующей нормальному тлеющему разряду области высоких давлений. Полученные расхождения могут объясняться недостаточной полнотой системы элементарных процессов, учитываемых при моделировании. Например, неучтенное присутствие в качестве примесей газов из трехатомных молекул может приводить к увеличению столкновительных потерь энергии электронов, что, в свою очередь, должно приводить к увеличению напряженности электрического поля в положительном столбе по сравнению с полученными в моделях значениями. Ключевые слова: тлеющий разряд, численное моделирование, положительный столб, прикатодный слой, нормальный разряд, аномальный разряд, полуэмпирическая модель, двумерная модель.

COMPUTATIONAL MODELING IN THE STUDY OF GLOW DISCHARGE PHYSICAL PROCESSES IN THE AIR AT LOW PRESSURES
M.V. Tchernychevaa, V.P. Marekb, A.S. Chirtsovc, D.A. Shvagerc
a Institut d'Electronique Fondamentale Paris, France, Maria.Tchernycheva@ief.u-psud.fr b Ecole Polytechnique, Paris, France, nika.marek@mail.ru с ITMO University, Saint Petersburg, Russia, alex_chirtsov@mail.ru
Simple numerical models of low-pressure air gas discharge were developed based on modern concepts of nonlocal plasma methods of the computational simulation of multi-component discharge environments. Voltage drop at the discharge gap was calculated as a sum of incidences in the cathode sheath and on the positive column of a glow discharge as part of the semiempirical approach. Electron impact ionization of the ground state, the recombination of charged particles on the walls of the discharge tube, ion-electron emission from the cathode, and the energy loss of electrons in collisions with ground state molecules have been taken into account in the models. The calculated voltage drop, depending on the pressure, has been compared with the results of another numerical discharge model based on the solutions of hydrodynamic equations for the particle densities in a two-dimensional case. Within these approximations, both models give satisfactory agreement between their results, which are strongly underestimated compared to the experiment corresponding to the normal high-pressure glow discharge. The resulting discrepancy can be explained by an incompleteness of elementary processes taken into account in the simulation. For example, an unaccounted presence of gases from triatomic molecules as an impurity can lead to an increase in the collision energy losses of electrons, which in turn should lead to an increased electric field in the positive column, as compared with the values obtained in the models. Keywords: glow discharge, computational modeling, positive column, cathode sheath, normal discharge, abnormal discharge, semi-empirical model, two-dimensional model.
Введение
Изучение газовых разрядов в воздушных смесях актуально с точки зрения многочисленных приложений (в том числе – интенсивно развивающихся современных плазменных технологий) и одновременно представляет фундаментальный интерес как удобный для эксперимента, но весьма сложный для анализа тестовый объект для приобретающих все большую популярность методов численного моделирования газовых разрядов [1–3], основанных на использовании недавно выдвинутой и интенсивно развивающейся концепции нелокальной плазмы [4]. В группе разрядов в воздухе на современном этапе наибольший интерес вызывают разряды при атмосферных давлениях (см., например, подборки публикаций [5, 6]). В связи с развитием нанотехнологий усилился интерес к микроразрядам, возникающим при нормальных давлениях в коротких межэлектродных промежутках [7–9]. Поэтапный переход к атмосферным микроразрядам от популярных для экспериментальных исследований середины ХХ века разрядов при пониженных давлениях [10] важен для работ по созданию газовых микродетекторов, основанных на принципах электронной столкновительной спектроскопии [11].

140

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics
2014, №3 (91)

М.В. Чернышева, В.П. Марек, А.С. Чирцов, Д.А. Швагер

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

Методы моделирования газоразрядных сред

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

его моделировании на «усредненно-микроскопическом» уровне. Подразумевается численное решение

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

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

стей чисто вычислительного характера (моделирование макросистемы с большим числом подансамблей

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

блем. К ним относятся выбор алгоритмов ограничения набора учитываемых процессов и числа кванто-

вых состояний частиц плазмы и задание значений констант скоростей kp+q = j+{j'} (или соответствующих сечений) плазменных реакций рождения частицы определенного сорта в квантовом состоянии j и прочих

продуктов {jʹ} в результате столкновения двух частиц p и q (в используемом описании одинаковые части-

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

Отдельную проблему составляет подбор данных о вероятностях элементарных процессов. В случае

столкновений атомов и молекул, как правило, безальтернативным источником являются банки экспери-

ментальных данных (например, [13]). Вероятности элементарных процессов с участием бесструктурных

частиц (фотонов и электронов) могут рассчитываться становящимися все более доступными по ресурсо-

емкости квантовомеханическими методами [14]. Однако в важнейших для плазменных процессов случаях

столкновений атомов с электронами традиционно используемое борновское приближение [15] и его мно-

гочисленные модификации [15–17] нередко дают расхождения с результатами измерений в два и более

раз [18, 19].

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

для концентраций частиц рассматриваемого типа в заданном квантовом состоянии nj:


   


n

( j

t

)



k () p  q  j { j}

np nq



k () j  p { j}

njnp



p,q, j

p, j'

,

J

() j

,

J

( j

)





D

j

n

( j

)





j

q

j

n

( j

)

E

,

(1)

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

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

E   ,   4 qj nj , j

(2)

и аналогичными по форме (1) уравнениями баланса для энергий.

Подходы к численному решению системы (1)–(2) во многом зависят от того, какая модель плазмы

(одно-, двух- или трехмерная) используется в рассматриваемой задаче. В одномерном случае для про-

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

тические подходы. Модели, претендующие на более адекватное и приближенное к реальности описание,

столь сложны, что для выполнения соответствующих расчетов приходится использовать профессиональ-

ные пакеты численного моделирования (например, COMSOL [20]) и системы высокопроизводительных

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

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

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

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

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

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

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

числе, трехатомных молекул. Перечисленные причины приближают задачи микромоделирования воз-

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics 2014, №3 (91)

141

КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ПРИ ИЗУЧЕНИИ ФИЗИЧЕСКИХ...

душных разрядов к границам современных возможностей физики низкотемпературной плазмы и вычислительных технологий. Построение адекватной численной модели такой системы можно рассматривать как одну из итоговых задач-максимум программы исследовательской активности в области изучения нелокальной плазмы на ближайшие годы. Это соображение, наряду с практической значимостью воздушных разрядов, предопределило его выбор в качестве тестового объекта для работ по развитию, основанных на концепции нелокальной плазмы, методов численного моделирования разрядов в многокомпонентных смесях.
Задача построения еще более приближенных к реальности (не использующих усреднение) численных микромоделей открытого неравновесного ансамбля частиц многокомпонентной плазмы (основанных, например, на методе Монте-Карло) на современном этапе представляется малоперспективной.
Вместе с тем, до сегодняшнего дня при решении прикладных задач остаются весьма популярными развитые в ХХ веке сравнительно простые полуэмпирические методы моделирования газоразрядных сред, основанные на аналитических подходах [10, 12, 21]. В их рамках специфика газовых сред учитывается введением эмпирических подгоночных параметров, получаемых из сравнения с экспериментом. Такие методы основываются на заведомо упрощенном и усредненном физическом описании различных областей разряда и, очевидно, способны давать оценки его усредненных характеристик, нередко получаемых в рамках прикладных экспериментальных исследований. Именно такие характеристики были получены в [22] при измерениях электрических характеристик тлеющего разряда в воздухе при пониженных давлениях. Это позволило осуществить оценку применимости полуаналитических моделей для интегрального описания тлеющего разряда в сложной многокомпонентной газовой смеси, необходимую на начальном этапе выполнения работ по развитию более совершенных моделей, подразумевающих переход к микроописанию совокупности элементарных процессов в разряде.

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

Основной вклад в измеряемое на опыте падение напряжения на межэлектродном зазоре традици-

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

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

Описание электрических свойств положительного столба строится на основе его макроскопиче-

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

ва энергий (получаемой электронами от электрического поля и отдаваемой ими частицам плазмы при

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

столбе напряженность электрического поля E с электронной температурой Te:

3 2

Te





 0,8 qe 

E N

.

(3)

Входящий в уравнение (3) параметр δ определяет долю энергии электрона, отдаваемую им при

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

элементарной теории упругих столкновений. Столкновение же с молекулами требует нетривиального

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

тате отдельно выполняемого численного моделирования. В нашем случае были использованы экспери-

ментальные данные [23–25].

В свою очередь, входящая в (3) электронная температура связана с частотой ионизации электрон-

ным ударом Zi(Te) и поперечными размерами трубки R. Эта связь следует из упрощенного (одномерного) варианта уравнения (1), в котором учитываются лишь доминирующий процесс прямой ионизации невоз-

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

d 2ne dr 2



1 r

dne dr



Zi

Тe 
D

ne



0

.

(4)

В соответствии с распределением Саха частота ионизации экспоненциально зависит от отношения

потенциала ионизации Ui к энергии электронов [26]:

Zi



N

6 102 

  

2kTe me

3/ 2  

me qe

1  

qeUi kTe

 

exp

 





qeUi kTe

  

.

Совместно с граничным условием ne(r=R) = 0, соответствующим гибели электронов на стенке, уравнение Бесселя (4) связывает Те с радиусом R разрядной трубки. Необходимая для вычисления напряженности поля электронная температура определяется либо в результате численного решения всей

описанной одномерной задачи (с возможным уточнением граничных условий на стенке [10]), либо при

решении содержащего зависящий от типа газа экспериментально подбираемый параметр C трансцен-

дентного уравнения

142

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics
2014, №3 (91)

М.В. Чернышева, В.П. Марек, А.С. Чирцов, Д.А. Швагер

kTe qeU i

1  

1 2

qeU i kTe

 

exp

 





qeU i kTe

1,16 107 

CpR 2

 1,

(5)

возникающего в случае использования известного для нулевых граничных условий решения (4) в форме

R (Zi / D)1/2 =2,4. Получаемые в результате описанных процедур зависимости электронной температуры и напря-

женности поля приведены на рис. 1, a, б. Падение напряжения на столбе вычисляется умножением рас-

считанной напряженности поля Е на экспериментально измеренную для тех же давлений длину положи-

тельного столба [22].

Те, эВ Е, В/см

1,5 1,0 0,5

2р, мм Hg 4 а

Аномальный разряд
2

u(p)

1

1,5 1,0 0,5
24 р, мм Hg б

8

6 d~ 4

~E ~V

2 е–1
10–2

100 в

102 ~j

Нормальный разряд

Аномальный разряд

Нормальный

u(ln p)

разряд

2

1

0

0

2

4

6 –6

–4 –2

0

2

гд

Рис. 1. К расчету падения напряжений u на положительном столбе и в катодном слое от давления p: полученная в результате численного решения (5) зависимость Te(p) (а); рассчитанная зависимость напряженности поля в положительном столбе от давления (б); полученная в [10] зависимость толщины
катодного слоя d , поля на катоде E и падения катодного напряжения V в безразмерных
величинах (в); полученные в рамках расчетов зависимости катодного падения напряжения (в кВ) в трубке от давления воздуха (в мбар) в ней в обычном (г) и логарифмическом (д) масштабах

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

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

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

дели однородного прикатодного слоя предложен в [10] на базе известного из элементарных курсов усло-

вия самоподдержания разряда:

 d  ln 1 1 ,

(6)

связывающего между собой толщину прикатодного слоя d, ионизационный коэффициент α и параметр γ, характеризующий вероятность выбивания электрона с катода в результате прихода на него иона. При использовании для α эмпирической формулы Таундсена модель позволяет связать плотность тока на поверхности катода с падением напряжения в прикатодном слое и оценить минимальное падение напряжения в слое и соответствующее ему «нормальное значение» плотности тока на катоде jn (рис. 1, в). Произведение этой нормальной плотности тока на площадь катода дает предельную величину разрядного тока, до достижения которой падение напряжения в прикатодной области остается постоянным («нормальный разряд»). При дальнейшем увеличении тока падение напряжения в слое начинает резко возрастать, что соответствует переходу к «аномальному» режиму разряда.
Для тестовой оценки падения напряжения в катодном слое взяты результаты, полученные в [10] в рамках простейших приближений (6). Выполненные с их использованием расчеты соответствующих плотностей тока на допускающем существование тлеющего разряда интервале давлений показали, что на реализованном в [22] интервале условий происходил переход от нормального типа разряда к аномальному (рис. 1, г, д).

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics 2014, №3 (91)

143

КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ПРИ ИЗУЧЕНИИ ФИЗИЧЕСКИХ...

Сделанные на основе [10] первоначальные тестовые оценки были уточнены в результате естест-

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

трического поля. В нем был использован очевидный аналог критерия самоподдержания разряда (6) для

случая зависимости напряженности поля от расстояния от катода:

d
  E  zdz  ln 1 1  .

(7)

0

Для E(z) зависимость выбиралась либо линейной, либо соответствующей данным зондовых изме-

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

значениям катодного падения напряжения были не более 10%.

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

моделирования разряда с параметрами, соответствующими эксперименту [22], вместе с эксперименталь-

но полученными данными.

6

Напряжение на трубке, кВ

5

4

Аномальный разряд

u(ln p)

Нормальный разряд

36

24

12

0 1 2 3 4 5 6 –6 Давление в трубке, мбар

–4

0 –2 0

2

аб
Рис. 2. Сравнение результатов численного моделирования при полуэмпирическом подходе (маркеры в виде ромбов) с экспериментально полученными зависимостями падений напряжений на трубке от давления (откачка трубки – маркеры в виде квадратов, напуск воздуха – треугольные маркеры) для питающего напряжения 5,1 кВ в обычном (а) и логарифмическом (б) масштабе

Существенно заниженные по сравнению с данными эксперимента результаты могут быть исправлены в результате варьирования подгоночных параметров полуэмпирической модели, полученных при измерениях [23–25] в существенно отличных от [22] условиях. Однако более привлекательным является подход, основанный на априорном улучшении физической модели путем явного учета существенно большего количества плазменных процессов, чем в (3)–(7). Такой подход требует перехода к микроскопическому описанию газоразрядных сред.

Переход от полуэмпирических макроскопических моделей газоразрядной плазмы к микромоделям

В качестве первого шага на пути построения сглаженной микроскопической модели разряда в воздушной смеси необходима апробация гидродинамического подхода на модельной системе, физические свойства которой максимально соответствуют приближениям, использованным в полуаналитических расчетах по методике (3)–(7).
В практике изучения воздушных разрядов и предварительных оценок их параметров весьма популярно использование его упрощенного аналога – разряда в аргоне. Возможность замены оправдывается близостью величины энергии ионизации аргона к мало различающимся между собой значениям для азота и кислорода. Учет в полуэмпирическом подходе только двух механизмов образования свободных электронов в разряде (ионизация молекул из основного состояния электронным ударом и выбивание электронов с катода ионами) делает удобным использование описанного в [20] стандартного примера моделирования разряда в аргоне для создания сглаженной микроскопической модели системы, соответствующей приближениям (3)–(7). В результате были выполнены расчеты распределения потенциала вдоль оси цилиндрической трубки и интегрального падения напряжения на межэлектродном участке при различных давлениях газа в рамках существенно двумерной численной модели разряда в аргоне, обладающего цилиндрической симметрией, параметры которого соответствовали использованным при полуэмпирических расчетах воздушного разряда (рис. 3).
Исходная модель разряда в аргоне учитывала следующие компоненты газоразрядной среды: атомы в невозбужденном состоянии (Ar), атомы в возбужденном долгоживущем метастабильном состоянии (Ar*), однократные ионы аргона в основном состоянии (Ar+) и электроны (е). В результате, помимо аналогов рассматриваемых в (3)–(7) процессов (прямой ионизации из основного состояния электронным ударом и ухода заряженных частиц на стенки), список реакций, протекающих в объеме,

144

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics
2014, №3 (91)

М.В. Чернышева, В.П. Марек, А.С. Чирцов, Д.А. Швагер

Ar  e  Ar  e , Ar  e  Ar * e , Ar * e  Ar  e ,

Ar  e  Ar  2e ,

(8)

Ar *e  Ar  2e ,

Ar *Ar*  Ar  Ar e ,

Ar *  Ar  Ar  Ar ,
содержал процессы с участием возбужденных атомов, которые, при прочих равных условиях, способны облегчать ионизацию и приводить к уменьшению эффективного сопротивления разрядного промежутка. Другой очевидной причиной, способной привести к аналогичным расхождениям результатов, получаемых по сопоставляемым моделям (3)–(7) и (8), является существенное различие в вероятностях γ выбивания с катода электронов ионами аргона и азота (соответственно 0,09 [27] и 0,005 [28]). В результате расчеты по схеме (8) дали еще более заниженный результат (рис. 3, а, кривые 2 и 3).
Исключение из рассмотрения реакций с участием атомов аргона в метастабильных состояниях путем искусственного уменьшения константы скорости их возбуждения на 1–2 порядка не привело к значительным изменениям интегральных электрических свойств модельного разряда, а корректировка эффективности выбивания электронов с катода повысила величину падения напряжения в прикатодном слое и привела к хорошему согласию между сопоставляемыми моделями в области аномального разряда, соответствующей падающему участку кривой U(p) (рис. 3, а, кривые 1 и 3). На соответствующем нормальному разряду интервале давлений результаты расчетов по описанной микромодели все еще остаются сильно заниженными.
На рис. 3, б, приведены зависимости распределения потенциала вдоль оси трубки, полученные при усредненном микромоделировании разряда в аргоне (8) в описанном модельном газе, имитирующем воздушную смесь. Как видно, доминирующий вклад принадлежит падению напряжения в прикатодной области, которое в области аномального разряда согласуется с полуэмпирическими расчетами. В области нормального разряда рост падения напряжения с давлением должен обусловливаться увеличением напряженности поля в положительном столбе, для которого микромодель дает заниженные значения. Последнее, по-видимому, связано с некорректным учетом механизмов потери энергии электронами при столкновениях, которая может быть грубо описана с помощью входящего в (3) параметра δ. В случае инертного газа он оценивается отношением масс электрона и атома и имеет величину порядка 10–4, для столкновений с молекулярным азотом значение δ оказывается примерно в 20 раз большим [24, 25]. При условии примерного равенства электронных температур (во всех трех случаях моделирования 1–1,5 эВ) в соответствии с (3) корректный учет потерь энергии электронами должен привести к существенному возрастанию напряженности поля (и, следовательно, падения напряжения) в положительном столбе примерно в 4 раза. В результате значения падения напряжения в области нормального разряда должны совпадать с результатами аналитической модели или даже превышать их. Последнее, разумеется, требует проверки с использованием микромодели тлеющего разряда в молекулярном азоте и его смеси с кислородом, а такжесопоставления результатов моделирования с экспериментом.

Напряжение на трубке, кВ Электрический потенциал, В

2000 1600 1200 800 400
0

1 3 2
1234567 Давление в трубке, мбар

700 600 500 400 300 200 100
0
0

Зависимость электрического потенциала от расстояния
0,1 0,2 0,3 0,4 0,5 0,6 Расстояние вдоль оси трубки, мм



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

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics 2014, №3 (91)

145

КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ПРИ ИЗУЧЕНИИ ФИЗИЧЕСКИХ...
Экспериментальное сопровождение работ по развитию методов микромоделирования разрядов в многокомпонентных газовых смесях
Первые оценочные сопоставления основанных на двух подходах тестовых моделей разрядов в воздухе с результатами экспериментальных измерений проводились по интегральным электрическим характеристикам разряда. Данные были получены на лабораторной установке, созданной на базе промышленно выпускаемого комплекта оборудования [29] и предназначенной для эксперимента по организации физического учебного практикума нового типа, включающего в себя в качестве обязательной составляющей научную деятельность обучаемых [30–32].
Тлеющий разряд зажигался в межэлектродном промежутке длиной 500 мм в трубке диаметром 50 мм, откачиваемой до давления от нескольких торр до сотых торра, от источника, подключаемого через балластное сопротивление в 5 Мом и регулируемого на интервале от 0 до 10 кВ. При питающих напряжениях 3–10 кВ наблюдался устойчивый тлеющий разряд с четко выраженной прикатодной областью и стратифицированным положительным столбом. Перечисленные параметры были использованы при построении рассмотренных моделей.
Полученные зависимости интегральных характеристик разряда от давления и питающих напряжений свидетельствуют о систематическом занижении рассчитываемых в рамках тестируемых методов полей и межэлектродных падений напряжения в тлеющем разряде в воздухе (рис. 4). Одной из главных причин может быть несоответствие реальных потерь электронами их энергии при столкновениях с используемым в моделях описанием из-за присутствия в реальной воздушной смеси примесей из трехатомных молекул (воды, углекислого газа) с богатыми колебательными и вращательными энергетическими спектрами. Для ответа на поставленный вопрос необходимы измерения в чистых газах, а также в смесях, содержащих калиброванные добавки. Кроме того, детальный характер описания газоразрядных сред, даваемый микромоделями, ставит вопрос о возможности получения экспериментальных данных с той же степенью полноты. Указанные соображения инициировали работы по созданию универсального измерительного стенда, включающего вакуумную камеру, допускающую зажигание разрядов с различной геометрией в газовых смесях с различными калиброванными составами и позволяющую осуществлять селективную по разрядному объему комплексную зондовую и оптическую (на следующем этапе) диагностику плазмы.
Разрядная трубка

Вакуумметр

аb

Высоковольтный источник питания

Система откачки
а

б

Напряжение на трубке, В

2000 1600 1200 800 400

0 0,5 1 1,5 2 2,5 3 3,5 4 Давление в трубке, мбар

Данные первого эксперимента Полуаналитическая модель

Разряд в аргоне

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

в
Рис. 4. Стандартная комплектация лабораторной установки для изучения газового разряда: схема учебной установки (а); общий вид установки (б); сопоставление данных, полученных на установке, с результатами
численного моделирования (в)

146

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics
2014, №3 (91)

М.В. Чернышева, В.П. Марек, А.С. Чирцов, Д.А. Швагер

Систематические измерения характеристик газовых разрядов в различных смесях и ресурсоемкие расчеты подразумевают выполнение больших объемов сравнительно однотипных исследовательских работ, к которым могут привлекаться наиболее подготовленные и мотивированные учащиеся из различных возрастных групп. Указанная деятельность может рассматриваться как своеобразный методический эксперимент по решению широко обсуждаемой сегодня задачи раннего привлечения молодежи к научным исследованиям в целях обеспечения непрерывного перехода от обучения к исследовательской деятельности.
Заключение
Результаты выполненного в рамках популярных сегодня полуэмпирических методов моделирования тлеющего разряда в длинной трубке с воздушной смесью при пониженных давлениях привели к существенно заниженным, по сравнению с данными эксперимента, значениям падения напряжения на межэлектродном промежутке. Полученные систематические расхождения были ожидаемы в связи с заведомо упрощенным характером подхода, основанного на аналитическом описании, и ставят вопрос о возможности построения более точных моделей разрядов в многокомпонентных газовых смесях. Выполненные на базе простейшей модели, использующей гипотетический воздухоподобный газ, тестовые расчеты, основанные на использовании усредненного микроскопического подхода, продемонстрировали качественно правильное поведение численной модели исследуемого разряда. Последнее открывает путь поэтапного развития модели с целью ее приближения к адекватному описанию разрядов в воздушной смеси. Детальное описание разрядов, получаемых в результате использования моделей апробированного типа, требует соответствующего развития используемого измерительного оборудования в направлении существенного расширения номенклатуры получаемых на нем экспериментальных данных о разряде. Соответствующие работы выполняются.
Литература
1. Bogdanov E.A., Kudryavtsev A.A., Ochikova Z.S. Main scenarios of spatial distribution of charged and neutral components in SF6 plasma // IEEE Transactions on Plasma Science. 2013. V. 41. N 12. P. 3254–3267.
2. Богданов Е.А., Капустин К.Д., Кудрявцев А.А., Чирцов А.С. Сопоставление различных вариантов гидродинамического (fluid) моделирования продольной структуры микроразряда атмосферного давления в гелии // ЖТФ. 2010. Т. 80. № 10. С. 41–53.
3. Bogdanov E.A., Chirtsov A.S., Kudryavtsev A.A. Fundamental nonambipolarity of electron fluxes in 2D plasmas // Physical Review Letters. 2011. V. 106. N 19. Art. N195001.
4. Kaganovich I.D., Demidov V.I., Adams S.F. and Raitses Y. Non-local collisionless and collisional electron transport in low-temperature plasma // Plasma Physics and Controlled Fusion. 2009. V. 51. N 12. Art N124003.
5. Современные проблемы физики газового разряда [Электронный ресурс]. Режим доступа: http://physelec.phys.msu.ru /science/discharge/discharge.html, свободный. Яз. рус. (дата обращения 12.11.2013).
6. Студентам и школьникам – книги – газовый разряд [Электронный ресурс]. Режим доступа: http://www.ph4s.ru/book_ph_razryad.html, свободный. Яз. рус. (дата обращения 12.11.2013).
7. Автаева С.В., Оторбаев Д.К., Скорняков А.В. Экспериментальное исследование характеристик тлеющего разряда в воздухе // Вестник КРСУ. 2002. Т 2. № 2. С. 4–6.
8. Mahadevan S., Raja L.L. Simulations of direct-current air glow discharge at pressures 1 Torr: Discharge model validation // Journal of Applied Physics. 2010. V. 107. N 9. Art. N093304.
9. Gambling W.A., Edels H. The high-pressure glow-discharge in air // British Journal of Applied Physics. 1954. V. 5. N 1. Р. 36–39.
10. Грановский В.Л. Электрический ток в газе. Установившийся ток. М.: Наука, 1971. 490 с. 11. Кудрявцев А.А., Мустафаев А.С., Цыганов А.Б., Чирцов А.С., Яковлева В.И. Спектры энергии элек-
тронов в гелии, наблюдаемые в микроплазменном детекторе CES // ЖТФ. 2012. Т. 82. № 10. С. 1–6. 12. Кудрявцев А.А., Смирнов А.С., Цендин Л.Д. Физика тлеющего разряда. СПб: Лань, 2010. 512 с. 13. Compilation of electron cross sections used by A.V. Phelps [Электронный ресурс]. Режим доступа:
http://jilawww.colorado.edu/~avp/collision_data/electronneutral/electron.txt/ свободный. Яз. англ. (дата обращения 14.09.2013). 14. Собельман И.И. Введение в теорию атомных спектров. М.: Физматгиз, 1963, 640 с. 15. Вайнштейн Л.А., Собельман И.И., Юков Е.А. Сечения возбуждения атомов и ионов электронами. М.: Наука, 1973. 142 с. 16. Очкур В.И. О методе Борна–Оппенгеймера в теории атомных столкновений // ЖЭТФ. 1963. Т. 45. С. 753. 17. Гордеев С.В., Чирцов А.С. Столкновительные переходы между различающимися по спину высоковозбужденными уровнями атомов второй группы // Вестник СПбГУ. Серия 4. Физика. Химия. 1991. № 1. С. 146–149.

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics 2014, №3 (91)

147

КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ПРИ ИЗУЧЕНИИ ФИЗИЧЕСКИХ...

18. Марек В.П., Чирцов А.С. Исследование столкновительных переходов с изменением спина между высоковозбужденными состояниями атомов гелия методом лазерной накачки // Известия Международной академии наук высшей школы. 2012. № 2 (60). С. 29–36.
19. Marek V., Chirtsov A. Research of collisional displacement of diverse highly exited states of Helium by means of the method of levels laser injection in cluster and plasma zone // Proc. of Int. Scientific Seminar Physics of Laser Processes and Applications. Ryazan, 2012. P. 96–101.
20. Russian Comsol Page. HUMUSOFT [Электронный ресурс]. Режим доступа: http://www.humusoft.com/produkty/comsol/ru свободный. Яз. рус. (дата обращения 14.09.2010).
21. Райзер Ю.П. Физика газового разряда. М.: Наука, Физматлит, 1987. 592 с. 22. Чирцов А.С., Марек В.П. Варианты использования компьютерных технологий для интенсификации
практикумов и приближения учебных работ к научным исследованиям // Информатика и образование. 2013. № 9. С. 22–34. 23. Ashihara O. The electron energy loss rates by polar molecules // Institute of Space and Aeronautical Science, Tokyo. 1975. V. 40. N 530. P. 257. 24. Dalidchik F.I., Sayasov Yu.S. Recombination of electrons in molecular gases // Journal of Experimental and Theoretic Physics, 1965. V. 49. P. 302–305. 25. Rudge M.R.H. Theory of ionization of atoms by electron impact // Reviews of Modern Physics. 1968. V. 40. P. 564. 26. Энгель А. Ионизированные газы / Пер. с англ. М., ГИФМЛ, 1959. 332 с. 27. Добрецов Л.Н., Гомоюнова М.В. Эмиссионная электроника. М.: Наука, 1966. 564 с. 28. Браун С. Элементарные процессы в плазме газового разряда / Пер. с англ. М.: Атомиздат, 1961. 322 с. 29. Catalogues – LD DIDACTIC [Электронный ресурс]. Режим доступа: http://www.lddidactic.de/en/catalogues.html свободный. Яз. англ. (дата обращения 12.05.2013). 30. Марек В.П., Чирцов А.С. Разработка мультимедийных описаний для нового лабораторного практикума по физике // Материалы XII Межд. конф. Физика в системе современного образования (ФССО-13). Петрозаводск: ПетрГУ, 2013. Т. 2. С. 217–220. 31. Чирцов А.С., Марек В.П. Новые подходы к созданию и использованию мультимедийных ресурсов // Компьютерные инструменты в образовании. 2011. №1. С. 58 –72. 32. Чирцов А.С., Абутин М.В., Марек В.П., Микушев С.В. Новые варианты использования информационных и мультимедийных технологий для реализации непрерывного высшего образования // Физическое образование в вузах. 2012. Т. 18. № 1. С. 109–125.

Чернышева Мария Владимировна Марек Вероника Петровна Чирцов Александр Сергеевич

– – –

Швагер Данила Александрович



работник-исследователь 1-го класса, Институт фундаментальной электроники, Париж, Франция, Maria.Tchernycheva@ief.u-psud.fr студент, магистр, Эколь - Политехник, Париж, Франция, nika.marek@mail.ru кандидат физ.-мат. наук, профессор, Университет ИТМО, СанктПетербург, Россия; доцент, Санкт-Петербургский государственный университет, физический факультет, alex_chirtsov@mail.ru аспирант, Университет ИТМО, Санкт-Петербург, Россия, shvagerdanila@gmail.com

Maria V. Tchernycheva
Veronica P. Marek Alexander S. Chirtsov
Danila A. Shvager

– Chargé de recherche CNRS 1 classe, PhD, Institut d'Electronique
Fondamentale, Paris, France, Maria.Tchernycheva@ief.u-psud.fr
– student, Ecole Polytechnique, Paris, France, nika.marek@mail.ru – Professor, PhD, Associate professor, ITMO University, Saint Petersburg,
Russia, alex_chirtsov@mail.ru
– postgraduate, ITMO University, Saint Petersburg, Russia,
shvagerdanila@gmail.com

Принято к печати 25.03.14 Accepted 25.03.14

148

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics
2014, №3 (91)