ЭНЕРГЕТИКА • ЭЛЕКТРОТЕХНИЧЕСКИЕ КОМПЛЕКСЫ И СИСТЕМЫ
УДК 537:621.3
В. М. КУЛЯПИН, И. М. АСЛЯМОВ КАТОДНЫЕ ПРОЦЕССЫ ЭЛЕКТРИЧЕСКИХ РАЗРЯДОВ
Предложена математическая модель нестационарного процесса фазовых превращений: плавления с уносом массы мишени под действием поверхностного источника энергии высокой плотности. Установлено, что исходные уравнения баланса энергии имеют непрерывные решения, но при выделении асимптотик возникает скачкообразность решения. Получены значения параметров области катодного падения потенциала. Процессы фазовых превращений; испарения; катодные процессы; электрический разряд
ВВЕДЕНИЕ
Моделирование катодных процессов позволяет вскрыть закономерности и параметры области катодного падения потенциала, в том числе такие, которые пока невозможно определить экспериментально: напряженность поля, плотность тока электронов, ионов и другие.
Предлагается математическая модель для нахождения нестационарного процесса фазовых превращений плавления с уносом массы мишени испарением под действием поверхностного источника энергии высокой плотности с учетом объемного тепловыделения. Получено аналитическое решение задачи теплопроводности с фазовыми превращениями. Установлено, что исходные уравнения баланса энергии имеют непрерывные решения, но при выделении асимптотик возникает скачкообразность решения, их критичность к малому изменению параметров.
В расчетах использовано решение задачи сверхбыстрого нагрева, плавления и разрушения материалов. Интегральный метод решения таких задач в ряде случае дает решение более точное, чем использование численного метода конечных разностей. На основе интегрального метода задача сводится к модели, позволяющей найти решение в замкнутой аналитической форме и провести сопряжение с газокинетическими процессами испарения.
При решении приняты допущения:
• теплофизические характеристики материала мишени постоянны и равны средним значениям для рассматриваемого диапазона температур;
• рассматриваемые процессы подчиняются закону Фурье.
В работе исследованы катодные процессы электрических разрядов различной длительности, проведено совместное рассмотрение эмиссионных и теплофизических процессов на границе «металл-плазма». Плотность тока электронов определена уравнением тока термоавто-электронной эмиссии Мерфи и Гуда, соотношение ионной и электронной составляющих подчиняется уравнению Маккоуна. В решение введены только значения работы выхода электронов, катодное падение потенциала и теплофизические постоянные материала.
Достоверность результатов подтверждена на основании сравнения с экспериментальными значениями, полученными Г. А. Месяцем. Расчетные величины хорошо согласуются с экспериментальными исследованиями.
1. МАТЕМАТИЧЕСКИЕ МОДЕЛИ ВЗАИМОСВЯЗАННЫХ ПРОЦЕССОВ ПРИ ВОЗДЕЙСТВИИ ЭЛЕКТРИЧЕСКОГО РАЗРЯДА
Эффекты в электрических элементах и устройствах систем управления связаны со сложными взаимодействиями электродинамических, термодинамических, гидродинамических процессов во всех агрегатных состояниях вещества - твердых телах, жидкостях, газах, плазме. Законченные физико-математические постановки здесь отсутствуют. Поэтому математическое моделирование является одним из действенных методов исследований.
В работе излагаются общие принципы анализа физического подобия и их приложения к процессам фазовых превращений. Широкое применение плазменных, ионных, электронных и искровых технологий выдвигает задачу исследования нелинейных процессов в электрических разрядах. Такие процессы характеризуются разнообразием разрядов в газе: неса-
Контактная информация: (347) 273-73-43
мостоятельный разряд под действием внешнего ионизатора и самостоятельные газовые разряды. Все виды разрядов: дуговой, искровой, коронный, тлеющий, возникающий при низких давлениях, имеют область катодного падения потенциала, которая обеспечивает механизм извлечения электронов из металла.
Такие процессы характеризуются высокой концентрацией мощности, малыми размерами области взаимодействия, малой продолжительностью, что затрудняет теоретическое и экспериментальное исследование традиционными методами. В этих условиях весьма перспективно применение методов математического моделирования.
В статье проведено рассмотрение сложных многопараметрических процессов в области катодного падения потенциала электрических разрядов, включающие теплофизические, эмиссионные и кинетические явления.
На основе решения уравнения баланса энергии вскрыты закономерности протекающих дискретных процессов, которые послужили основой разработки экспериментальных методов исследования и методик расчета ионноплазменных и других нанотехнологий с устройствами, в которых исполнительным органом является электрический разряд.
Показана действенность применения математического моделирования при исследовании многопараметрических взаимосвязанных нелинейных процессов. Анализ результатов проведен с широким привлечением экспериментальных данных.
2. МОДЕЛЬ ПРОЦЕССОВ НА КАТОДЕ ЭЛЕКТРИЧЕСКИХ РАЗРЯДОВ
Рассмотрим физические явления, возникающие на границе плазмы электрического разряда и поверхности катода. Модель физических явлений, возникающих при воздействии излучения электрического разряда, схематично представлена на рис. 1.
Рис. 1. Модель переходных процессов фазовых переходов
Тепловой поток бомбардировки катода ионным током плотностью дг(0 с учетом ионной температуры плазмы и энергии возбуждения рекомбинированных ионов воздействует на поверхность. Температура поверхности начинает повышаться. После того, как она достигает температуры первого фазового превращения, начинается движение границы фазового перехода Х(0 с затратами на покрытие скрытой теплоты фазового перехода Ь.
3. ВЗАИМОСВЯЗАННЫЕ ПРОЦЕССЫ НА КАТОДЕ ЭЛЕКТРОРАЗРЯДНЫХ УСТРОЙСТВ
3.1. Параметры микропроцессов на границе «металл-плазма»
Проведем совместное рассмотрение эмиссионных и теплофизических процессов в элементарном неподвижном катодном пятне на основе решения уравнения баланса энергии на катоде. В решение вводим только экспериментальные значения работы выхода ф, катодное падение потенциала ик и теплофизические постоянные катода - скрытая теплота парообразования Ь0 и плавления Ь, удельная теплоемкость с и температура плавления Тм.
Температурные режимы катодного пятна рассчитываем с учетом плавления и испарения материала катода поверхностным и объемным источниками теплового потока. Элементарные процессы рассматриваем в динамике. Так как инерционность теплофизических процессов намного выше инерционности эмиссионных процессов, учитываем зависимости от времени только теплофизических процессов. Учитываем энергию резистивного нагрева только расплавленной зоны, размер которой определяется величиной и длительностью воздействия теплового потока. С учетом принятых допущений решается уравнение баланса энергии в элементарном катодном пятне:
д + (х - хо) = до + де + дг. (1)
В уравнении учитываются источники энергии - бомбардировка катода ионным током дь резистивный нагрев расплавленной зоны металла ду и пути рассеяния энергии - термоэмиссию электронов де, испарение, плавление и теплопередачу в катод д0. Радиационные потери дг в рассматриваемом диапазоне температур незначительны и не учитываются при расчетах. Тепловой поток бомбардировки катода ионным током равен
дг = т(1 + Л + п) Ще ик , (2)
где ц - коэффициент аккомодации; п - коэффициент, учитывающий влияние ионной температуры плазмы; V - коэффициент, учитывающий энергию возбуждения рекомбинированного иона.
Соотношение между ионной / и электронной /е составляющими тока разряда у = / / /е определяется уравнением Маккоуна для тока амбиполярной эмиссии
/ \
Е 2
1 + -
у = -
1
т,
1
12
V те
1
7,57 • 105 • ик2 • ]е
(3)
где
тг / те - отношение массы иона к массе электрона; Е - напряженность поля у катода.
Плотность тока электронной эмиссии ]е для высокотемпературной области катодного пятна описывается уравнением термоавтоэлектрон-ной эмиссии Мерфи и Гуда
]е _ 12070'
1,64 • 10-2 Е3/4 / Т0 біп(1, 64 • 10-2 Е3/4 / Т0)
X
х ехр
^4,39 Е1/2
1,16 • 104 фЛ
(4)
70
V и У
где 70 - температура поверхности катодного пятна.
Тепловой поток резистивного нагрева расплавленной зоны металла, отнесенный к единице площади катодного пятн, равен
2
ду _ 0,48(4^ гк 3)^-;
3 а
2
ду = 0,482^- (X - X0)3;
(5)
3а
где а_
1 + 3 ат (Т0 - Т)
электропроводность
катода; ат - электропроводность при температуре плавления; ат - температурный коэффициент; у = ]е+]г =(1+у)/е - плотность тока в катодном пятне.
Глубина плавления Х - Х0 определяется из [1-5]. Тепловой поток, отводимый термоэмиссией электронов, равен
д = (ф - 3,810”4 -Еш + 1,724-10-4 То)]е. (6) Температура катодного пятна в зависимости от плотности теплового потока определяется уравнением
_11(70 - тт) ,
д0 _----------------+
0 X - Х0
[ І0 + ь+с(7) - т )1
м л в ч
— ехр(А - 4,234--------),
т0 т0
(7)
где М - атомная масса материала катода, А, В -табулированные постоянные скорости испарения.
По приведенным уравнениям рассчитаны значения параметров области катодного падения потенциала, удовлетворяющие условия равновесия теплофизических и эмиссионных процессов на катодах из серебра, меди, железа, никеля и золота. Результаты расчетов представлены на рис. 2-5.
10( -2 5 10‘ 2 Е Е/см
Рис. 2. Равновесные параметры разрядов на серебряном электроде
на медном катоде
г
на стальном катоде
510- 10е 210е 510е 10' 210' Е В/см
Рис. 5. Равновесные параметры разрядов на никелевом катоде
Выбраны материалы с однократной ионизацией, чтобы полнее выявить механизм взаимосвязанных процессов..
Из рассмотрения результатов расчета можно выделить три характерные области катодного падения потенциала электрических разрядов. Рассмотрим их на примере серебра.
При напряженности Е > 2,6-107В/см (рис. 2) длительность существования катодных пятен меньше времени нагрева поверхности катода до температуры Т0 и начала испарения, термоэмиссия электронов ослабляет поверхностный источник тепла и, по-видимому, основным источником энергии является энергия резистивного нагрева. При напряженности Е > > 2,6-107 В /см доля энергии поверхностного теплового потока снижается, и основным источником энергии становится энергия резистивного нагрева материала катода. Процессы
также дискретны. Температура поверхности снижается, извлечение электронов подчиняется автоэлектронному механизму эмиссии.
Рис. 6. Равновесные параметры разрядов на катоде из золота
Длительность существования элементарных катодных пятен в функции напряженности поля у катода представлены на рис. 7.
Рис. 7. Длительность существования элементарных катодных пятен
Следовательно, процессы на катоде в зависимости от времени определяются различными источниками энергии и подразделяются на три вида. Дискретные процессы, когда основную роль играет тепло резистивного нагрева, дискретные процессы, определяемые энергией бомбардировки катода ионным током, а тепло Джоуля-Ленца приводит к прекращению разряда. Стационарные процессы, когда баланс энергии выполняется независимо от длительности разряда и величины расплавленной зоны.
3.2. Классификация катодных процессов
В соответствии с особенностями протекания теплофизических и эмиссионных процессов катодные процессы электрического разряда можно разделить на три типа. Следует отметить относительную условность границ разбиения, так как они определяют только преобладающий
механизм развития пятен, которые существуют на катоде одновременно и постоянно вырождаются одно в другое, возникают вновь.
Катодные пятна в области напряженности поля (для серебра) Е > 2,6-107 В/см (рис. 2) дискретны, характеризуются преобладающей ролью объемного источника нагрева расплавленного объема, малой продолжительностью существования (меньше 1-30 мкс) и высокой плотностью электронного тока(/ > 2 106 А/см2). Катодные пятна в этой области напряженности поля относятся к первому типу пятен и характерны для искровых разрядов и начальной стадии электрического пробоя. Примем за максимальную продолжительность существования пятен первого типа время нагрева поверхности катодного пятна до температуры Т0 и начала интенсивного испарения
т =__________2_________ (9)
'0 - 2. у ’
Ц + Ц + с(Т0 ~ т)
_ Ц0 + с(Тт - Т) _
Минимальная продолжительность пятен первого типа определяется начальными условиями разряда и максимальной напряженностью поля у катода. Примем за верхнюю границу пятен первого типа время '1 , при котором нормированная функция распределения длительности существования дискретных элементарных катодных пятен равна единице. То есть нижний предел интеграла от функции распределения
2
| £ (' )4Т = 1. (10)
к
Интегрирование ведется до времени ' =2, при котором прекращается влияние объемного источника тепла, и катодные пятна вырождаются в неподвижные. Оценим предельные значения параметров пятен первого типа. Используем функцию распределения продолжительности существования элементарного катодного пятна
£( К) = 0,1/ Г (11)
Тогда из условия нормирования функции распределения продолжительности существования катодных пятен определяется нижний предел
к
І!= ехр(-9,306853) = 9,1 10-5.
Напряженность поля при данной продолжительности существование определяется
Е = Ег + 9,3 07 - 106 В/см. (12)
Из решения уравнения баланса энергии определяем значения параметров катодных пятен ]е, ], Е для времени = 9.1-10-5. Верхние
значения параметров катодных пятен первого типа представлены в табл. 1.
Т аблица 1
Параметры процессов первого типа
Парамет- ры Металл Е-107, В/см je■106, А/ см2 ;.-106, А/см2 4-10-9, с
Серебро 3,43 10,0 1,0 1,13
Медь 2,43 6,36 0,775 94,2
Золото 2,33 1,75 0,285 250,5
Железо 2,93 4,75 0,86 0,6
Никель 2,83 2,34 0,807 1,7
Дискретные катодные пятна второго типа соответствуют длительности существования от '0 до ' = 2 и области напряженности поля 2-107 < Е < 2,6-107 В/см (для серебра). Эмиссионные процессы катодных пятен 2-го типа определяются в основном энергией поверхностного источника бомбардировки катода ионным током. Резистивный нагрев расплавленной зоны приводит к распаду неподвижного элементарного катодного пятна. По мере увеличения расплавленной зоны резистивный нагрев приводит к тому, что через время 'ь суммарный тепловой поток, подводимый к пятну, начинает превышать величину, которую способны отвести эмиссия электронов и теплофизические процессы. Наличие избыточной энергии вызывает перегрев расплавленного металла, взрыв и распад элементарного катодного пятна. Катодный процесс перемещается в новую точку привязки разряда. Катодные пятна 2го типа в основном определяют эрозию катода и самопроизвольное перемещение катодных пятен. Граница пятен 1-го и 2-го типов довольно четко определена как по длительности существования катодных пятен, так и по величине напряженности электрического поля у катода. Минимальное время существования катодных пятен 2-го типа определяется временем нагрева материала катода до температуры Т0 и начала испарения с поверхности. Минимальное время существования элементарного катодного пятна второго типа определяется уравнением
t , _ 2[L0 + с(Т0 - Тт )]Р£tb min с (13)
lbmin _ о ’С, ( )
qo
Т аблица 2 Параметры процессов второго типа
Мате- риал Ag Cu Au Fe N1
Е • 107, B/см 2,5 1,5 1,4 2,0 1,9
Т0, К 3901,5 3900,0 4079,6 4392,1 4504,6
je-106, А/см2 2,11 1,16 0,47 1,73 0,88
У = j/je 0,22 0,19 0,23 0,24 0,42
Q v - 106, Вт/см2 9,67 4,74 2,41 8,61 7,33
qe-106, Вт/см2 6,51 3,71 1,75 6,16 3,46
qv (X-X0), Вт /см2 0,84 0,61 0,34 1,40 1,87
j - 106, А/см2 0,47 0,22 0,66 0,43 0,37
tb* 0,03 0.03 0,019 0,051 0,05
tb -106, с 3,78 31,4 52,3 0,34 0,94
fmax , кГц 265 31,8 19 290 1060
В табл. 2 приведены значения минимального времени существования катодных пятен 2-го типа и соответствующие значения параметров области катодного падения потенциала. Как видно из табл. 2, для различных катодов минимальная продолжительность существования катодных пятен 2-го типа (и, соответственно, максимальная продолжительность существования пятен 1-го типа) находится в пределах
0,3-30 мкс.
Нижняя граница пятен 2-го типа определяется временем достижения расплавленной зоной металла установившегося значения, а по напряженности соответствует минимальной напряженности поля, при которой наблюдается распад элементарного катодного пятна. Установившееся значение расплавленной зоны достигает за время tв тах = 2, или
4тах = 2[Ьо + с(То - Тт)\ р£ /^о^ С. (14) Значения параметров области катодного падения потенциала для пяте второго типа при 4мах = 2 приведены в табл. 3.
Таблица 3
Минимальные значения параметров дискретных процессов
Материал катода Ag Cu AU Fe N1
Е -107 , В/см 2,0 1,0 0,90 1,5 1,4
То, K 3895,3 3778,1 3906,2 4300,0 4356,2
je-106 , А/см2 1,02 0,34 0,12 0,71 0,31
У _j/je 0,30 0,28 0,37 0,33 0,65
6, 2 2 | Н q В 6,17 2,10 0,99 4,82 3,97
qe- 106,Вт/см2 3,34 1,19 0,48 2,69 1,30
qv (Х-Х0И06, Вт/см2 1,11 0,30 0,14 0,98 1,39
jг 106, А/см2 0,30 0,09 0,04 0,24 0,20
Т BMAX '10 ,С 0,20 1,97 2,02 0,003 0,03
f, кГц 4,9 0,50 0,49 39 27,5
Из сравнения данных табл. 2 и 3 следует, что продолжительность существования катодных пятен второго типа находится в пределах от 0,3-50 мкс до 0,03-2 мс. Границы пятен второго типа нанесены на рис. 2-6 пунктирной линией. На основании расчетов получено, что для пятен второго типа характерно постоянство результирующего потока, расходуемого на теплофизические процессы q0 » const, и, соответственно, температуры поверхности Т0.
При напряженности поля Е < 2-107 В /см на катоде образуются неподвижные пятна третьего типа. Даже при максимальной величине расплавленной зоны металла в пятне наблюдается равновесие теплофизических и эмиссионных процессов. Эрозия пятен происходит только за счет испарения материала катода. Максимальная плотность тока электронов находится в пределах (0,12-1,022)-106 А/см 2. В зависимости от напряженности поля изменяется суммарный тепловой поток q0 = qг - qe, при этом изменяется и температура катодных пятен. Определим нижнюю границу существования катодных пятен третьего типа по напряженности как точку, при которой функция распределения плотности тока электронов равна нулю (при постоянстве температуры катодного пятна)
g (je) = 0. (15)
Нижние границы существования стационарных пятен 3-го типа приведены в табл. 4.
Таким образом, из расчета баланса энергии на катоде определены пределы различных механизмов взаимодействия материала катода с плазмой, на основании чего проведена классификация типов катодных процессов.
Т аблица 4
Параметры неподвижных точек привязки разряда
Материал катода АВ си Аи Ге N1
Е -107, В / см 0,6 0,275 0,4 0,6 1,1
]е ■ 103, А/см2 2,697 3,65 6,0 25,4 87,0
./'■ ■ 103, А/см2 2,5 7,6 8.7 35,5 133
Границы пятен всех трех видов нанесены на рис. 2-6 вертикальными линиями. В стационарных процессах перенос тока осуществляется пятнами всех трех типов, образующих квази-стационарные макропятна. Для характеристики интегральных параметров процесса необходимо получить функции распределения.
ЗАКЛЮЧЕНИЕ
1. Из совместного рассмотрения в динамике процессов на катоде установлено, что катодные процессы дискретны. Распад элементарных катодных пятен вызван превышением скорости выделения энергии над скоростью теплопередачи испарением, теплопроводностью и термоэмиссией электронов, дискретность процессов является причиной высокочастотных колебаний катодного потенциала.
2. Установлено, что характеристики процессов в аппаратах с электрическими разрядами могут быть определены на основе решения уравнения баланса энергии на катоде. Учитывается энергия бомбардировки катода ионным током и резистивный нагрев расплавленной зоны катода. Отвод энергии осуществляется термоэмиссией электронов фазовыми преобразованиями.
3. Установлено, что различные виды электрических разрядов определяются соотношением катодных пятен трех типов. Катодные процессы дискретных пятен первого типа характерны для импульсных разрядов с длительностью меньше времени нагрева катода до температуры испарения, когда основным источником является энергия резистивного нагрева.
4. Катодные процессы в дискретных пятнах второго типа определяются энергией бомбардировки катода ионным током, резистивный нагрев не превышает 40% суммарной энергии. Но резистивный нагрев приводит к развитию неустойчивости катодных процессов, когда
скорость подвода энергии начинает превышать скорость отвода термоэмиссией электронов, теплофизическими процессами происходит перегрев и взрыв расплавленного катодного пятна. Неподвижные пятна третьего рода характерны для процессов с низкой напряженностью поля, в них соблюдается равновесие теплофизических и эмиссионных процессов.
5. Установлено, что неустойчивость катодных процессов определяется микропроцессами в элементарном катодном пятне, перегревом и распадом точки привязки разряда
6. Показано, что при исследовании катодных процессов необходимо применять статистические методы.
СПИСОК ЛИТЕРАТУРЫ
1. Куляпин В. М. Некоторые задачи теплопроводности с фазовыми превращениями // Инж.-физ. журнал. 1971. Т. 20, № 3. С. 497-504.
2. Куляпин В. М., Аслямов И. М. Электрический разряд в системах электрооборудования летательных аппаратов // Электрические контакты и электроды: Тр. ин-та проблем материаловедения НАН Украины. Киев, 2008. С. 23-31.
3. Месяц Г. А. Электронная эмиссия из сегнето-электрических плазменных катодов // Успехи физических наук. 2008. Т. 178, № 1. С. 85-108.
4/ Месяц Г. А., Яландин М. И. Пикосекундная электроника больших мощностей // Успехи физических наук. 2005. Т. 175, № 3. С. 226-245.
5. Месяц Г. А. Эктоны в вакуумном разряде: пробой, искра, дуга. М.: Наука, 2000.
ОБ АВТОРАХ
Куляпин Владислав Максимилианович, проф., каф. электрооборуд. ЛА. Дипл. инж.-электромех. (Новосиб. электротехн. ин-т, 1958).
Д-р техн. наук по по эл-там и уст-вам систем упр-я (УГАТУ, 2003). Иссл. в обл. систем упр-я с электрич. разрядами.
Аслямов Ирек Мирзафа-тихович, асп. той же каф. Дипл. инж. по электрооборуд. ЛА (УГАТУ, 2005). Иссл. в обл. систем регулир. электрооборуд. ЛА.