О зависимости параметров вольт-амперных характеристик коронного разряда от температуры

О зависимости параметров вольт-амперных характеристик коронного разряда от температуры

Ф. П. Гросуа, Ан. М. Бологаь, М. К. Бологаа, О. В. Моторина

aИнститут прикладной физики АН Молдовы, ул. Академическая, 5, г. Кишинев, MD-2028, Республика Молдова, e-mail:mbologa@phys.asm.md ьТехнологический институт Карлсруэ, Институт технической химии, Герман фон Гельмгольц Пл.1, Эггенштайн - Леопольдсхафен, 76344, Германия, e-mail: andrei.bologa@kit.edu

Анализируются зависимости параметров вольт-амперных характеристик коронного разряда от температуры газа. Показано, что напряжение зажигания коронного разряда Uc(T) для гелия и азота убывает с ростом температуры. Для гелия зависимость параметра А(Т) носит сложный характер, возрастая на начальном и конечном участках исследуемого интервала температур (20°С^369°С) и убывая в средней его части с двумя экстремумами. Ввиду нестабильности разряда в азоте определенной закономерности получить не удалось. Представлены обобщенные зависимости, позволяющие судить о наличии того или иного типа разряда, отличающегося от коронного. Выдвинуто предположение о существенной роли в наблюдаемых процессах двух эффектов: абсорбции ионов поверхностью электродов и непрерывного изменения подвижности ионов из-за изменения массы в связи с процессами кластеризации и декластеризации ионов, между которыми в стационарном состоянии предположительно существует динамическое равновесие.

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

УДК 537.527.3

ВВЕДЕНИЕ

Таунсендовский тип уравнения вольт-амперных характеристик (ВАХ) коронного разряда (КР) имеет вид [1-3]:

I = ЛП

где коэффициент пропорциональности

Л = а  е к,

в и к - относительная диэлектрическая проницаемость газа и подвижность ионов знака корони-рующего электрода; величина а зависит от геометрии электродов и их размеров; ис - напряжение зажигания коронного разряда. Для цилиндрической коаксиальной системы электродов (ЦСЭ) известно теоретическое значение а [2]:

8 пе„

a =

R2 • ln(R/го У

где г0, К - радиусы коронирующего и заземленного электродов.

Вольт-амперные характеристики коронного разряда зависят от параметров Л и ис, которые можно определить на основе статистической обработки экспериментальных ВАХ, например методом наименьших квадратов [4]:

A = I -U • Y; и = I •и -Y •U2

U2 - U

I - U • Y

I / U = Y.

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

туры за счет электрического ветра [5] или использованием КР в лазерной технике [6, 7], коронный разряд часто протекает в условиях широкого варьирования давления р и температуры Т газа. Поэтому важным представляется знание зависимости указанных выше параметров от этих термодинамических величин, то есть А(р,Т) и ис(р,Т). Решение задач предполагает нахождение параметров Л и ис статистическим путем. Первоначально устанавливаются вариационные ряды 1(и) при различных фиксированных давлениях и температурах и рассчитываются А(р,Т) и ис(р, Т) по формулам (4). Далее находят непосредственно сами зависимости по соответствующим уравнениям регрессии.

Ранее [8] была исследована зависимость параметров Л(р) и ис(р) от давления при заданных значениях температуры для гелия и азота. Установлено, что в гелии зависимость А(р) соответствует классическим представлениям об обратно пропорциональной зависимости коэффициента подвижности ионов от давления. Для напряжения зажигания короны ис подтвердилась формула Пика. В случае азота исследования показали наличие минимума в зависимости А(р), который можно предположительно конкретно объяснить на основе гипотезы о применимости уравнения Клаузиуса-Мосотти для диэлектрической проницаемости азота.

Цель данной работы - проведение аналогичных исследований применительно к определению соответствующих зависимостей А(Т) и

© Гросу Ф.П., Болога Ан.М., Болога М.К., Моторин О.В., Электронная обработка материалов, 2016, 52(1), 84-90.

ис(Т), к тем же газам и, по возможности, при аналогичных условиях.

ПАРАМЕТР ЩТ)

