CC BY
12
равномерно сходится к решению уравнения (7), причем
К (х) - <2 (х)|ф! ,Ь] - + , 0 < ^ = СОШ*■
Следствие 2. При выполнении условий теоремы 2 решение уравнения (7) единственно в пространстве СГ\ЪХ, Ь], 0 </ < 1.
Доказательство теоремы 2 проведено в работе [2, с. 49-51] и используются аналогичные выкладки как в доказательстве теоремы 1.
Вышеизложенные теоремы и следствия позволяют заключить, что имеет место равномерная сходимость функции срЕ(х) к решению р(х) уравнения (1) для всех
х е [0,Ь], при £ ^ 0, где функции рЕ(х) и р(х) определены по правилу
р (х) = (х), х е [0, Ь1 ф) = <1(х), х е [0, Ь11
£ \(х), х е [Ь1,Ь] < (x), х е [ьl, Ь\
причем выполняются условия согласования \(Ь) = р2г(Ь), <р(Ь) = Рг(Ь)■
Литература
1. Каракеев Т. Т., Рустамова Д. К. Регуляризация и метод квадратур для линейных интегральных уравнений Вольтерра третьего рода // Исслед. по интегро-дифференц. уравнениям. - Бишкек: Илим, 2009 - Вып. 40. - Стр. 127-132.
2. Омуров Т. Д., Каракеев Т. Т. Регуляризация и численные методы решения обратных и нелокальных краевых задач. - Бишкек, «Илим» 2006 г. - 164 с.
3. Демидович Б. П. Лекции по математической теории устойчивости. - Москва: Наука, 1967. - 472 с.
Расчет электрической дуги в канале в аксиальном магнитном поле Урусова И. Р.1, Урусова Т. Э.2
'Урусова Индира Руслановна / Urusova Indira Ruslanovna - кандидат физико-математических
наук, старший научный сотрудник; 2Урусова Толкун Эсеновна / Urusova Tolkun Esenovna - кандидат физико-математических наук,
ведущий научный сотрудник, Институт физико-технических проблем и материаловедения, Национальная академия наук Кыргызской Республики, г. Бишкек, Кыргызская Республика
Аннотация: в рамках трехмерной нестационарной математической модели в приближении частичного локального термодинамического равновесия плазмы выполнен расчет электрической дуги постоянного тока, горящей в канале во внешнем аксиальном магнитном поле. По результатам расчета получена винтовая форма столба дуги.
Ключевые слова: электрическая дуга, численное моделирование, трехмерная нестационарная математическая модель, внешнее магнитное поле, винтовая форма дуги.
Введение. Известно, что цилиндрическая осевая симметрия протяженной электрической дуги при наложении внешнего аксиального магнитного поля (ВАМП) может нарушаться, и дуга принимает винтовую пространственную форму [1-3]. Принято считать, что причиной перехода столба дуги от цилиндрической формы к винтовой являются малые флуктуации параметров дуговой плазмы.
Экспериментальные и теоретические исследования винтовой формы дуги сравнительно немногочисленны, практически отсутствуют численные исследования в рамках трехмерных математических моделей. В связи с этим изучение физических процессов формирования винтовой формы дуги представляет определенный научный интерес.
Постановка задачи и математическая модель. В цилиндрическом канале длиной Ь, диаметром Б в аргоне атмосферного давления горит электрическая дуга постоянного тока с межэлектродным расстоянием I. Условная схема расчетной области показана на рис.1. Катодом «-» и анодом «+» являются графитовые цилиндрические стержни с плоской токоведущей поверхностью. Внешнее аксиальное однородное магнитное поле Н совпадает с направлением оси х декартовой системы координат х, у, г.
Нестационарная трехмерная система уравнений в приближении частичного локального термодинамического равновесия (ЧЛТР) плазмы может быть записана в следующем виде [4-6]:
Рис. 1. Условная схема расчетной области дуги в канале во внешнем аксиальном магнитном поле Н
уравнение непрерывности газа тяжелых частиц:
др/дг + ¿МрУ) = 0 (1) уравнение непрерывности газа электронов:
дые/дг + ¿МЖеУе)=Яе (2) уравнение сохранения энергии газа электронов:
д/дг(Ъ/2кТе+ЩМе+$м[(5/2кТе+ЩИеУе] = ¿гу(Я^гайТе)+|2/ст-^-£(Те -Т) (3) уравнение сохранения энергии газа тяжелых частиц:
д/дг[3/2кТ(М+М,)] +div[5/2kT(NI+Na)У)]=div(ЯgradT)+£(Te -Т) (4) уравнение баланса импульса газа вдоль осей координат х, у, г:
дри/дг + ¿гу(рУи) = ¿гу(цгайи) - дР/дх + цо(|хН)х + 5х+(р-рсо)§ (5) дрм.де + вшм(рМм) = вшм(цпкфвм) - дЗ.дн +ц0(охР)н + ын (6) дрц.де + вшм(рМц) = вшм(цпкфвц) - дЗ.дя+ц0(охР)я + ы (7) уравнения Максвелла:
кще У = Об кще Р = об вшм Р = Об (8) закон Ома в обобщенной форме:
цо (МхР) + У = о.ст+ (ц охР-пкфв Зу).йуТу (9) закон парциальных давлений:
Р/кТ = N +Ыа + ЫеТД. (10)
При записи системы уравнений приняты обозначения: t - время, р - плотность дуговой плазмы, рт - плотность окружающего холодного газа, 1е - коэффициент теплопроводности газа электронов, 1 - коэффициент теплопроводности газа, ц -коэффициент вязкости, ст - коэффициент электропроводности, у - коэффициент излучения, Ыа, Ые - концентрация ионов, атомов и электронов, Яе = (К Ые Ыа - Кг
Щ - скорость генерации электронов, К - константа ударной ионизации, Кг -константа трехчастичной рекомбинации, - ионизационный потенциал
плазмообразующего газа, Рг = ЫекТе - парциальное давление электронного газа, к -постоянная Больцмана, В - коэффициент энергообмена между электронами и тяжелыми частицами (атомы, ионы) при соударениях, g - ускорение свободного падения, де - элементарный электрический заряд (электрона), ц0 - магнитная константа, V - вектор скорости газа, Е, Н, ] - соответственно векторы напряженности электрического поля, собственного магнитного поля и плотности тока, Т -температура газа тяжелых частиц, Тг - температура электронного газа, Р - давление, и, V, V - компоненты вектора скорости V в направлении осей х, у, I, Vd = ]/(де№е) -вектор скорости дрейфа электронов, V = - -Оа/Те^а^е - вектор скорости термодиффузии, V., = - Вс/Ы&хадЫе - вектор скорости амбиполярной диффузии, Ба -коэффициент амбиполярной диффузии, V,, = V+Vd+Vt+Va суммарная скорость электронов, sz, sy, sx - вязкие слагаемые.
Принято, что дуговая плазма является однократно ионизованной, квазинейтральной, течение ламинарное, дозвуковое, излучение объемное; вязкой диссипацией энергии, индукционными токами пренебрегается [4]. Приэлектродные процессы не рассматриваются. Коэффициенты переноса и свойства плазмы аргона являются функциями температуры электронов и тяжелых частиц и рассчитываются в соответствии с известной методикой [4].
Электромагнитная часть задачи решается с использованием скалярного потенциала электрического поля р и векторного магнитного потенциала А. Используя известные соотношения Е = - gradр, го!А=Н, закон Ома, закон неразрывности электрического тока divj = 0 и уравнения Максвелла, получим уравнения для расчета скалярного потенциалар и компонент векторного потенциала Аг, Лу, Ах, которые имеют следующий вид:
div(стgradр)=div[стцo (УхН) -ст (цojxH-gradPe)/qeNe], (11) div(gradЛx)= - ]х, div(gradЛy)= - ¡у, div(gradAг)= - (12)
Отметим, что в работе [7] предложена математическая модель, которая позволяет описывать процессы в канале МГД-устройства с коническим осесимметричным каналом, не прибегая при этом к решению сложных дифференциальных уравнений.
Исходная система уравнений (1-12) для рассчитываемых переменных после несложных преобразований может быть записана согласно известной методике [8] в виде обобщенного дифференциального уравнения:
да рФ/дt+div(PрVФ) = div(ygradФ) + д (13)
где Ф - одна из переменных: Аг, Ау, Ах, р, V, V, и, Т а , Т, Те, Ые, коэффициенты а, Д
у, д зависят от смысла переменной Ф.
Дискретизация нестационарного обобщенного дифференциального уравнения (13) проводится по неявной разностной схеме методом контрольного объема [8], численное решение конечно-разностного аналога проводится методом Зейделя-Гаусса с применением нижней релаксации. Используется метод фиктивных областей [9], адаптированный для расчета характеристик электрической дуги [10].
Граничные и начальные условия. Во входном и выходном круговых сечениях канала (см. рис.1) для расчетных характеристик дуги задаются условия дФ/дх = 0 гладкого сопряжения с внешней средой. Потенциал электрического поля р
рассчитывается из условия протекания тока по нормали к токоведущим торцевым поверхностям электродов.
Значения температуры и концентрации электронов равны значениям температуры Тетт и концентрации Жетт «холодного» не ионизованного газа: Те = Тетт = 3 кК, Ие = Жетт = 1017 м-3.
На боковой поверхности канала течение электрического тока отсутствует, скорость газа равна нулю, температура полагается равной 500 К. Температура и концентрации электронов равны Те = Тет1П, N = Жетт
При постановке начальных условий при г = 0 полагается, что между электродами существует токопроводящая высокотемпературная (Т = 10 кК) зона в форме цилиндра с неподвижным газом.
Обсуждение результатов. Выполнен расчет дуги силой тока I = 40 А, длиной I = 10 см в канале длиной Ь = 12 см и диаметром В = 10 мм. Величина ВАМП равна Н = 6 кА/м. Сеточный шаг разностной задачи в направлениях осей х, у, г одинаковый Д = 0,5 мм, временной шаг т = 10-4 с.
На рис. 2а приведена эволюция расчетных распределений температуры дуги в вертикальном среднем сечении Х-У при г/2 (качественно такие же расчетные распределения наблюдаются в вертикальном среднем сечении Х-2 при у/2). Укажем, что в течение первых 100 мс после инициировании разряда, численный расчет характеристик дуги проводится без ВАМП.
В дальнейшем, при наложении ВАМП, на интервале времени примерно 100 -г- 280 мс наблюдается переход от цилиндрической формы столба дуги к винтовой форме. С момента времени г > 300 мс изменения характеристик и формы столба дуги не значительны, и к моменту г = 400 мс численный расчет был остановлен. Обращает на себя внимание, что винтовые возмущения зарождаются вблизи анода «+» и распространяются с течением времени в сторону катода «-». Причины этого требуют дополнительных исследований. Следует сказать, что физика процессов, обусловливающих винтовую форму столба дуги, является достаточно сложной и до конца не изученной [1- 3].
В канале наблюдается система вихрей, скалярное поле скорости газа,
определяемое по формуле V = л/и2 + V2 , показано на рис. 3б. Отметим, что в отсутствие осевой симметрии протекающих процессов характер течения дуговой плазмы является достаточно сложным, и без трехмерной графики представить общую картину течения весьма затруднительно.
Рис. 2. Распределения температуры газа Т и скалярного поля скорости газа V в различные моменты времени г. I = 40 А, I = 10 см, Н = 6 кА/м.
Заключение. В рамках нестационарной трехмерной математической модели в приближении частичного локального термодинамического равновесия плазмы численно реализована винтовая форма электрической дуги во внешнем аксиальном магнитном поле, горящая в цилиндрическом канале. Результаты численного расчета удовлетворительно согласуются с результатами эксперимента [3]. Предложенная математическая модель может быть использована в исследовании причин, обусловливающих переход от цилиндрической формы столба дуги к винтовой форме при наличии ВАМП.
Литература
1. Maecker H. Principles of arc motion and displacement, Proc. IEEE, 1971, Vol. 59, No. 4, P. 439-449.
2. Ментель Ю. Магнитная неустойчивость электрической дуги // Теория электрической дуги в условиях вынужденного теплообмена. Новосибирск. Наука, Сибирское отделение, 1977. С. 182-204.
3. Асиновский Э. И., Кузьмин А. К., Пахомов Е. П. Измерение геометрических параметров винтовой дуги // М.: Теплофизика высоких температур, 1980. т. 18. № 1. С. 9-15.
4. Энгельшт В. С., Гурович В. Ц., Десятков Г. А. и др. Низкотемпературная плазма, т. 1. Теория столба электрической дуги. Новосибирск: Наука, 1990. 374 с.
5. Чередниченко В. С., Аньшаков А. С., Кузьмин М. Г. Плазменные электротехнологические установки. Новосибирск: НГТУ, 2005. 508 с.
6. Урусов Р. М., Урусова И. Р. Нестационарная трехмерная модель электрической дуги, ч. 1. Математическая модель и результаты тестирования // Теплофизика и аэромеханика, 2014. т. 21, № 1. С. 121-134.
7. Хайруллин И. Х., Камалов Ф. А. Математическое моделирование процессов в канале МГД-устройства с коническим осесимметричным каналом // Инженерный вестник Дона, 2012, № 4, [Электронный ресурс]. Режим доступа: URL: ivdon.ru/ru/magazine/archive/n4p2y2012/1444.
8. Patankar S. V. Numerical Heat Transfer and Fluid Flow, Hemisphere Publ. Corp., New York, 1980.
9. Смагулов Ш., Сироченко В. П., Орунханов М. К. Численное исследование течений жидкости в нерегулярных областях. Алматы, 2001. 276 с.
10. Урусов Р. М., Урусова Т. Э. Применение метода фиктивных областей для расчета характеристик электрической дуги // М.: Теплофизика высоких температур, 2004. т. 42. № 3. С. 374-382.
Численное моделирование процесса гашения электрической дуги
во внешнем поперечном магнитном поле 1 2 Урусова Т. Э. , Урусова И. Р.
1Урусова Толкун Эсеновна / Urusova Tolkun Esenovna - кандидат физико-математических наук,
ведущий научный сотрудник; 2Урусова Индира Руслановна / Urusova Indira Ruslanovna - кандидат физико-математических
наук, старший научный сотрудник, Институт физико-технических проблем и материаловедения, Национальная академия наук Кыргызской Республики, г. Бишкек, Кыргызская Республика
Аннотация: в рамках трехмерной нестационарной математической модели выполнен расчет электрической дуги во внешнем поперечном магнитном поле. Установлено, что при некотором критическом значении внешнего магнитного поля наблюдается разрыв столба дуги и ее гашение. Результаты качественно согласуются с опытными данными.
Ключевые слова: электрическая дуга, численное моделирование, трехмерная нестационарная модель, внешнее поперечное магнитное поле, гашение дуги.
Введение. Низкотемпературная плазма, в том числе электродуговая, находит широкое применение в различных технологических аппаратах и процессах [1-3], и в целях оптимизации режимов работы необходимы дальнейшие исследования. Так, в [4] рассмотрены проблемы построения фрактальной модели магнитоплазменного электродинамического ускорителя. В настоящей статье представлены результаты моделирования гашения электрической дуги во внешнем поперечном однородном магнитном поле.
Постановка задачи и математическая модель. В декартовой системе координат х, y, z рассчитывается открытая электрическая дуга постоянного тока, горящая в аргоне атмосферного давления.
