CC BY
37
НЛУЧНОЁ ИЗДАНИЕ МГТУ ИМ. Н.Э. ЬАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС 77-48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Математическое моделирование селективно излучающих разрядов с сильной неравновесностью
# 05, май 2012
Б01: 10.7463/0512.0364256
Градов В. М.
УДК 537.523+535.23
Россия, МГТУ им. Н.Э. Баумана [email protected]
Неравновесные селективно излучающие разряды, к которым относятся, например, разряды низкого давления в смеси паров ртути с инертным газом, являются объектом многочисленных фундаментальных и прикладных исследований. Библиография по физике данных разрядов весьма обширна и отражена, в частности, в монографиях [1-4]. Существующие теории основываются на ряде допущений, связанных с приближенным рассмотрением радиационных процессов (метод эффективного времени жизни с постоянным по радиусу коэффициентом поглощения), использованием усредненных по радиусу величин (концентраций и температур электронов и тяжелых частиц), упрощением схемы переходов между уровнями атомов, введением различных аппроксимаций для сечений упругих и неупругих столкновений частиц компонент и т.д. В настоящей работе строится вычислительная модель разряда, свободная от большинства указанных допущений и более реалистично учитывающая истинную картину процессов в плазме с сильной неравновесностью. Модель реализована в соответствующем программном обеспечении.
Модель разряда формулируется в гидродинамическом приближении [5]. В качестве различных сортов частиц рассматриваются электроны, ионы, атомы плазмообразующей компоненты в основном «0» и различных возбужденных «п»
состояниях с объемными концентрациями П0 и Пп , а также частицы буферного газа. Кинетическая часть модели строится на основе индивидуального учета группы низко расположенных уровней излучающих атомов с детальным рассмотрением столкновительных и радиационных механизмов их заселения и опустошения, и объединения большого числа уровней, примыкающих к континууму,
в квазинепрерывный спектр со свободными электронами. В этом спектре переходы описываются моделью диффузии электронов в пространстве энергии [6]. Скорости неупругих процессов возбуждения а п_к и ионизации а п атомов при соударениях с электронами и обратных им процессов в к _ п, в п рассчитываются с
помощью функции распределения электронов по энергиям, которая находится из уравнения Больцмана [5]. В нем в качестве параметров фигурируют температура
Те и концентрация Пе электронов, а также заселенности основного По и резонансного Пр уровней излучающих атомов. Перенос излучения рассматривается в
нескольких оптически плотных резонансных линиях с фойхтовским контуром, уширенных доплеровским, штарковским, резонансным и вандерваальсовским механизмами. Остальные линии учитываются в приближении объемного высвечивания.
Вычислительная модель включает уравнения энергии для электронов и тяжелых частиц (атомов и ионов) с температурами Те и Т, непрерывности и движения для заряженных частиц и возбужденных атомов и переноса излучения. В уравнении энергии для электронов учитываются неупругие потери на излучение и возбуждение частиц
А- А (Л - АА 4(5 кТе +1) ]+ав2 -
Я2х ах ах Ях ах 2
3 ш
- Пе Е -к (Те - т)2-*- V ек-Е ВпП п - ШуР* = 0 к 2 шк П
Уравнения непрерывности для возбужденных частиц, в которых учтено радиационное заселение из основного состояния, относятся к типу Бибермана- Хол-стейна
і а с
— -Т- 2.ІП = Е Пк(ПеРк-п + Акп) + П-Р П + 7- I к0- П^-
Ях ах к> П П! Д!
- Пп(а пПе + Е Пе«п-к + Е Пе Р п-к + Е А„к)
к> п к <п к <п
Система дополняется граничными условиями и рядом соотношений, отражающих протекание процессов в замкнутом объеме, и решается численными методами. В уравнениях обозначено: X = Г / Я- безразмерная радиальная координата в цилиндре радиуса Я, 1п - радиальный поток атомов, находящихся на
уровне п, % е, С - коэффициенты теплопроводности электронов и электропроводности плазмы, Еп,1- энергии возбуждения уровня п и ионизации атома, А пк - вероятность спонтанного распада уровня, ^ - интегральный по спектру радиальный поток излучения, к0_ п,и0_п - спектральные коэффициент поглощения и объемная плотность излучения в резонансной линии на частоте V, с-скорость света, V ек - частота упругих столкновений электронов с тяжелыми частицами.
Количество входных параметров модели ограничено: давления (масса) компонент, диаметр разрядной трубки, температура стенки, электрический ток. Применение описанной модели требует обширной базы данных по характеристикам элементарных процессов в плазме, включая многочисленные сечения упругих и неупругих столкновений, сечения фотоионизации, вероятности спонтанных радиационных переходов, параметры уширения линий. Эти данные частично заимствуются из различных источников или рассчитываются согласно [7, 8].
Ниже приведены некоторые результаты моделирования ртутных стационарных разрядов низкого давления с аргоновым буфером. В модели учитывались эффекты, связанные с изотопным составом природной ртути [2].
В области энергий электронов до порога возбуждения атома ртути на уровни
6р 3Р
0,1,2 (около 4.7 эВ) функция распределения электронов по энергиям мало отличается от максвелловской в силу доминирующего значения межэлектронных столкновений. В запороговом диапазоне возрастает роль неупругих процессов, и отличие функции распределения от максвелловской становится заметным (как в [2]).
Е, эВ
Рис. 1. Отношение скоростей возбуждения электронным ударом различных уров-
3
ней из состояния 6рР1 , рассчитанных при истинной и максвелловской функциях распределения электронов, в зависимости от энергии возбуждаемых
уровней.
Те = 12 • 103 К, п0 = 2 • 1014см- 3 ,заселенность резонансного уровня равна 0.3
от больцмановской.
11 11 12 3
1 -3 - пе = 10“,5 • 10“,2 • 10“см ' , соответственно. Маркеры на кривых соответствуют следующим уровням в порядке увеличения их энергий:
6р 3Р2, 7$381, 7 р 3Р(Ш, 6а 3О,,2 3, 8$ 381,8р 3Р(Ш,
9 $3 81, 9 р 3 Р0,1,2, 10 $3 81,10 р 3 Р0,1,2 '
3 3
Из рис.1 следует, что уровень 6р Р2 , ближайший к уровню 6р 3Р1 (пороговая энергия равна 0.57 эВ), возбуждается в основном максвелловскими электронами. Для уровня 7б 81 (порог 2.84 эВ) скорость возбуждения составляет
3
уже 0.7-0.92 от равновесной, а для уровня 10р ^0,1,2 (порог 5.14 эВ) - 0.4 - 0.8.
Видно, что степень влияния фактора неравновесности функции распределения на скорости возбуждения и ионизации зависит от порога процесса. Сопоставление данных теории и эксперимента проведено на рис. 2, 3. Нарастание температуры электронов к стенке связано с поглощением излучения, выходящего из центральных зон разряда с повышенной концентрацией излучающих частиц.
1.6Е+04
1.2Е+04
1_, 8.0Е+03 ш I-
4.0Е+03 0.0Е+00
0 0.2 0.4 0.6 0.8 1
г=г / К
Рис. 2. Распределение температур электронов (1,3) и тяжелых частиц (2) в разряде. 1,2-настоящая работа, R=19 мм, ток равен 0.42 Л, давление ^ и Лг 0.85 и 465 Па; 3-данные эксперимента[4], R=18 мм, ток равен 0.4 А, давление ^ и Лг 0.85 и
390 Па.
Перепад температуры тяжелых частиц между стенкой и центром не превышает 70К, температура тяжелых частиц во много раз отличается от температуры электронов (рис.2). Тонкой физической и важнейшей эксплуатационной характеристикой разрядов является поток излучения в резонансных линиях (рис.3). В рассмотренных условиях в линии 253.7 нм излучается от 50 до 70% подводимой электрической мощности. С ростом электрического тока мощность излучения нарастает линейно, а затем имеется тенденция к насыщению из-за возрастающей роли тушащих электронных ударов.
0.8
0 --------1--------------1--------------1-------------1-------
0.3 0.6 0.9 1.2 1.5
P, Вт/см
Рис.3. Зависимость потока излучения в линии 253.7 нм от электрической мощности. Давление Hg 0.85 Па. 1-данные эксперимента [4], R=18 мм, давление Ar 390 Па; 2- результаты настоящей работы, R=19 мм, давление Ar 465 Па.
Предложенная теоретическая модель удовлетворительно (в пределах 15 %) согласуются с результатами экспериментов по всем основным параметрам разрядов (сравнение проводилось также с данными других авторов, например, [2, 3, 9]) и открывает новые возможности в исследовании физических эффектов в плазме и разработке перспективных источников селективного излучения.
СПИСОК ЛИТЕРАТУРЫ
1. Грановский В.Л. Электрический ток в газе. Установившийся ток. - М.: Наука, 1971.- 543 с.
2. Миленин В.М., Тимофеев Н.А. Плазма газоразрядных источников света низкого давления. - Л.: Изд-во Ленинтрад. ун-та, 1991 - 239 с.
3. Рохлин Г.Н. Разрядные источники света. - 2-е изд., перераб. и доп. - М.: Энергоатомиздат, 1991.- 720 с.
4. Уэймаус Д. Газоразрядные лампы / Пер. с англ.; Под ред. Г.Н. Рохлина и М.И. Фугенфирова. - М.: Энергия, 1977. - 367 с.
5. Митчнер М., Кругер Ч. Частично ионизованные газы / Пер. с англ.; Под ред. А.А. Иванова. - М.: Мир, 1976. - 496 с. (Mitchner M., Kruger Ch. Partially Ionized Gases. Department of mechanical engineering Stanford university. New York London Sydney Toronto, 1973).
6. Биберман Л.М., Воробьев В.С., Якубов И.Т. Кинетика неравновесной низкотемпературной плазмы. - М.: Наука, 1982. - 375 с.
7. Sinyarev G.B., Gradov V.M. Complex Calculations of High Pressure Discharge Irradiation Characteristics // Beitrage aus der Plasmaphysik. - 1984. - Bd.24, №5. -p.499
- 528
8. Gradov V.M., Gavrish S.V. Mathematical modeling of selective emitting nonequilibrium plasma in complex optical systems // Light & Engineering.-1997. - v.5, №3.-p.16-19.
9. Kreher J., Stern W. Increased Power Concentration and its Effects on the Discharge Parameters of the Low Pressure Hg-rare Gas Positive Column. II. Variation of Tube Radius // Contrib. Plasma Phys. - 1989.- v.29, №2. - p.181-196.
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS 77 - 4821 1. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Mathematical simulation of selectively radiating charges with high non-equilibrium
# 05, May 2012
DOI: 10.7463/0512.0364256
Gradov V.M.
Russia, Bauman Moscow State Technical University
With the framework of multi-liquid hydrodynamic approximation, the author proposes a computational model of non-equilibrium space-inhomogeneous plasma of discharge in metallic vapor in the compound with inert gas which has different temperatures of light and heavy particles. All the functions of energy and state particle distribution were considered to be non-equilibrium and calculated in the course of evaluation. Beside collision processes in the kinetics of energy level population of plasma-generating components, radiative excitation by resonant lines expanded with different mechanisms is taken into account. Theory and experimental data are in good agreement (5-15%) with temperature fields and electron and heavy plasma particle concentration, and also with basic radiative and electro-physical characteristics of discharges in the mercury vapor with argon buffer.
Publications with keywords: computing model, non-equilibrium plasma, discharge Publications with words: computing model, non-equilibrium plasma, discharge
References
1. Granovskii V.L. Elektricheskii tok v gaze. Ustanovivshiisia tok [The electric current in the gas. Sustained current]. Moscow, Nauka, 1971. 543 p.
2. Milenin V.M., Timofeev N.A. Plazma gazorazriadnykh istochnikov sveta nizkogo davleniia [Plasma of gas discharge light sources of low pressure]. Leningrad, LU Publ., 1991. 239 p.
3. Rokhlin G.N. Razriadnye istochniki sveta [Discharge light sources]. Moscow, En-ergoatomizdat, 1991. 720 p.
4. Waymouth J. Electric Discharge Lamps. Cambridge, MA, The M.I.T. Press, 1971. (Russ. ed.: Ueimaus D. Gazorazriadnye lampy. Moscow, Energiia, 1977. 367 p.).
5. MitchnerM., KrugerCh. Partially Ionized Gases. Department of mechanical engineering Stanford university, New York- London -Sydney- Toronto, 1973. (Russ. ed.: Mitchner M., Kruger Ch. Chastichno ionizovannye gazy. Moscow, Mir, 1976. 496 p.).
6. Biberman L.M., Vorob'ev V.S., Iakubov I.T. Kinetika neravnovesnoi nizkotemper-aturnoiplazmy [Kinetics of nonequilibrium low-temperature plasma]. Moscow, Nauka, 1982. 375 p.
7. Sinyarev G.B., Gradov V.M. Complex Calculations of High Pressure Discharge Irradiation Characteristics. Beitrage aus der Plasmaphysik, 1984, vol. 24, no. 5, pp.499
- 528.
8. Gradov V.M., Gavrish S.V. Mathematical modeling of selective emitting nonequilibrium plasma in complex optical systems. Light & Engineering, 1997, vol. 5, no. 3, pp.16-19.
9. Kreher J., Stern W. Increased Power Concentration and its Effects on the Discharge Parameters of the Low Pressure Hg-rare Gas Positive Column. II. Variation of Tube Radius. Contributions to Plasma Physics, 1989, vol. 29, no. 2, pp.181-196.
