CC BY
10Journal of Siberian Federal University. Engineering & Technologies 2 (2013 6) 183-191 УДК 533.95
Расчет границ радиационного и теплопроводного режимов дугового разряда
Е.Н. Васильев*
Институт вычислительного моделирования СО РАН, Россия 660036, Красноярск, Академгородок, 50/44
Received 15.03.2013, received in revised form 22.03.2013, accepted 31.03.2013
По результатам расчета в осесимметричной постановке структуры, удельных и интегральных мощностей энергетических механизмов стационарного равновесного дугового разряда определены границы радиационного и теплопроводного режимов горения разряда в воздухе, углекислом газе и аргоне.
Ключевые слова: дуговой разряд, энергетический баланс, излучение, теплопроводность.
Введение
При постановке и выполнении исследовательских задач, практическом использовании дуговых разрядов может быть полезна информация о параметрах газоразрядной плазмы и степени воздействия различных физических факторов на ее состояние. На параметры плазмы наибольшее влияние оказывают энергетические механизмы, от вклада которых зависят распределения температуры и электропроводности, другие характеристики разряда. В свою очередь, соотношение мощности энергетических механизмов определяется параметрами процесса: силой тока, давлением, геометрическими размерами, транспортными и радиационными свойствами рабочего газа. Для случая стационарной дуги, ограниченной цилиндрической стенкой, определенного газа и заданного давления характеристики разряда зависят прежде всего от силы тока I и радиуса стенки R0. В данной статье по результатам расчетов структуры разряда и характеристик энергетического баланса в широком диапазоне значений силы тока и радиуса стенки определены на плоскости I-R0 границы радиационного и теплопроводного режимов и соответствующие им значения максимальной температуры для дугового разряда в воздухе, аргоне и углекислом газе. Работа является развитием исследований [1, 2], в которых рассмотрены особенности формирования структуры дугового разряда в воздухе и определены границы радиационного режима.
Постановка задачи
Рассмотрим электродуговую плазму, находящуюся в состоянии локального термодинамического равновесия, при атмосферном давлении данное приближение применимо при силе тока I > 10 А [3]. Тепловыделение в электропроводном газе обеспечивается за счет джоулевой
© Siberian Federal University. All rights reserved
* Corresponding author E-mail address: [email protected]
диссипации, основными механизмами отвода теплоты от столба разряда являются излучение и теплопроводность. Конвекция, как известно, не оказывает существенного влияния на энергобаланс сильноточной дуги, ограниченной стенкой. Как показали данные измерений, вклад конвекции в ртутной дуге ат мо сферного давления не превышал 1 % [4] .
Распределение температуры в разряде, горящем в цилиндрическом объеме и имеющем осевую симметрию, описано уравнением Эленбааса-Хеллера следующего вида:
1 д С ^дт\ 0
1 + ^ -Чп = 0 (!)
г дг\ дг )
где г - ртдиус, Г аа температуиа, ¿н = оЕ2 - мощность джоулевой диссипнци и, дк - объемная мощность радиационных энергопотерь, К и а -коэффицитнты твплопроаодности и электро-проводноати раза, Е - напряженность электричеокого поля.
Ураененио (1) дополняетоя граничными сословиями в центре разрядное области
сТ
= 0 (2)
r=0
дг
и на внешней границе
[т]нЩ=т°. (3)
Напряженность электрического поля зависит от силы тока и распределения электропроводности по радиусу
Ro
E = I Пп. J ardr. (4)
о
Элоктрические дуги, как правило, имеют поперечный размер, не превышающий нескольких сантиметров, и малую оптическую толщину. В этом случае излучение имеет преимущественно объемныай характер. Так, в рзботе [5] бышо показано, что учут перзноса излучения в виде лучистой теплопроводно сти начинает оказытвать какоа-либо заметнот влияние на характеристики разряду е аргоне приз атмосферном давлении только при температурах свыше 15 000 К С учотом итого в данной математической модели для расчета радиауио нных энергопотерь было исполь зовано приближение обвамного излучателя [1].
Раальныае аазы имеют сложные нелинейною темпераоурныое зависимо сти теплофизических свойств, что обусловливает необходимоств применения численшгх меоодов. Решение краевой задачи (н)-(4) прсводилось с помощью конечноразносаного алгоритма, величина пространственно го шага в расчетах принималась равной 0.5 мм. Коэффициенты черноты аезду^ха взяты в [6],
[7], аргона к [5], углекислого газа в [7] , транспортные свойства воздуте в [8]. аргона в [9], угле-нислого газо в [10]. Температурные завис;имо сги сво йств газов вводились в программу в виде тоулиц для степеней черноты е( 7), коэффициентов эоектропроводности c(T) и 'теплопроводности Х(ТУ с максимальным значением температуры Г =20 кК (для аргона 17 кК) и шагом AT = 1 кК.
Результаты расчетов
Вычислительная модель позволяет рассчитать структуру разряда для разных газов при различных значениях силы тока I и радиуса расчетной области R0. Из решения уравнений
(1) - (4) определялось распределение температуры и рассчитывались напряженность электрического поля Е, удельные и интегральные энергетические характеристики дугового столба. Напряженность электрического поля является важным параметром, характеризующим энергобаланс в разряде, поскольку она определяет величину мощности джоулевой диссипации, компенсирующей теплопроводные и радиационные энергопотери. На рис. 1 и 2 приведено сопоставление вольтамперных характеристик (ВАХ), полученных из экспериментов и с помощью данной математической модели. На рис.1 представлены расчетные и экспериментальные [11] значения Е при атмосферном давлении в воздухе (свободно горящая дуга) и углекислом газе (радиус канала 2 см). На рис. 2 для аргона при атмосферном давлении и двух значений радиуса канала (2.5 и 3 мм) сопоставлены данные измерений [12] и расчетов, которые в этом случае были проведены при величине разностного шага 0.1 мм. И эксперимен-
Рис. 1. Экспериментальные ВАХ и расчетные значення Е (квадратики для воздуха, крестики для СО2)
15
5
о т то но гаи
1,8
Рис. 2. Экспериментальные ВАХ и расчетные значения Е для аргона (кривая 1 и крестики для К0 = 2,5 мм, кривая 2 и квадратики для К0 = 3 мм)
тальные и расчетные данные свидетельствуют о достаточно сильной зависимости величины Е от рода газа и параметров разряда.
Для определения границ режимов дугового разряда рассчитывался энергетический баланс разряда в локальных и интегральных величинах. В первом случае сравнивалось соотношение удельных мощностей механизмов излучения qR, теплопроводности и джоулевой диссипации qJ в центре разряда. Во втором - определялись суммарные по всему объему разряда мощности лучистых энергопотерь QR и джоулева тепловыделения QJ, а также тепловая мощность Qт, уходящая из разряда за счет механизма теплопроводности через площадь боковой стенки. Следует отметить, что в расчетах с высокой точностью выполнялись балансовые соотношения как интегральные QJ = QR + Qт, так и локальные qJ = qR + qт во всех расчетных точках. По удельным и интегральным энергетическим характеристикам рассчитывались отношения кь = qR/qJ и к = QR/QJ, в качестве критериев при определении границ приняты значения к = 0,1, 0,5 и 0,9. При значении к < 0,1 влияние излучения в энергобалансе пренебрежимо мало, при к > 0,9, наоборот, доминируют радиационные энергопотери. Значение к = 0,5 соответствует равенству вкладов теплопроводности и излучения. Из решения уравнений (1)-(4) для указанных значений к определены величины силы тока Г и максимальной температуры в центре разряда Ттах.
Результаты расчетов для воздуха при атмосферном давлении приведены на рис. 3 и 4. На рис. 3 для значений к = 0,1, 0,5 и 0,9 отображены зависимости Г^о), каждая из которых является геометрическим местом точек, соответствующих режимам с силой тока Г в диапазоне значений R0 = 5 ^ 50 мм. Сплошные линии 1, 2 и 3 на рис. 3-10 соответствуют значениям к: = 0,9, 0,5 и 0,1, штриховые линии 4, 5 и 6 относятся к аналогичным значениям кь. Кривые 2 и 5, для которых выполняются равенства QR = QT и qR = qт, делят плоскость R0 - Г на две области, выше кривых находятся режимы, в энергобалансе которых вклад излучения больше, ниже - режимы с более сильным влиянием теплопроводности. Выше границ 1 и 4 расположены режимы, в энергобалансе которых доминирует излучение QR >> Qт и qR >> qт, соответственно, теплопро-
4.5
4
3
2
5 10 15 20 25 30 35 40 45 50
R, мм
Рис. 3. Зависимости lgI*(R0) для воздуха приp = 0,1 МПа
°5 10 15 20 25 30 35 40 45 50
R, мм
Рис. 4. Зависимости Tmax(R0) для воздуха при p = 0,1 МПа
4 X ■
V S -
&
1
5 10 15 20 25
R, мм
Рис. 5. Зависимости ^Г(Г0) для воздуха прир = 1 МПа
водные режимы, когда QT >> QR и qT >> qR, находятся ниже кривых 3 и 6. Графики значений максимальной температуры Ттах(Г0) в центре разряда, соответствующие полученным зависимостям Г(Г0), приведены на рис. 4. Зависимости Ттах(Г0) для всех значений к имеют монотонно падающий характер. Явно выраженный излом кривых 6 на рис. 3 и 4 связан с наличием в зависимости Х(Т) для воздуха локального максимума Хтах = 4,2 Вт/(м*К) при Т = 7 кК (при Т = 6 кК и 8 кК величина X = 2,4 Вт/(м*К)).
На рис. 5 и 6 приведены зависимости Г(Г0) и Ттах(Г0) для воздуха при давлении р = 1 МПа. Кривые на графиках построены в диапазоне значений Г0 = 5 ^ 25 мм, при расчете структуры было установлено, что увеличение давления с 0,1 МПа до 1 МПа приводит к уменьшению по- 187 -
R, мм
Рис. 6. Зависимости Tmax(Ro) для воздуха при p = 1 МТПа
5 10 15 20 25 30 35 40 45 50
R, мм
Рис. 7. Зависимости (Г0) для аргона
перечного размера разрядного столба примерно в 2 раза. Изменение давления в наибольшей степени влияет на излучательную способность газа, при повышении давления пропорционально растет мощность радиационных энергопотерь, по этой причине переход к радиационному режиму наблюдается при существенно меньших значениях тока и температуры. Изломы на кривых 1 и 5 также являются следствием локальных максимума (Хтах = 3,9 Вт/(м*К) при Т = 8 кК) и минимума (ХтП = 1,7 Вт/(м*К) Т = 10 кК) в зависимости Х(Т) для воздуха при давлении р = 1 МПа.
В отличие от молекулярных газов аргон имеет сравнительно невысокую излучательную способность, коэффициент теплопроводности также существенно ниже (особенно для темпе-
- 188 -
2>С|-------------------------1-----------<------------1-----------г
U5 10 15 20 25 30 35 40 45 50
R, мм
Рис. 8. Зависимости Tmax(R0) для аргона
UJ5 10 15 20 25 30 35 40 45 50
R, мм
Рис. 9. Зависимости Igl(R0) для углекислого газа
ратур T < 10 кК) и зависимость Х(Т) имеет более ровный характер. Границы режимов для аргона при атмосферном давлении приведены на рис. 7 и 8. Отличительной особенностью графиков (за исключением кривой 1 на рис. 8) является достаточно слабая зависимость от радиуса.
Для углекислого газа следует отметить близкую с воздухом по характеру и значениям зависимость ЦТ) и существенно более высокие значения коэффициента черноты, коэффициент электропроводности у СО2 выше при температурах T < 11 кК, далее соотношение становится обратным. Зависимости Г (R0) и Tmax(R0) для углекислого газа при давлении p = 0,1 МПа отображены на рис. 9 и 10. К особенностям интегрального энергобаланса для СО2 можно отнести то, что значениям R0 > 32 мм соответствуют только радиационные режимы.
5 10 15 20 25 30 35 40 45 50
R, мм
Рис. 10. Зависимости Tmax(^Q) для углекислого газа
Заключение
Таким образом, рассчитанные зависимости дают возможность оценивать степень влияния энергетических механизмов и величину максимальной температуры в дуговом разряде при различных условиях. В статье приведены два варианта расчета границ радиационного и теплопроводного режимов: в удельных и интегральных величинах. Сложно однозначно определить, какой из них является более предпочтительным для анализа энергобаланса. С одной стороны, соотношение интегральных мощностей наиболее полно отражает энергобаланс всего разряда в целом. С другой - именно соотношение удельной мощности энергетических механизмов в центре разряда устанавливает масштаб температуры, от которого в значительной степени зависят многие, в том числе и интегральные характеристики разряда. По-видимому, в каждом конкретном случае при анализе энергетического баланса следует исходить из цели исследования и специфики поставленной задачи.
Список литературы
[1] Васильев Е.Н. // Теплофизика и аэромеханика. 2010. Т. 17. № 3. С. 441.
[2] Васильев Е.Н. // Журнал технической физики. 2012. Т. 82. № 12. С. 47.
[3] Грановский В.Л. Электрический ток в газе. Установившийся ток. М.: Наука, 1971. 543 с.
[4] Kenty C. // Journal of Applied Physics. 1939. V.10. № 10. P. 714.
[5] Ветлуцкий В.Н., Онуфриев А.Т., Севастьяненко В.Г. // Журнал прикладной механики и технической физики. 1965. № 4. С. 71.
[6] Авилова И.В., Биберман Л.М., Воробьев В.С. и др. Оптические свойства горячего воздуха. М.: Наука, 1970. 320 с.
[7] Каменщиков В.А., Пластинин Ю.А., Николаев В.М., Новицкий Л.А. Радиационные свойства газов при высоких температурах. М.: Машиностроение, 1971. 440 с.
[8] Соколова И.А. // Журнал прикладной механики и технической физики. 1973. № 2. С. 80.
- 190 -
[9] Devoto R.S. // Physics of Fluids. 1967. V. 10, № 2. P. 354.
[10] АндриатисА.В., Соколова И.А. // Математическое моделирование. 1995. Т. 7. № 2. С. 89.
[11] Энгель А. Ионизованные газы. М.: Физматгиз, 1959. 332 с.
[12] Шейндлин А.Е., Асиновский Э.И., Батурин В.А., Батенин В.М. // Журнал технической физики. 1963. Т. XXXIII. № 10. С. 1169.
Calculation of Division Boards of Radiation and Heat Conduction Regimes of the Arc Discharge
Evgeniy N. Vasil’ev
Institute of Computational Modeling SB RAS 5Q/44 Akademgorodok, Krasnoyarsk, 66QQ36 Russia
By results of calculation in axisymmetric formulation of structure, specific and integrated capacities of energy mechanisms of the stationary equilibrium arc discharge division boards of the radiation and heat conduction regimes of the discharge burning in air, carbonic gas and argon are defined.
Keywords: the arc discharge, power balance, radiation, heat conduction.
