CC BY
72
S. A. Levchuk
MODELING OF SYMMETRICAL STRAINLY-DEFORM STATE OF TWO PLATE BODY COMPOUND JOINTED AT THE RIGHT ANGLE WITH GREEN’S TYPE MATRIX
Задача о статическом деформировании составного тела решается с использованием специально построенных матриц типа Грина. На примере составного тела из двух пластин при специальных краевых условиях продемонстрирована эффективность метода расчета.
Ключевые слова: статическое деформирование, матрицы типа Грина, составное тело из двух пластин, специальные краевые условия.
The problem of compound body static deformation is solved with use of spesially constructed Green’s type matrix. The calculation method effictiveness had been shown for the compound body with two plates with special border conditions.
Key words: static deformation, Green’s type matrix, two plate body compound, special border conditions.
УДК 519.63:533.9.07
Е. В. Цегельник, канд. физ.-мат. наук Е. К. Островский, В. О. Гарин Национальный аэрокосмический университет им. Н. Е. Жуковского «ХАИ», г. Харьков
ГАЗОДИНАМИЧЕСКИЙ ЗАТВОР АТМОСФЕРНОГО ВОЗДУХА В ИСТОЧНИКАХ НИЗКОТЕМПЕРАТУРНОЙ ПЛАЗМЫ -ПЛАЗМОТРОНАХ
Проведено математическое моделирование процесса молекулярной диффузии атмосферного воздуха во встречном потоке защитного инертного газа на участке газодинамического затвора в канале слива плазмы из плазмотрона. Показано, что при массовом расходе инертного газа аргон, вплоть до 10кг/с, течение в канале остается ламинарным. Установлено, что основным источником поступления отравляющего катод атмосферного воздуха в плазмотрон является его диффузия во встречном потоке инертного газа в сливном канале плазмы на участке газодинамического затвора атмосферного воздуха. Наиболее интенсивно молекулярная диффузия воздуха происходит в пограничном слое у стенок сливного канала плазмы, который на своем начальном участке выполняет функции газодинамического затвора.
Ключевые слова: низкотемпературная плазма, катод, термоэлектронная эмиссия, дуговой разряд, газовый поток, ламинарное течение, плазмотрон.
Введение
В промышленности высокоразвитых стран отмечается тенденция к расширению применения плазмотронов как источников низкотемпературной плазмы. Разработаны технологии получения чистых металлов и сплавов при прямом восстановлении их из руд и рудных концентратов, при подогреве в ходе внепечной обработки и разливе металлов [1, 2].
Особенностью плазмотронов, работающих в металлургии в условиях интенсивного испарения металлов с поверхности расплавов, является необходимость использования малых рабочих напряжений, исключающих появление неуправляемых пробоев. Такой режим
работы возможно реализовать только в низковольтном сильноточном дуговом разряде [3].
Большая мощность в таком разряде достигается за счет больших токов электронной эмиссии с развитой поверхности в полом катоде. Такие токи возможно реализовать переходом от электрической дуги с пятном на катоде к диффузному равномерно распределенному разряду по всей поверхности катода. Это возможно реализовать выбором как Фомы катода (полый катод), так и переходу к новому классу «горячих» высокотемпературных электродов - термоэмиссионным катодам с малой работой выхода электрона Ф = 2,0...2,4 эВ в диапазоне температур 1200 .1600 К [4].
© Е. В. Цегельник, Е. К. Островский, В. О. Гарин, 2010 120
Результаты исследования и их обсуждение
В исследуемом дуговом плазмотроне с выбросом низкотемпературной плазмы в атмосферу для защиты термоэмиссионного катода от отравления атмосферным воздухом используется газодинамический затвор, роль которого выполняет поток защитного и одновременно плазмообразующего инертного газа аргон (Аг), продуваемого через сливной канал плазмы плазмотрона со стороны катода. Основным источником поступления отравляющего катод атмосферного воздуха является его молекулярная диффузия во встречном потоке инертного газа на участке газодинамического затвора, расположенного в начале сливного канала плазмы. Г азодинамические характеристики встречного потока аргона определяют степень отравления катода атмосферным воздухом и его ресурс [5].
Конструкция плазмотрона в разрезе приведена на рис. 1.
Рис. 1. Разрез плазмотрона с показом линий тока в потоке инертного газа:
0-Ь - газодинамический затвор; Ь-К - полый катод
В канале плазмотрона изображены линии газового потока, указывающие на ламинарный характер течения потока газа на участке газодинамического затвора (0-Ь). Картина течения получена численным моделированием течения в программном продукте Со8то8¥1оч>отк8 [6]. Для цилиндрического канала ламинарный режим течения будет сохраняться при значениях числа Рейнольдса Яе < 2300. Число Рейнольдса Яе = 2М/жцЯ зависит от массового расхода газа М [кг/с], радиуса канала плазмотрона Я [м] и коэффициента вязкости Аг пАг [кг/См-с)], который зависит от температуры [7]. Расчеты температурной зависимости числа Яе для плазмотрона с Я = 0,02 м для всех возможных температурных режимах газового потока в канале плазмотрона Т = 300.20000 К, в диапазоне массовых расходов Аг М = 10-5...10-3 кг/с показывают, что число Яе не превышает своего критического значения 2300. Это указывает на то, что течение инертного газа в широком диапазоне рабочих параметров остается ламинарным.
Описание взаимной молекулярной диффузии атмосферного воздуха во встречном потоке Аг при его ламинарном течении описывается в общем случае нестационарным дифференциальным уравнением диффузии, полученного из уравнения баланса молекул воздуха в потоке инертного газа в сливном канале на участке газодинамического затвора [8]
дП2
дґ
= Б-
д 2 п
д2
- + Уг ^
2 2Г дг
(1)
где п2г - концентрация молекул (атомов) воздуха, как функция координат г иг ; У2г - скорость потока Аг, как функция координат г и г ; о - коэффициент взаимной молекулярной диффузии воздуха в аргоне.
Если распределение концентрации во времени установилось в плазмотроне дп2/ д/ = 0, то уравнение (1) преобразуется к виду
Откуда
^ ІЇП7Г
Б—— + У2Гп2Г = соші
(2)
Преобразуем запись дифференциального уравнения (2) к стандартному виду
П2Т + РП2Т Q,
(3)
где Р = У2г/О и сош1/О = Q .
Уравнение (3) - линейное дифференциальное уравнение первого порядка. Решение такого уравнения записывается в виде [8]
| &
,1
-| Рёг
(4)
Приближенное аналитическое решение уравнения стационарной диффузии возможно получить, если
функциональные зависимости Р{т,т) и &(гг) в уравнении (4) заменить их средними постоянными значениями Р и & на длине интегрирования 0-Ь в газодинамическом затворе.
После подстановки Р и & и интегрирования урав -нения (4) решение преобразуется
п = п ’ + Се ~Рг,
где п' = соші/У - концентрация молекул воздуха в инертном газе.
По ГОСТ 10157-79 на промышленный жидкий и газообразный инертный газ Аг концентрация примесей атмосферных газов п колеблется в пределах
п =
е
1021.1022 м-3 в зависимости от сорта Аг.
Определим константу интегрирования С в последнем уравнении из начальных условий. При г = 0, п = П0 = 2,5 -1025 м-3 - концентрация молекул воздуха при нормальных условиях атмосферы на выходе сливного канала плазмотрона, тогда С = П0 - п'. Учитывая, что П0 >> п’ и С = П0 уравнение (4) запишется в виде
У
п = п' + по ехр| --=г |.
(5)
Умножая на кТ уравнение (5) перейдем к записи решения через парциальное давление воздуха в среде аргон + воздух
Р = Р' + Р0 ехр
(6)
Эксперименты, как и теоретические расчеты, указывают на наличие значительного градиента температуры в поперечном сечении в потоке нагретого газа и плазмы в канале слива плазмы [9].
В пристеночной области охлаждаемого сливного
канала плазмы, где г/Я = 0,8..1, температура в потоке Аг 600.3000 К и диффузия атмосферного воздуха в сторону катода проходит в молекулярном состоянии.
В горячей зоне, вблизи оси потока, где г/Я = 0.0,8, температура газа 3000.25000 К. При такой температуре происходит полная диссоциация газовых молекул воздуха. В центре потока имеет место также концентрационная диффузия, но только атомов азота, кислорода и углекислого газа.
Входящий в уравнение (6) коэффициент взаимной диффузии записывается в виде
(кТ )2 (тАг + тв )
п тАг тв р(гАг + гв У
(7)
пЦв
Аг
где для недиссоциированного воздуха: тв = 49-10"27 кг -средняя масса молекулы воздуха, гв = 1,82-10"10 м -
средний радиус молекул воздуха, ц в = 0,029 кг/моль -молярная масса молекулярного воздуха; для диссоциированного воздуха: тв = 25-10-27 кг - средняя масса атома воздуха, гв = 1,2-10-10 м - средний радиус атома воздуха, ц в = 0,0145 кг/моль - молярная масса атомарного воздуха; тАГ = 67-10-27 кг - масса атома Аг, гаг= 1,83-10-10кг - радиус атома Аг цаг = 0,040 кг/моль -молярная масса Аг; Р = 1,01 -105 Па - давление в кана-
ле; к = 1,38-10-23 Дж/К - постоянная Больцмана; ^ = 8,314 Дж/(моль-К) - газовая постоянная; Т - температура газа.
Температурные зависимости коэффициента взаимной диффузии для диссоциированного атомарного и молекулярного воздуха при ламинарном течении газа представлены на рис. 2.
Рис. 2. Зависимость коэффициента взаимной диффузии атомного и молекулярного воздуха в аргоне оттемперату-ры при Р = 105 Па:
______атомарный кислород, • • • молекулярный кислород,
+ + + молекулярный воздух
Как следует из графиков коэффициенты взаимной диффузии для диссоциированного атомарного и молекулярного воздуха отличаются в 1,5 раза в рабочем диапазоне температуры газа 300.25000 К. Наибольшее изменение Б наблюдается в диапазоне температуры 300.5000 К, т.е. в пристеночной области сливного канала в плазмотроне.
В длинных каналах слива плазмы, длиной более 7.8 см, и небольших массовых расходах вязкого инертного газа, менее 10-3 кг/с, в канале устанавливается ламинарное течение. Распределение скоростей по радиусу в таком потоке описывается формулой Хагеля-Пуазейля
У =■
2М
.2 А
1 --
Я2
(8)
где М - массовый расход инертного газа, Я - радиус сливного канала на участке газодинамического затвора, р - плотность газа.
В установившемся ламинарном течении в сливном канале продольная скорость не зависит от продольной координаты г , что позволяет принять ее как постоянную У в уравнении (5).
+
Среднее уравнение коэффициента диффузии Б в том же уравнении (5) находится для каждого значения
отношений г/Я как среднее арифметическое значений Б на расстояниях г от торца газодинамического затвора в направлении к катоду плазмотрона.
Решение уравнения (6) с граничным условием по
воздуху Р0 = 1,01 -105 Па позволило рассчитать распределение парциального давления атмосферного воздуха в потоке инертного газа Аг в сечениях г = 0,6.. .0,9 м от торца газодинамического затвора в направлении к катоду плазмотрона.
В пристеночном пограничном слое, где г = (0,8.1)Я , расчет проведен для молекулярного воздуха. В ядре дуги, где г = (0.0,8) Я - для осреднен-ных параметров диссоциированного воздуха в атомном состоянии. Распределение парциального давления воздуха в поперечных сечениях сливного канала на участке газодинамического затвора, представленное на рис. 3, указывает на решающую роль продольной концентрационной диффузии воздуха у стенок газодинамического затвора в сливном канале плазмотрона на проникновение натекания его с торца плазмотрона в область катода.
1 0,8 0,6 0,4 0,2 0 0,2 0,4 0,6 0,8 1
г/П
О— г= 0,01 м; -о- т= 0,04 м; т= 0,05 м; —— т= 0,07 м; —*— т= 0,08 м; —С— г= 0,09 м
Рис. 3. Распределение парциального давления атмосферного воздуха в разных сечениях канала плазмотрона
Предельно допустимое парциальное давление атмосферного воздуха в область расположения термоэмиссионного катода, установлено в работе [5] и равно Рв ііт= 1-10-2 Па. Учитывая, что по основным ком -понентам воздуха и парам воды в составе плазмообразующего инертного газа Аг технической очистки парциальное давление равно Р = 0,6-10-2 Па, то допустимая диффузионная добавка не должна превышать Рвдоб = 0,4-10-2 Па.
Выбор длины канала газодинамического затвора г в плазмотроне с термоэмиссионным катодом определяется из условия, что парциальное давление диффундирующего отравляющего атмосферного воздуха не должно превышать предельного значения Рпр = 0,4-10-2 Па на расстоянии г в месте расположения катода в плазмотроне.
Из результатов расчета (рис. 3) следует, что эффективно защищен от отравления лишь центр катода, где скорость встречного потока инертного газа максимальна. Хуже защищены участки катода по периферии при г/Я > 0,8. Эти участки катода попадают в пристеночный пограничный слой в потоке защитного инертного газа, где скорость течения, согласно уравнения (8), мала. Для усиления защиты катода по периферии длина газодинамического затвора должна быть увеличена до значения г > 0,9 м.
Заключение
1. Методами математического и компьютерного моделирования уравнений газодинамики исследовано течение плазмообразующего инертного газа аргон в плазмотроне
2. Установлен ламинарный характер течения аргона на участке газодинамического затвора вплоть до расхода аргона М = 10-3 кг/с.
3. Изучен вклад продольной взаимной диффузии отравляющих катод атмосферных газов через встречный ламинарный поток плазмообразующего инертного газа аргон в общее парциальное давление атмосферы в области катода.
4. Проведенные расчеты подтверждают, что основной диффузионный поток атмосферы в сторону термоэмиссионного катода формируется у стенок сливного канала плазмотрона на участке газодинамического затвора, где скорость потока аргона Угг = 0.
Перечень ссылок
1. Мельник Г. А. Перспективы использования плазменных источников теплоты в агрегатах внепечной обработки стали / Г. А. Мельник, О. С. Забарило,
A. А. Ждановский и др. // Проблемы специальной электрометаллургии. - 1991. - Вып. 2. - С. 60-66.
2. Жадкевич М. Л. Создание и освоение плазменного ковша-печи / М. Л. Жадкевич, Г. А. Мельник, О. С. Забарило и др. // Проблемы специальной электрометаллургии.- 1998.- Вып. 1. - С. 42-49.
3. Плазменная металлургия / Цветков Ю. В., Николаев А. В., Панфилов С. А. и др. - Новосибирск : Наука. Сиб. отд-ние, 1992. - 265 с. - (Низкотемпературная плазма. Т. 8).
4. Кривцов В. С. Проблемы создания высокоресурсных сильноточных электродуговых плазмотронов /
B. С. Кривцов, С. И. Планковский // Авиационнокосмическая техника и технология. - 2005. - № 7 (23). -
C. 7-21.
5. Кривцов В.С. Состав газа в плазмотроне с ламинарным потоком плазмообразующего инертного газа / В. С. Кривцов, С. И. Планковский, Е. В. Цегельник, Е. К. Островский, А. А. Таран, А. В. Лоян // Авиационно-космическая техника и технология. - 2005. - № 4 (20). - С. 19-24.
6. Планковский С. И. Методика расчета основных газодинамических параметров потока инертного газа в плазмотроне / С. И. Планковский, Е. В. Цегельник, Е. К. Островский, В. А. Максименко // Авиационно-космическая техника и технология. - 2007. - № 7 (43). - С. 94-98.
7. Ваграфтик Н. Б. Справочник по теплофизическим свойствам газов и жидкостей / Н. Б. Варгафтик. - М. : Наука, 1972. - 721 с.
8. Выгодский М. Я. Справочник по высшей математикера / М. Я. Выгодский. - М. : Наука, 1977. - 872 с.
9. Жуков М. Ф. Электродуговые нагреватели газа (плазмотроны) / М. Ф. Жуков, В. Я. Смоляков, Б. А.Урюков. -М. : Наука, 1973. - 232 с.
Одержано 01.10.2009
E. V. Tsegelnyk, E. K. Ostrovsky, V. O. Garin A GAS-DYNAMIC SHUTTER OF ATMOSPHERIC AIR IN LOW TEMPERATURE PLASMA SOURCES - PLASMOTRONS
Проведено математичне моделювання процесу молекулярної дифузії атмосферного повітря в зустрічному потоці захисного інертного газу на ділянці газодинамічного затвору в каналі викиду плазми із плазмотрона. Показано, що при масовій витраті інертного газу аргон, аж до 10-3 кг/с, течія у каналі залишається ламінарною. Встановлено, що основним джерелом надходження атмосферного повітря в плазмотрон, яке отруює катод, є його дифузія в зустрічному потоці інертного газу в зливному каналі плазми на ділянці газодинамічного затвору атмосферного повітря. Найбільш інтенсивно молекулярна дифузія повітря відбувається в прикордонному шарі біля стінок зливного каналу плазми, що на своїй початковій ділянці виконує функції газодинамічного затвору.
Ключові слова: низькотемператрна плазма, катод, термоелектронна емісія, дуговий розряд, газовий потік, ламінарна течія, плазмотрон.
The mathematical modeling of atmospheric air’s molecular diffusion process in oncoming rare gas stream on the gas-dynamic lock region of the plasmatron is carried out. It is shown, that during the rare gas’s argon mass-consump-tion up to 10-3 kg/s, the laminar flow in channel remains. It was determined, that the main source of an atmospheric air’s entry in the plasmatron is the air’s diffusion in oncoming rare gas stream on the gas-dynamic lock region. The molecular air diffusion takes place with the most intensity in the boundary layer near the walls of the confluent plasma’s channel, which is the above mentioned gas-dynamical lock.
Key words: low temperature plasma, cathode, termo-electrical emission, arc digit, gas stream, laminar flow, plas-motron.
УДК 621.793
П. И. Можаев, канд. техн. наук Е. Я. Швец, д-р техн. наук В. Н. Михайлин,
Л. Л. Веревкин, канд. техн. наук Н. В. Светанько
Государственная инженерная академия, г. Запорожье
МОДЕЛЬ РАСЧЕТА ЭФФЕКТИВНОСТИ ЭКРАНИРОВАНИЯ ЖЕЛЕЗОЦИНКОВЫХ ФАЗ В ДИФФУЗИОННЫХ ЦИНКОВЫХ ПОКРЫТИЯХ
Предложен новый метод расчета эффективности экранирования железоцинковых фаз в термодиффузионных цинковых покрытиях от электромагнитного излучения.
Ключевые слова: железоцинковые фазы, диффузия, ферромагнетик, диамагнетик, бинарная система, гомогенность, эвтектоидное превращение.
Анализ достижений , „
вых фаз и слоя почти чистого цинка. При этом строе-
Диффузионное цинковое покрытие, полученное в ние покрытия повторяется даже при самой кратковре-
расплаве цинка, состоит из нескольких железоцинко- менной выдержке стали в расплаве цинка.
© П. И. Можаев, Е. Я. Швец, В. Н. Михайлин, Л. Л. Веревкин, Н. В. Светанько, 2010
124