Теоретические предпосылки. Зависимость начального напряжения зажигания коронного разряда от давления и температуры качественно можно объяснить на основе формулы Пика, которая с точностью до числовых коэффициентов предположительно соблюдается не только для ЦСЭ, но и для ряда других электродных систем [1, 2]. Полагаем:

ис ~ Ес ~ 5-[ь + ^=г|| = Ь 5 + с->/5,

где Ь и с - числовые коэффициенты, зависящие от геометрии электродов.

Далее

5 Y . YRT о = —. p =

5 =

T0 • p

Yo M Po •T

здесь у и р - плотность и давление газа; R - универсальная газовая постоянная; T - абсолютная температура; M - молярная масса газа; индекс «0» указывает на состояние газа, близкое к нормальному (Т0 = 25°C; р0 = 1 атм) [2]. В результате с точностью до постоянных коэффициентов получим зависимость

Uc = b •

При p = const зависимость от температуры приобретает вид:

Uc = bjT + cj Jr,

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

Обсуждение экспериментальных данных. Применительно к гелию исходим из экспериментальной выборки ВАХ для различных температур при постоянном давлении в 20 бар, представленной на рис. 1. На рис. 2 приведено соответствующее обобщение согласно формулам [4]:

У3 = I/

Целый ряд точек выходит за пределы ошибок эксперимента (± 20%) относительно обобщающей зависимости. С целью более наглядного представления этих точек построены обобщенные графики (рис. 3) для каждой зависимости в отдельности, которые, включая общий (рис. 2), могут служить критериями принадлежности разряда к типу коронного. Исключая точки, не подчиняющиеся закономерностям коронного разря-

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

Так, при комнатной температуре (рис. 3а) коронный разряд проявляет себя в интервале безразмерных напряжений и*е(1;2), а при температуре 100°С (рис. 3б) - в интервале и*е(1,1;2,8). Для температуры 369°С (рис. 3з) можно прогнозировать интервал коронного разряда и*е(2; 3,5). Если, например, технологический процесс охватывает диапазон температур от 20°С до 369°С, то следует ориентироваться на обобщающую зависимость, представленную на рис. 2, из которой следует и*е(1,8;2,4). Близость точек к отрезку прямой является признаком наличия коронного разряда, подчиняющегося формуле (1). По обобщенной зависимости на рис. 2 можно предсказать, что в пределах ошибок эксперимента в интервале и*е(1,8;2,4) будет иметь коронный разряд. Там, где зависимость (1) не подтверждается, имеют место иные формы газового разряда. Графические зависимости на рис. 3 могут рассматриваться как «идентификаторы» коронного разряда. В подрисуночных данных рис. 3 указаны параметры А и ис, по которым можно найти их аппроксимации ис (Т) и А(Т).

На рис. 4 приведен график экспериментальной зависимости ис(Т) для обобщенных частных зависимостей от температуры, представленных на рис. 3, которая качественно соответствует уменьшению напряжения зажигания короны с ростом температуры. Однако в области Т = 500К наблюдается «аномалия» этой зависимости по сравнению с ожидаемой (8): появляется перегиб зависимости, не соответствующий формуле (8). В окрестности предполагаемой точки перегиба параметр ис почти не меняется с ростом температуры, после чего при температуре около 600К следует его спад (рис. 4). Можно предположить, что в основе данного физического явления лежит способность поглощения газов электродами, особенно таких легких газов, как водород и гелий [9]. Однако ограничимся констатацией экспериментальных фактов, когда наблюдаются уменьшение напряжения зажигания короны с ростом температуры и наличие плато в этой зависимости, после которого следует быстрый спад ис(Т) (рис. 4). Чтобы наглядно представить, насколько экспериментальные данные отклоня-

10 9 8 7

V 20 °С

+ 100 °с

•Сг 150 °С Д 200°С О 250 °С О 300 °с + 350 °С  369°С

2 4

и, кВ

Рис. 1. Экспериментальные ВАХ.

1,2 1,4 1,6 1,8 2 2,2 2,4 2,6

и/ис

(а) Г = 20°С,Л = 0,02931 мА/кВ2, ис = 7,8206 кВ. 5

4 3,5 3

К-

2.5 2 1.5 О,

V 20 "С + НЮ "С й 150 °С Л 200"С О 250"С О 300 °с X 350°С  369 "С

Л у » л

иг «6* &

1,5

2,5 У/11

3,5

Рис. 2. Обобщение данных рис. 1.

1,5 2 2,5

и/и,

(б) Т= 100°С,Л = 0,04943 мА/кВ2, ис = 5,6981 кВ.

(в) Г= 150°С, А = 0,03562 мА/кВ2, ис = 4,3050 кВ.

1 1,5 2 2,5 3 3,5 4 4,5 5

и/и,

(г)Г=200°С,Л = 0,04038 мА/кВ2, ис = 4,2411 кВ.

(д) г = 250°С, Л = 0,05142 мА/кВ2, ис = 4,0655 кВ.

(е) г = 300°С, Л = 0,07246 мА/кВ2, ис = 3.9796 кВ.

4,5 4

3,5 3

2,5 ">

1

+ 350 °С

2,5 3

и/и,.

3,5

4,5

4.5 4

3.5 3 2:5 2 1.5 1

369 °С

1

1,5

2,5

3,5

4,5

и/и,.

(ж) г = 350°С, А = 0,08511 мА/кВ2, ис = 3,0158 кВ.

(з) г = 369°С, А = 0,09990 мА/кВ2, ис = 2,7816 кВ.

с;

Рис. 3. Обобщающие графики для различных температур,

1«

Т, К

Рис. 4. Зависимость ис(Т). Не, р = 20 бар.

ются от зависимости типа Пика, на рис. 4 приведена кривая регрессии, аппроксимирующая опытные данные, согласно формуле:

2,094 -103

тт 2,094-103 2,613-10"

ис = —-+-

400

500

600

т, к

700

800

^ - .

" Т у/Т Т

Исследования, проведенные в азоте, показывают нестабильность протекания коронного разряда. На рис. 5 приведена зависимость, которая имеет преимущественно качественный характер. Наблюдается определенное сходство с коронным разрядом в гелии (плато на рис. 4), однако в азоте вблизи точки 470К наблюдается не плато, а локальный минимум. Что же касается воздуха, то, согласно данным [10], наблюдаемые зависимости хорошо согласуются с формулой Пика.

На рис. 6 приведены вольт-амперные характеристики, снятые при КР в коаксиальной модели электродов [10], а на рис. 7 - обобщения по предлагаемой методике аналогично рис. 1 и 2. «Критерий» (рис. 7) принадлежности к разряду типа КР соблюдается с высокой точностью. Следует ожидать и соблюдения формулы Пика (8), о чем свидетельствует график (рис. 8), построенный согласно уравнению регрессии, аналогичному (8):

Рис.5. Зависимость ис(Т). Азот р = 7 бар.

„ 1,175 -103 73,86

ис = --+—^.

с Т ТТ

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

Неотъемлемой частью КР является электрический ветер (ЭВ), как следствие наличия в разряде свободных нескомпенсированных зарядов объемной плотностью р, которые, будучи обусловлены [11] множителем (и-ис) ~ р в (1), отражены в вольт-амперных характеристиках (в виде биссектрис). Следовательно, эти характеристики содержат информацию и об электрическом ветре.

ПАРАМЕТР А

Теоретические предпосылки. Считаем диэлектрическую проницаемость гелия постоянной. Тогда, согласно формуле (2), следует ожидать, что зависимость параметра А от температуры можно объяснить зависимостью подвижности активных ионов от температуры к = к(Т). Обратимся к классическим зависимостям подвижности газовых ионов от температуры (давление газа считаем постоянным). Так, согласно теории Ланжевена, подвижность ионов должна быть обратно пропорциональной среднеквадратичной

2,8 2,6 2,4 2.2 2 1,8 1,6 1,4 1,2 1

АЗШ к к

К.

10 12

и, кВ

Рис. 6. Исходные ВАХ для воздуха; «+» корона.

1 1.2 1,4 1,6 1,8 2 2.2 2,4 2,6 2,8

и/ис

Рис. 7. Обобщающая ВАХ.

сс

а

300

400

500 Т, К

600

700

Рис. 8. Зависимость ис(Т). Воздухр =1 бар.

а

аз

< 0,06

< 0,06

1,2 1,4 1,6 1,8

Рис. 9. Аппроксимация экспериментальных данных А(Т) полиномом третьего порядка. Точки - эксперимент.

скорости и и прямо пропорциональной длине свободного пробега иона Я, которая, в свою очередь, обратно пропорциональна плотности газа. Следовательно, подвижность ионов должна быть пропорциональна корню квадратному из температуры, в то время как, согласно другим исследованиям, в частности теоретическим, температурная зависимость подвижности обратно пропорциональна корню квадратному из температуры [12]. Учет зависимости эффективного поперечного сечения столкновений частиц от температуры, по формуле Сюзерленда, также не привел к достаточной ясности [13] в рассматриваемом вопросе. Согласно работам [14, 15],

Рис. 10. Аппроксимация данных эксперимента А(Т) полиномом пятого порядка (11). Точки - эксперимент.

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

т-1/2

ры, как Т в случае малых ионов, например для гелия, что качественно согласуется со справочными и экспериментальными данными [16], но возрастает с повышением температуры для крупных ионов кластерной структуры. В этих случаях проявляются две противоположные тенденции: кластеризации и декластеризации ионов, что может приводить к росту подвижности с температурой, убыванию или постоянству подвижности [14].

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

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

А

где температура нормирована на Т0 = 293К (20°С), то есть Т* = Т/Т0. Что касается зависимости А(Т) для азота, то в данном случае возникли определенные проблемы с воспроизводимостью опытных данных, что свидетельствует о необходимости проведения дополнительных экспериментальных исследований. Применительно к обобщениям на рис. 6-7 для воздуха, по данным работы [10], укажем на их хорошую воспроизводимость. Соответствующий график А(Т) приведен на рис. 11 (положительная полярность коронного разряда).

х КГ3

1

0,9

0,8

0,7

tN сс 0,6

1 0,5

S

< 0,4

0,3

0,2

0,1

0

300 350 400 450 500 550 600 650 т, К

Рис. 11. Экспериментальная зависимость А(Т) для воздуха при атмосферном давлении.

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

ВЫВОДЫ

1. В работе приведены общие и раздельные обобщения ВАХ КР при различных температурах газа, позволяющие «идентифицировать» коронный разряд.

2. Зависимость Uc(T) для гелия и азота носит сложно убывающий характер. Для воздуха соблюдается монотонно убывающий характер зависимости в соответствии с формулой Пика.

3. В зависимости А(Т) для гелия в интервале температур от 20 до 369°С наблюдается возрастающая ветвь до максимума, затем спад до минимума с дальнейшим подъемом. Эта зависимость хорошо аппроксимируется полиномом пятой степени. Для азота такая зависимость не установлена. В случае воздуха наблюдается слабая зависимость рассматриваемого параметра от температуры.

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

Работа выполнена при финансовой поддержке по двустороннему проекту A§M-BMBF: А§М: 13.823.15.09/GA, Институциональному проекту: 15.817.02.07А и BMBF:FKZ-Nr. 01DK13014.

ЛИТЕРАТУРА

1. Верещагин И.П., Левитов В.И., Мирзабекян Г.З., Пашин М. Основы электродинамики дисперсных систем. М.: Энергия, 1974. 420 с.

2. Райзер Ю.П. Физика газового разряда. М.: Наука,

1992. 536 с.

3. Самусенко А.В., Стишков Ю.К. Электрофизические процессы в газах при воздейстии сильных электрических полей. Санкт-Петербург: Изд-во Госуниверситета, 2011. 567 с.

4. Grosu F.P., Bologa An.M., Bologa M.K., Motorin O.V.

Similarity Methods and Generalizations in Study of Corona Discharge. Proceedings of 11th International Conference on Modern Problems of Electrophysics and Electrohydrodynamics (MPEE-2015). June 29 -July 3, 2015, Peterhof, Russia, pp. 129-135.

5. Болога М.К., Гросу Ф.П., Кожухарь И.А.

Электроконвекция и теплообмен. Кишинев: Штиинца, 1977. 320 с.

6. Salvermoser M., Murnick D.E. Efficient, Stable, Corona Discharge 172 nm Xenon Excimer Light Source. JAppl Phys. 2003, 94(6), 3722.

7. Lo Shui-Yin, Lobo Julio D., Blumberg Seth, Dibble Theodore S., Zhang Hu, Tsao Chun-Cheng, Okumura Mitchio. Generation of Energetic He Atom Beams by a Pulsed Positive Corona Discharge. J Appl Phys. 1997, 81(9), 5896-5905.

8. Гросу Ф.П., Болога Ан.М., Болога М.К., Моторин О.В. Влияние давления газа на вольт-амперную характеристику коронного разряда. Сборник докладов XI конференции «Волновая электрогидродинамика проводящей жидкости. Долго-

живущие плазменные образования и малоизученные формы естественных электрических разрядов в атмосфере». 04.07. - 07.07. 2015 года, Ярославль, Россия, с. 56-63.

9. Черданцев Ю.П., Чернов И.П., Тюрин Ю.И. Мето-

ды исследования систем металл-водород. Учебное пособие. Томск: Из-во ТПУ, 2008. 286 с.

10. Токарев А.В. Коронный разряд и его применение. Бишкек: КРСУ, 2009. 138 с.

11. Гросу Ф.П., Болога Ан.М., Болога М.К., Моторин О.В. О моделировании коронного разряда методами теории подобия. ЭОМ. 2014, 50(2), 41-48.

12. Смирнов Б.М. Диффузия и подвижность ионов в газе. УФН. 1967, 92(1), 75-103.

13. Леб Л. Основные процессы электрических разрядов в газах. Москва-Ленинград: Госиздат технико-теоретической литературы, 1950. 672 с.

14. Karpas Z., Eiceman G.A., Krylov E.V., Krylova N. Models on Ion Heating and Mobility in Linear Field Drift Tubes and in Differential Mobility Spectrometers. Int J Ion Mobility Spectrom. 2004, 7, 8-18.

15. Zeev Karpas, Zvi Berant, Oded Shahal. Effect of Temperature on the Mobility of Ions. J Am Chem Soc. 1989, 111(16), 6015-6018.

16. Таблицы физических величин. Справочник. Под

ред. акад. И.К. Кикоина. М.: Атомиздат, 1976.

1008 с.

Поступила 20.05.15

Summary

The dependences of parameters of current-voltage characteristics of corona discharge on temperature are analyzed. It is shown that the inception voltage Uc(T) for both helium and nitrogen decreases as the temperature increases. The dependence of the parameter A (T) for helium is of a complex character, it increases both in the initial and final parts of the investigated range of temperatures (20°C - 369°C) and decreases in its middle part with two extrema. It was not possible to receive such dependence for nitrogen due to the corona discharge instability. Presented generalized dependences give possibility to distinguish the corona discharge from other kinds of discharge. It is assumed that two factors play a significant role in the investigated process: ion absorption by the surface of electrodes and continuous change of ion mobility due to the ion mass change in the processes of ions clasterization and declasterization. In a steady slate there is a dynamic equilibrium between those processes.

Keywords: corona discharge, current-voltage characteristics, current, voltage, ion mobility, temperature, pressure, cluster.