CC BY
107
Водород выделяется при облучении из металлов преимущественно из точки воздействия электронного пучка.
"Сильный" уровень неравновесности в водородной (дейтериевой) подсистеме металлов при облучении возможен благодаря квазиизолированности ее колебательного спектра от фононного спектра металлов и относительно высокой энергии неравновесных колебаний, например, энергия плазменных колебаний в протонной подсистеме металла достигает « 0,2 эВ. При таких условиях возбуждения обнаружена эмиссия положительных ионов (вероятно Н+, Б+) из Рё (Б,Н), облучаемого электронами.
"Красная граница" эффекта находится в облас-
ти энергии электронов ~6 кэВ, что соответствует передаче при электронном ударе протонам энергии 13...15 эВ.
Ток эмиссии положительных ионов (10-3...10-2мкА) растет квадратично с током первичного электронного пучка. Это связано с двухстадийностыо процесса возбуждения эмиссии - стимуляцией излучением неравновесного выхода водорода из объема на поверхность Рё и "выбиванием" протонов (дейтронов) с поверхности при электронном ударе адсорбированных атомов водорода (дейтерия).
Наблюдаемые эффекты раскрывают новые каналы релаксации энергии первичного излучения в металлах насыщенных водородом.
СПИСОК ЛИТЕРАТУРЫ
1. Тажибаева И.Л., Шестаков В.И., Руденко И.В., Чуж-ко Р.К. Изучение возможности стимулирования процесса диффузии водорода в металлах // Журн. физ. хим. - 1984. - Т. 58. - № 5. - С. 1208-1211.
2. Dobrozemsky R., Schwarzinger G., Stratova C. Effect of Radiation on Hydrogen Penetration Stainless Steel // Sup. Atom and Surface Phys. - Salzburg: Contr. Junsb. 1980. - P. 64-69.
3. Polushkin B.G. The effect of Radiation on the Diffusion of Hydrogen in Metals // J. Nucl. Mater. - 1992. - Vol. 94. - P. 219-223.
4. Тюрин Ю.И., Чернов И.П. Радиационно-стимули-рованный выход водорода из металлов и сплавов // Физ. ХОМ. - 1998. - № 2. - С. 9-13.
5. Чернов И.П., Черданцев Ю.П., Тюрин Ю.И. Миграция водорода в стали и сплавах, стимулированная излучением // Изв. вузов. Физика. - 1994. - № 11. - С. 72-79.
6. Chernov I.P., Tyurin Yu.I. Hydrogen migration in stainless steel and titanium alloys, stimulated by ionizing radiation // J. Nucl. Mater. - 1996. - Vol. 233-237. - P. 1118-1122.
7. Chernov I.P., Tyurin Yu.I., Cherdantzev Ju.P. Hydrogen migration and release in metals and alloys at heating and radiation effects // J. Hydrogen Energy. - 1999. - Vol. 24.
- P. 359-362.
8. Тюрин Ю.И., Чернов И.П., Смекалина Т.В. Механизмы радиационно-стимулированной диффузии и выхода водорода из металлов // Физ. ХОМ. - 1998. -№ 3. - С. 11-16.
9. А. с. 1695181 СССР. Способы определения атомарного водорода // Л.В. Сериков, Л.Н. Шиян, А.Ф. Горбачев, Ю.И. Тюрин. 1988. - Бюл. изобр. - 1991. - № 44.
- С. 44.
10. Зандберг Э.Я., Ионов Н.И. Поверхностная ионизация. - М.: Наука, 1969.
11. Тюрин Ю.И., Чернов И.П. Неравновесный выход атомарного водорода из металлов при облучении // Доклады АН СССР. - 1999. - Т. 367. - № 3. - С. 328332.
12. Тюрин Ю.И., Чернов И.П., Афанасьев Б.В., Смекалина Т.В. Эмиссия положительных ионов при облучении металлов электронами // РФХ-10. Материалы Международной конф. - Томск, ТПУ. 1999. - С. 318319.
УДК 533.9
ЭЛЕКТРОДИНАМИКА ВЫСОКОЧАСТОТНОГО ФАКЕЛЬНОГО РАЗРЯДА
И.А. Тихомиров, В.А. Власов, Ю.Ю. Луценко, А.А. Зорин
Проведен анализ электродинамики высокочастотного (ВЧ) факельного разряда. Предложена и обоснована электродинамическая модель ВЧ факельного разряда.
Введение
ВЧ разряды, в соответствии с установившейся классификацией [1-4] называются ВЧ-индукцион-ными (ВЧИ), ВЧ-емкостными (ВЧЕ) и ВЧ-дуговы-ми (ВЧД) разрядами. ВЧ-индукционный, так называемый безэлектродный Н-разряд, был описан еще в 1884 году [5]. В ВЧ-дуговых разрядах между
электродами внутри разрядной камеры возбуждается дуговой Е-разряд. Обычный ВЧЕ разряд возбуждается в разрядной камере при наложении электродов на внешние диэлектрические стенки камеры. Среди ВЧ разрядов особое место занимает, так называемый, ВЧ факельный разряд. Разряд был открыт и так назван С.И. Зилитинкевичем в 1928 г. [6], т.к. разряд (рис. 1) с электрода - 1 горит в виде
факела, в котором высокотемпературный канал - 2 окружен светящейся диффузной оболочкой - 3, диаметр которой значительно больше, чем у канала. Энергия к разряду подводится от ВЧ генератора через электрод - 1.
По своей природе и свойствам ВЧ факельный разряд является одноэлектродным Е-разрядом, который горит с учетом емкостной связи факел-земля [7]. Вид, форма, режимы горения факельного разряда, его специфика описаны в работах [1-3, 819]. Разряд возбуждается в любой (включая и водород) газовой среде при давлениях 10...5-105 Па. Исследования показали, что плазма ВЧ факельного разряда обладает ярко выраженной неравновесностью. В зависимости от условий электронная температура в плазме разряда изменяется от 5000 до 35000 К и выше, а газовая (молекулярная) - от 1500 до 5000.9000 К. В каждом конкретном случае электронная температура в несколько раз (иногда на порядок), выше, чем газовая. Степень ионизации плазмы может изменяться с учетом электронной плотности 1015...1022 м-3 в диапазоне значений 10-8 ...10"1.
В свое время [18, 20] была предложена приближенная электродинамическая модель разряда. Канал разряда рассматривается в этой модели как однородный проводник с распространяющейся по нему основной ТМ-волной. Выделяемая в разряде энергия определяется омическим сопротивлением, а ее рассеяние - радиальной теплопроводностью. Моделирование для более мощных ВЧ факельных разрядов, когда их геометрические размеры становятся соизмеримыми с длиной двухпроводной линии анодного контура ВЧ генератора, рассмотрено в работах [1, 21, 22]. С учетом слабого скин-эффекта диссипация электромагнитных волн различного типа в разряде, а не только основной ТМ волны проанализирована в [17]. В плазме канала, со слабо изменяющейся в осевом направлении величиной
комплексной диэлектрической проницаемости, могут возбуждаться симметричные поперечно-электрические (ТЕ) и поперечно-магнитные (ТМ) волны, а также несимметричные волны, являющиеся гибридными (рис. 2). Характер электромагнитного поля разряда изучался (рис. 1) с помощью од-новитковой рамки - 4, присоединенной к осциллографу. Численные расчеты и эксперименты [17] показали, что при средних мощностях определяющий вклад в энергобаланс разряда вносит основная ТМ волна. Вклад побочных ТМ волн незначителен и слабо зависит от мощности и частоты. Следует учитывать, что вклад ТЕ волн при средних мощностях (до 10 кВт) меньше, чем от основной ТМ волны, но с ростом мощности разряда и частоты доля вклада от ТЕ волн возрастает [17].
В работе [23] изучено распределение (рис. 1 и рис. 2) компонент электромагнитного поля Н^ и Ег и их затухание в канале разряда; определено радиальное распределение удельной электропроводности плазмы в разряде. Перегиб профиля удельной электропроводности наблюдается на границе раздела "канал-диффузная оболочка". При этом максимум магнитного поля Н^ наблюдается на перегибе, а максимум Ег - в пределах диффузной оболочки разряда - 3.
Данные о свойствах и особенностях ВЧ факельного разряда, специфике его электромагнитного поля, логично позволяют нам перейти к анализу электродинамики разряда с обоснованием новой электродинамической модели разряда в отличии от модели, предложенной в работах [18, 20].
1. Распределение электрических токов
в канале высокочастотного факельного разряда
В отличие от работ [18, 20] будем рассматривать канал разряда в виде излучателя конечной длины. В этом случае канал разряда можно предс-
Иъ
/
И
2 4
Рис. 1. ВЧ факельный разряд
ТМ
Еъ
Е,
-0
и
ТЕ
Иъ
Ет
И,
г
Рис. 2. Распределение компонент электромагнитного поля в канале разряда
ъ
г
0
тавить [25] в виде тонкого линейного излучателя конечной длины. Распределение токов в таком излучателе в соответствии с интегралом Кирхгофа-Гюйгенса [25] будет зависеть от Н^.
НР = 4п
1 2 -Ч
4п J
jkyl (z') е
-jkr'
2
yI (z') e
-jkr'
dz',
(1)
Здесь. /(г') - распределение электрических токов вдоль линейного излучателя; Г - расстояние от элемента излучателя с координатой г' до точки наблюдения Р, имеющей координаты х, у, г; к -постоянная распространения электромагнитной волны в непроводящей среде; Н„- аксиальная компонента магнитного поля (рис. 1).
Электрическую составляющую электромагнитного поля можно легко определить через известную величину магнитной составляющей посредством следующего выражения [25].
Е = (1/>£) Ух И.
Прежде, чем проводить анализ экспериментальных результатов, рассмотрим распределение компонент поля для случая "электродинамической модели" разряда по работе [18]. Согласно этой "электродинамической модели" вдоль канала ВЧФР распространяется затухающая поперечно-магнитная волна. При этом предполагается, что отражение ее на конце канала разряда несущественно, и им можно пренебречь. В этом случае распределение электрических токов в канале ВЧФР запишется в виде [18].
I = /о^; у = а + ] в, (2) где у - комплексный коэффициент распространения электромагнитной волны (волновое число) вдоль канала разряда.
Преобразуем интеграл (1) с учетом выражения
(2).
Получим.
1гу10и ку10¥'
4%нр =
,'2
dz' +
+ j
-i
yI oU + kyI V
2
dz'
(3)
где:
U = e^'cosKßz' + кГ)]; V = е~а' sin[-(ßz' + кГ)]. Результаты расчетов величин интеграла (3) и, соответственно, распределения аксиальной компоненты магнитного поля вдоль оси ВЧФР для раз-
„I I /2 , ß2 ß
личных значений \у\ = ^а +ß ; arg у= —
а
при y =310-2 м приведены на рис. 3. Расчет прово-
дился методом Симпсона.
Как видно из рисунка, постоянство аксиальной компоненты магнитного поля для случая "электродинамической модели" [18] может наблюдаться только на расстояниях от электрода, не превышающих нескольких сантиметров, что составляет величину не более 15...20 % от длины канала ВЧФР. Заметим также, что и характер изменения фазового сдвига вдоль оси ВЧФР, рассчитанный на основе "электродинамической модели" не соответствует полученной нами экспериментальной кривой.
Таким образом, полученные нами экспериментальные результаты не могут быть интерпретированы соответствующим образом в рамках "электродинамической модели" [18].
Перейдем теперь непосредственно к задаче определения распределения электрических токов в ВЧФР по результатам измерений аксиальной компоненты магнитного поля в ближней зоне его излучения.
В общем случае задача нахождения распределения электрических токов в излучателе по данным измерений создаваемого им поля достаточно сложна, и ее решение математически разработано [25] только для анализа компонент поля в дальней зоне излучения. Однако ее приближенное решение несложно получить, основываясь на том или ином численном методе расчета интеграла Кирхгофа-Гюйгенса.
В соответствии с методом трапеций приближенное значение интеграла (2) запишется в виде [26].
H,
Р
(z)=±]
4п i
jkyl(z')e jkr'
yI(z')e jkr'
dz' =
:J F (z '; r'; y )dz' >
1
-[ F ( z') + F ( z'n)] + +F (z') +... + F ( z'n_i)
ly=const
Подынтегральную функцию F(r ;z,';y') можно представить в виде: F(r ;z';y') =a(z;Z )|y=constI(z'),
где: a(z,z r) = -i-4n
jk 1
2
- jkr .
Выбрав соответствующим образом п экспериментально полученных значений аксиальной компоненты магнитного поля мы можем записать п линейных уравнений на п неизвестных величин. В матричном виде данная система уравнений запишется следующим образом.
[/ • [А] = [В], (4)
где [I] - вектор-столбец неизвестных электрических токов, соответствующих определенным участ-
1
кам канала разряда; [А] - матрица коэффициентов а(г;г'); [В] - вектор-столбец экспериментально полученных значений аксиальной компоненты магнитного поля.
В качестве примера рассмотрим решение задачи о нахождении распределения электрических токов излучателя по известным значениям аксиальной компоненты магнитного поля для случая ВЧФР с длиной канала Х=18 см, что соответствует мощности, вкладываемой в разряд Ж «2 кВт.
Разобьем канал разряда на 6 отрезков длиной 3 см. В этом случае матрица коэффициентов а(1;1') при у =30 мм будет иметь следующий вид:
А =
18524 13103 3318 1175 531 282 83
6551 37049 13103 3318 1175 531 140
1659 13103 37049 13103 3318 1175 265
587 3318 13103 37049 13103 3318 587
265 1175 3318 13103 37049 13103 1659
141 531 1175 3318 13103 37049 6551
83 282 531 1175 3318 13103 18524
Здесь с целью упрощения мы пренебрегаем запаздыванием сигнала, что выражается в отбрасывании мнимых частей коэффициентов, величина которых для ближней зоны излучения сравнительно мала.
Вектор-столбец значений аксиальной компоненты магнитного поля составим, пользуясь имеющимися экспериментальными данными. В данном случае он будет иметь следующий вид:
В =
1,00 + 7 0,00 0,99 + 7 0,09 0,90 + 7 0,19 0,96 + у 0,26
В =
0,94 + 7 0,32 0,82 + 7 0,57 0,69 + 7 0,72
Здесь с целью упрощения анализа мы используем относительные, а не абсолютные значения аксиальной компоненты магнитного поля. Переход к абсолютным значениям несложно осуществить умножением полученных значений амплитуд токов на соответствующий множитель.
Решение системы уравнений (4) найдем путем использования метода Гаусса. После преобразований получим значения величин электрических токов, соответствующих определенным участкам канала ВЧФР. Вектор-столбец электрических токов будет иметь вид:
I =
3,88-10-5 • 73,000 е
2,16 40 -5 7 3,80 е
1,21 •Ю -5 е710,10
2,24 •Ю -5 714,00 ее
I =
1,38 • 10~ 1,11 • 10
7 5,10
е
/ 55,10
-5 7 55,00 4,50 • 10 5 • е-1
"Выбросы" амплитуды электрического тока в крайних точках связаны с краевыми эффектами, а также с приближенностью метода решения задачи. Однако в реальности краевые эффекты не столь существенны для случая факельного разряда, так как с одной стороны канал разряда электрически связан и является своего рода продолжением электрода, с другой стороны канал разряда не имеет резкой границы. Расчет, проведенный с учетом диффузной части разряда, подтверждает вышеприведенные соображения.
На рис. 4 приведены распределения электрических токов вдоль оси свободного ВЧФР с учетом и без учета диффузной части разряда, а также для случая экранированного ВЧФР.
Как видно из рис. 4, постоянству амплитуды аксиальной компоненты магнитного поля в ближней зоне излучения разряда соответствует распределение электрического тока с неэкспоненциальным характером изменения амплитуды вдоль оси разряда.
Заметим, что отсутствие существенных изменений электрических параметров плазмы разряда вдоль его оси можно объяснить наличием:
- слабого затухания электромагнитной волны, за-питывающей разряд;
- отраженной электромагнитной волны в конце разряда.
При отсутствии затухания электромагнитной волны горение разряда невозможно, так как в этом случае не происходит передачи энергии электромагнитной волны к плазме разряда.
Такими образом, с целью интерпретации полученных расчетных и экспериментальных результатов мы должны предположить наличие отраженной электромагнитной волны в канале факельного разряда [27-29].
Дф,
рад
1,2 - - 2,0
1,0 1 —уСХ -1 1,6
0,8 \ / 0,8
0,6 0
Гг^ \/
0,4 \ - 0,8
0,2 2 1,6
1 1 1 -----
9 12 Ъ 102, м
15
18
Рис. 3. Расчетное распределение аксиальной компоненты магнитного поля вдоль оси ВЧФР: 1) относительная амплитуда; 2) фазовый сдвиг. (1, 2) - ¡у! = 45 м-1; агду=л/4 (11, 2) - ¡у! = 30 м1; агду=л/4
5
0
3
6
Дф,
град
40 20 0 20 40
0,25 0,50 0,75 Z = ъ/Ь
б)
1,00
Дф,
град
40 20 0 -20 -40
0,25 0,50 0,75 Z = ъ/Ь
в)
1,00
Дф,
град 40
20 0 20 ■40
Рис. 4.
0,25 0,50 0,75 Z = ъ/Ь
Распределение амплитуды (1) и фазового сдвига (2) высокочастотного тока в канале ВЧФР: без учета (а) и с учетом диффузной части разряда (б); (в) -экранированный разряд
2. Электродинамическая модель высокочастотного факельного разряда с учетом отраженной электромагнитной волны
Впервые высокочастотный факельный разряд был рассмотрен в виде линии с распределенными параметрами конечной длины в работе [30]. Авторы вышеуказанной работы справедливо указали на недостатки электродинамической модели ВЧФР [18], основанной на представлении разряда в виде цилиндра бесконечной длины, которые выражались, прежде всего, в расхождении теоретических расчетов с экспериментальными результатами. Заметим, что модель бесконечного цилиндра применима только в случае отсутствия в канале ВЧ факельного разряда отраженной электромагнитной волны. Другими словами, учет конечности длины канала ВЧФР при теоретическом рассмотрении выражается прежде всего в учете наличия в канале разряда отраженной электромагнитной волны, что и сделано в предлагаемой модели.
Рассмотрим волну электрического тока, распространяющуюся вдоль линии конечной длины (см. рис. 5). В общем виде ток, протекающий в точке г, определяется [26] следующим образом:
=
1е
ух
-у2
+ се ^ е '
(5)
Здесь: I, 10 - комплексные амплитуды ВЧ-то-ков, соответственно в точках г и г = 0; у = а + ]в -комплексный коэффициент распространения (волновое число) электромагнитной волны в канале ВЧФР; се К - комплексный коэффициент отражения электромагнитной волны на конце линии.
Выясним далее, при каких значениях параметров а, в, С,, с, рассчитанное на основе экспериментальных результатов распределение электрического тока в канале ВЧФР может быть удовлетворительно описано выражением (5).
Результаты расчетов распределения амплитуды и фазового сдвига электрического тока вдоль канала ВЧФР, проведенные согласно выражению (5) для различных значений а, в представлены на рис. 6. Сопоставляя кривые распределений электрического тока, изображенного на рис. 4, можно видеть, что наибольшее совпадение распределений электрического тока, полученных, соответственно, на основе экспериментальных данных и в результате расчетов по формуле (5) наблюдается в случае а = 3 м-1, в< 30 м-1, | = 0, с = 1.
Ъ
О-
I
Ъ=0
Рис. 5. Схема линии конечной длины
I
I
0
Дф,
1/10 1,0 град 1/10 1,0
\ Г \ л 2 40 20
0,5 1 0 -20 0,5
0 | | \У_ | -40 0
0,25 0,50 0,75 1,00
1/1о
1,5 1,0 0,5 0
1/1о 1,0
0,5 0
1/10
ча=1 м"1 р=40 м"
X
X
0,25
0,50 0,75
а=5 м"1 в=40 м
-1
X
X
0,25
0,75
0,50
а=1м-1 р=30м-1
0,25 0,50 0,75 а=5 м"1 в= 30 м-1
Дф,
град
60 40 20 0 -20 -40 -60 -80 1,00
Дф,
град
40 20 0 -20 -40
1,0 -
2
0,5 \ V// _
1
0 1 1 1
Дф,
град 1/10
45
1,0
30
15
0
-15 0,5
-30
-45
-60 0
1,00
1/10 1,5 1,0 0,5 0
1/10 1,0
0,5 0
_1_
0,25 0,50
а=3 м-1 в=40 м-1
_1_
X
0,25 0,50 0,75 а=6 м-1 в=40 м-1
X
0,25
0,50
0,75
а=3м-1 р=30м-1
X
X
Дф,
град
40 20 0 -20 -40
0,75 1,00
Дф,
град
0,25 0,50 0,75 а=6 м-1 в=30 м-1
60 40 20 0 -20 -40 -60 -80
1,00
Дф,
град
40 20 0 20 40
Дф,
град 45 30 15 0
-15
-30 -45
-60
1,00
Рис. 6. Расчетные кривые распределения амплитуды (1) и фазового сдвига (2) электрического тока в канале ВЧФР
1
Малую величину коэффициента затухания электромагнитной волны, распространяющейся вдоль канала ВЧФР, подтверждают также измерения фазового сдвига между радиальной компонентой магнитного поля Er и аксиальной компонентой магнитного поля Иф. Заметим, что для цилиндрических структур в непроводящей зоне отношение комплексных амплитуд Er и Иф непосредственно определяется [31] величиной волнового числа h:
Er Ицсо
— =-^г- (6)
H„ kl
Учитывая, что величина - действительная,
kl
из выражения (6) получим:
J
Er H,
P
= Ce^,
(7)
0,3
0,2
W, кВт
Рис. 7. Зависимость соотношения в/а от мощности разряда W
1=еЕ +еН ^г пр •
Учитывая, что геометрические размеры рамки и частота электромагнитного поля постоянны, можно записать:
JYr JVE jY
e T = Ae r + Be r sin p ,
H,
P
где А и В - постоянные действительные величины, характеризующие величины сигналов, поступающих на вход измерительного прибора соответственно от электрической и магнитной компонент поля.
Преобразуя выражение (7), получим:
где С - действительная константа; в - коэффициент затухания; а - фазовый коэффициент (к = а + _/в).
Таким образом, согласно выражению (7) по величине фазового сдвига между радиальной компонентой электрического поля и аксиальной компонентой магнитного поля однозначно определяется отношение коэффициента затухания к коэффициенту фазы электромагнитной волны, распространяющейся в плазме ВЧФР.
Эксперимент по определению характера магнитного поля ВЧФР позволяет также провести оценку разности фаз между компонентами поля Ег
и И-
Рассмотрим более подробно смещение кривой сигнала I в зависимости от угла вращения рамки р (рис. 1), обусловленное наводками компоненты поля Ег. Сигнал от рамки, поступающий на вход измерительного прибора определяется как сумма ЭДС, наводимых магнитным и электрическим полями
р/а
0,4
A2 + B2 sin2 p + +2 AB cos
ÍYe
-Y
H.
P
sin ф
(8)
^sínY E + BsinT H sinф Er Нф
T J = arctg---. (9)
J ^cosT E + BcosT H sinф
Er Нф
Пусть для определенности Y E = 0 .
r^r
Тогда выражение (9) примет следующий вид: B sin Tн slnф
T j = arctg -
P
A + B cos Y тт sin p
HP
При p= 0 очевидно, что Y7 = 0. Получим:
ctg Y j =
С
sin Y tj sin p Hp
(10)
Здесь С, = —.
в п п
Выражение (10) в случае р = — и р =--
2 2
запишется в виде системы уравнений:
ctg Y T
С
Р=П "sin Y H
+ ctg Y
H
(11)
P
P
ctg Y,
С
= n = ctgY H - Y
— np sin Yh
P=
•• (12)
P
Для случая ВЧФР, горящего при атмосферном
давлении мощностью W=1 кВт на расстоянии от электрода l =60 мм, имеем:
р=П
= 7
р=-П
= -40
Соответственно, в результате решения систему уравнений (11-12) получим, ^е н = 16;
С = 1,22. в
Используя (7), будем иметь — = 0,287 .
а
в
Результаты измерений — в зависимости от а
мощности разряда приведены на рис. 7.
Как видно из рис. 7, величина в находится в
а
пределах 0,25...0,35.
Учитывая, что длина электромагнитной волны имеет величину большую, чем длина канала, т.е. 2п/а > Ь, получим, что падение интенсивности поля в е раз происходит на расстоянии от электрода, которое сопоставимо или превышает длину канала разряда, так как в ^ 1,8/Ь. Таким образом, полученная нами оценка величины коэффициента затухания для ВЧФР мощностью 0,7.1 кВт в несколько раз меньше соответствующей величины, приведенной в работах [18, 20], что говорит в пользу предлагаемой новой модели разряда.
Рассмотрим распределение мощности диссипи-рующей в плазме ВЧФР электромагнитной энергии.
Мощность джоулевых потерь, возникающих в результате диссипации электромагнитной волны, запишется следующим образом [20]:
W = 2k\)r\I |2 rdrdz.
(13)
00
Здесь \1 | - действующее значение высокочастотного тока; г - радиус канала ВЧФР; Я - величина активной составляющей сопротивления плазмы
s
н
и
I
1 -
0,25 0,5 0,75 Z=z/L
Рис. 8. Зависимость dW/dz от длины канала разряда
1,0
ВЧФР, приходящейся на единицу длины, которая в случае слабого скин-эффекта определится следующим образом:
2\-1
R = (anr2)
(14)
Используя выражение (5) для высокочастотного тока, распространяющегося по линии конечной длины, получим выражение для действующего значения тока:
<а? ~ I2 I2
= 10
cos Pz + +Ce~az cos(y-Pz)
eaz sin Pz +
+Ce
-az
sin(p-Pz)
(15)
С учетом выражений (14) и (15) интеграл (13) для случая C = 1, р = 0 запишется следующим образом:
2 z
W = 2а"1 |I0| J(ch2az-sin2 pz)dz . (16) 0
Мощность, диссипирующую в канале ВЧФР, найдем путем интегрирования выражения (16).
Получим:
W (L) =
2а"1 |I0
ch 2aL
+
2a
+—— sin BL cos BL - — 4P 2
(17)
Найдем распределение плотности источников диссипации электромагнитной энергии вдоль канала ВЧФР. Мощность, приходящаяся на единицу длины канала ВЧФР, характеризуется величиной
ёШ
-. Из выражения (17) следует, что:
dW dz
2а 1 |l012 (ch2az - sin2 Pz).
Результаты расчетов величины
dW dz
для случая
разряда с длиной канала Ь приведены на рис. 8. Как видно из рис. 8, по длине разряда можно выделить три зоны интенсивного энерговыделения - на краях и в центральной части канала разряда.
Подводя итог вышесказанному, делаем такие заключения:
1. Полученное на основе теоретических расчетов и экспериментальных результатов распределение электрического тока вдоль канала ВЧФР в отличии от работ [18, 20] носит неэкспоненциальный характер.
2
3
2
0
2. Расчеты, проведенные на основе предложенной модели разряда, представленного в виде линии конечной длины с распределенными параметрами, подтверждают предположение о наличии в канале ВЧФР отраженной электромагнитной волны.
3. Отражение электромагнитной волны и тока в конце канала разряда носит характер, близкий к синфазному.
4. Длина электромагнитной волны, распространяющейся вдоль канала ВЧФР, как правило, превышает длину канала разряда, т.е. выполняется
условие Ь < 2пал. 5. Затухание электромагнитной волны в плазме ВЧФР в 5.10 раз меньше, чем это следует из результатов работ [18, 20].
Все это позволяет сделать вывод, что в предлагаемой электродинамической модели лучше учитываются свойства и особенности ВЧ факельного разряда, чем в модели авторов работ [18, 20].
Результаты проведенных исследований и расчетов направлены на более эффективное исследование ВЧ факельного разряда в технологии и технике.
СПИСОК ЛИТЕРАТУРЫ
1. Ъжмиров ИА. // Изв. СО AН СССР. - 1980. - Вып. 2. - № 8. - С. 3-13.
2. Физика и техника низкотемпературной плазмы / Под ред. С.В. Дресвина. - M.: Aтомиздат, 1972. - 352 с.
3. Сорокин ЛЖ // Tеория электрической дуги. - Новосибирск: Наука, 1977. - T. 46, вып. 4. - С. 15-26.
4. Русанов В.Д., Фридман A.A. Физика химически активной плазмы. - M.: Наука, 1984. - 416 с.
5. Hittorf W. // Annal. Phys. Chem. - 1884. - P. 71-74.
6. Зилитинкевич С.И. // Tелеграфия и телефония без проводов. - 1928. - № 9. - С. 20-27.
7. Нейман M.C // Известия электропромышленности слабого тока. - 1935. - № 7. - С. 38-48.
8. Aппаратура и методы исследований ВЧ-разрядов (по материалам семинара ВДНХ 16.01.1973 г.) // Под ред. проф. ИА. ^хомирова. - Ъомск: Изд-во Try, 1976. -64 с.
9. Капцов НА. Электроника. - M.: ЩГГЛ, 1954. - 320 с.
10. Mатвеев Я.Н. // ЖГФ. - 1941. - T. 11. - С. 63-80.
11. Бамберг АС., Дресвин С.В. // ЖГФ. - 1962. - T. 32. -С. 31-40.
12. Trunecek V. // Folia Fac. Sci. - Brno: Nat. Universit. -1971. - V. 12. - Р. 51-54.
13. Mарусин В.В., Ъжмиров ИА // Применение низкотемпературной плазмы в технологии: Tруды Всес. семинара. - Киев: ОНШ AK УССР, 1968. - С. 16-18.
14. Ъжмиров ИА, Лисицын ВЖ // Tруды IV Всес. конф. по физике и генераторам низкотемпературной плазмы. - Aлма-Aта: КазПИ, 1970. - С. 11-12.
15. Ъжмиров ИА, Ъжмиров В.В., Федянин Б.Я. // Изв. ТОИ. - 1972. - T. 225. - С. 110-113.
16. Ъжмиров ИА и др. // Tруды VI Всес. конф. по генераторам низкотемпературной плазмы. - Фрунзе: KraM, 1974. - С. 8-10.
17. Ътаомиров ИА. Луценко Ю.Ю. // Изв. СО AK
СССР, сер. техн. наук. - 1989, вып. 1. - С. 81-84.
18. Качанов А.В., Трехов Е.С., Фетисов Е.П. // ЖТФ. -1970. - Т. 15. - С. 340-345.
19. Тихомиров И.А., Власов В.А. // Изв. вузов. Физика. -2000. - № 4. - С. 143-150.
20. Качанов А.В., Трехов Е.С. // Физика газоразрядной плазмы. - М.: Атомиздат, 1968. - Вып. 1. - С. 60-67.
21. Тихомиров И.А. и др. // Изв. ТПИ. - Томск. - 1976. -Т. 276. - С. 40-43.
22. Тихомиров И.А. и др. // Труды V Всес. конф. по физике низкотемпературной плазмы. Т. 2. - Киев: ФОЛ АН УССР, 1979. - С. 11-15.
23. Тихомиров И.А., Луценко Ю.Ю. // Изв. СО АН СССР, сер. техн. наук. - 1989. - Вып. 3. - С. 109-115.
24. Тихомиров И.А., Луценко Ю.Ю. // ЖТФ. - 1989. - Т. 59. - Вып. 11. - С. 128-130.
25. Уолтер К. Антенны бегущей волны. - М.: Связь, 1970.
- 372 с.
26. Анго А. Математика для электро- и радиоинженеров.
- М.: Наука, 1955. - 778 с.
27. Луценко Ю.Ю., Тихомиров И.А., Шайхеев А.Г., Ко-рючкин А.В. // Плазменная техника, технология. -Казань, 1988. - С. 5-6.
28. Тихомиров И.А., Луценко Ю.Ю. // Тез. докл. V Всес. конф. по физике газового разряда. - Омск, 1990. - С. 81.
29. Луценко Ю.Ю., Тихомиров И.А. // Тез. докл. Всес. семинара по высокочастотному пробою газов. - Тарту, 1989. - С. 74.
30. Тихомиров И.А., Тихомиров В.В., Левашов В.С. // Изв. ТПИ. - 1976. - Т. 276. - С. 12-15.
31. Стрэттон Д. Теория электромагнитизма. - М.-Л.: Гос-техиздат, 1948. - 548 с.
32. Тихомиров И.А., Власов В.А., Луценко Ю.Ю. Физика и электрофизика высокочастотного факельного разряда и плазмотроны на его основе. - М.: Энергоа-томиздат, 2002. - 200 с.
