CC BY
4Электронный журнал «Труды МАИ». Выпуск № 44
www.mai.ru/science/trudv/
УДК 533
Процессы переноса вблизи цилиндрического зонда в потоке плотной плазмы в магнитном поле.
Котельников В.А., Котельников М.В., Мальцев-Горский Д.А.
Аннотация: Сформулирована математическая и численная модели процессов переноса электрического заряда вблизи цилиндрического зонда в поперечном потоке плотной плазмы при наличии осевого магнитного поля.
Обнаружена нелинейность в поведении средней плотности электронного тока в зависимости от параметра Холла. В ограниченном интервале изменения числа Рейнольдса кривые проходят через максимум, который, который сохраняется при изменении безразмерного радиуса зонда, его потенциала из отношения температур ионов и электронов. Дано физическое объяснение обнаруженному нелинейному эффекту.
Ключевые слова: процессы переноса, уравнение эйлера, уравнение пуассона, метод крупных частиц, магнитное поле, плотная слабоионизованная плазма, параметр холла.
Введение
Исследования процессов переноса в потоках плотной плазмы в магнитном поле актуально для развития методов зондовой диагностики, для понимания процессов в МГД-генераторах, некоторых типов плазмотронов, плазмохимических реакторов и других устройств, использующих потоки плотной плазмы в магнитном поле.
В настоящей работе основное внимание уделяется переносу заряда на положительно заряженный цилиндрический электрод, находящийся в поперечном потоке слабоионизованной плотной плазмы в продольном магнитном поле (рис. 1)
Радиус цилиндра г , его потенциал р . Предполагается, что длина зонда много больше его радиуса.
Постановка задачи
Математическая модель задачи включает уравнения неразрывности для ионов и электронов, уравнения движения заряженных компонент, уравнение Пуассона для самосогласованного электрического поля, а так же систему уравнений Эйлера для нейтральной компоненты. [1]
Рис. 1 Расположение зонда в потоке.
3п
а
3п
- + —IV (п и ) = 0
Ьё^Киe) = 0
а
с1и кТ
т =--L Уп1 + ге(Е + и х В) - V;а (и; - иа )
т„
1 -1 —и.
п
—1
кТ
п
Упе - е(Е + ие Х В) - Цеа Vea (ие - иа )
Аф = е(пе -2п1)/80,Е = -Уф
^ + —IV (раПа )= 0
3(Раиа )
3(РаЕа )
+ (Ра иа иа )== -УР
+ (РаиаЕа ) == (иаРа )
Ра =(Еа - и2/2)Ра (У - 1), Еа = CvTa + ^ / 2
(1)
В приведенной системе уравнений п, и, Т, р, л, о, Р, р, Е, 5 - концентрация, скорость, температура, плотность, приведенная масса, частота столкновений, давление, потенциал, напряженность электрического поля и индукция магнитного поля соответственно. Остальные обозначений общепринятые. Индексы г, е, а относятся к ионам, электронам и нейтральным частицам. Кроме системы (1) математическая модель включает также систему начальных и граничных условий. На границе втекания задавались параметры плазмы, на границе вытекания ставились условия свободного вытекания. Потенциал цилиндра считался постоянным, потенциал на удаленных границах полагался равным нулю. На электроде ставится условие непротекания для нейтральных частиц и условие идеальной каталитичности для заряженных компонент плазмы.
В системе (1) отсутствуют уравнения энергии для ионов и электронов. Это связано с тем, что в условиях слабой степени ионизации температура ионов полагается равной температуре нейтральных частиц, а температура электронов Т варьировалась в виде
Т
параметра е = —.
Те
Учтены также следующие допущения:
- химические реакции заморожены;
- вязкостью пренебрегаем;
- собственное магнитное поле мало в сравнении с внешним;
Сформулированная задача оказывается многопараметрической. Среди основных параметров укажем следующие:
Г =— (гв - радиус Дебая) - безразмерный радиус цилиндра;
р0 = - безразмерный потенциал цилиндра;
кТ
Р 1
• Р = —^ = -,= Р 42
гт\
V те J
гТ \
_е
Т
V Тг J
- отношение коэффициентов диффузии электронов и
ионов;
• и0 = 7—-—ч - безразмерная направленная скорость;
р 1ГР )
Т
• в = — - отношение температур ионов и электронов;
Те
гр
2
2ю1е еВ Р1е = —— - параметр Холла; (ще=--ларморовская частота вращения)
г,е г г 'V г ,е
и<„.„ т
г,е
Сформулированная задача решалась численно. Условие слабой степени ионизации позволяет отдельно решать уравнения для нейтральных частиц. Решение осуществлялось с использованием метода крупных частиц Давыдова [2,3]. Полученые поля скоростей, концентраций, температур нейтральных атомов использовалось как фон для решения электродинамической части задачи. Уравнения неразрывности и движения для ионов и электронов также решались методом крупных частиц. При этом на каждом временном слое напряженность электрического поля Е получалась путем решения уравнения Пуассона. Последнее решалось с использованием метода разделения переменных Фурье [4].
Результаты вычислительных экспериментов
Результаты вычислительных экспериментов представлены на рис. 2 ^ 11. Экспериментальные данные удобно представлять в зависимости от безразмерных комплексов, таких как число Маха М, равное отношению направленной скорости и0 к
скорости звука и электрическое число Рейнольдса Яеэ, равное отношению конвективных членов уравнения к диффузионным членам. Оно может быть выражено произведением ионного числа Шмидта на число Рейнольдса потока ( Яеэ = Яе- 8е1).
В случае положительно заряженного цилиндра > 0) на его поверхность
осуществляется перенос из окружающей плазмы отрицательного заряда, т.е., в зондовой цепи протекает электронный ток. Наличие направленной скорости приводит к появлению с теневой стороны цилиндра вытянутой возмущенной области (следа). Вследствие этого наряду с радиальным электрическим полем Е возникает также азимутальное электрическое поле . При отсутствии магнитного поля (В = 0) имеет место плоскость симметрии, проходящая через вектор скорости и и ось цилиндра. При включенном магнитном поле (В > 0) плоскость симметрии отсутствует. Это связано с тем, что в скрещенных полях Ев и
В возникает дрейфовое движение, которое направлено с одной стороны цилиндра от него, а с другой - по направлению к цилиндру. В теневой области в зависимости от параметров задачи может возникать вихревое движение, приводящее к возникновению локальных максимумов электронного тока.
На рис.2 приведена зависимость безразмерной плотности электронного тока ]е на поверхность цилиндра от двух параметров: параметра Холла Д и электрического числа Рейнольдса Яе3. Значения остальных параметров указаны на рисунке.
(Л ср = , м' = —)
р MJ гв
Рис. 2. Зависимость средней плотности электронного тока от параметров ße и Re3
Г = 10, M = 0.05 - 0.6, <р0= 5,s = 0,3)
1 - Re3=10; 2 - 25; 3 - 50; 4 - 100; 5 - 200; 6 - 500;7 - 1000.
Кривые, приведенные на рис. 2, в определенном интервале изменения параметров ßi и Re3 существенно нелинейны. При Re3 < 100 средняя плотность электронного тока j с ростом магнитного поля (оно входит в параметр ß) уменьшается, что согласуется с известными в литературе данными [5]: уменьшением поперечной диффузии заряженных частиц в магнитном поле. Однако в области Re3 > 100 в интервале ß < 10 2 jcp растет с
ростом индукции B, а при ß > 2 -10 2 начинает уменьшаться, как и при Re3 < 100.
Обнаруженный эффект проявляется при различных значениях числа Маха в интервале 0,05 <M < 0,6. На рис. 3-5 приведены аналогичные рис.2 зависимости, но при других
значениях параметра r0 (3 < r0 < 100)
)е ср
4,00 3,50 3,00 2,50 2,00 1,50 1,00 0,50
0,00
/
/
/
~ - - . - - _ _
0,00 0,01 0,02 0,03 0,04 0,05
-Р!е=25--(Чей 00-------!?е=200 ---Р!е=500
Рис.3 Зависимость средней плотности электронного тока от параметров Р; и
Яеэ (г0 = 3, М = 0.6, % = 5, е = 0,3)
Рис.4 Зависимость средней плотности электронного тока от параметров Р; и
Яеэ (г0 = 30, М = 0.6, % = 5, е = 0,3)
Рис. 5 Зависимость средней плотности электронного тока от параметров ре и Яеэ
(г0 = 100, М = 0.6, % = 5, е = 0,3).
Исследовалось также влияние на среднюю плотность электронного тока параметра е (0,1 <е< 2). Соответствующие кривые приведены на рис. 6,7.
Рис.6 Зависимость средней плотности электронного тока от параметров ßi и s (Re3 = 500, r0 = 3, M = 0.6, % = 5)
Рис.7 Зависимость средней плотности электронного тока от параметров ßi и s (Re3 = 500, r = 30, M = 0.6, % = 5)
На рис.8 ^ 10 даны аналогичные зависимости от параметра % (l <%< 7)
Je ср
Je ср
40,00
35,00 ■ / _ _ _ " N N
30,00 ■ 25,00 ■ ///> л--'// ///У
20,00 ■ //// 7/
15,00 ■ ff
10,00 ■
5,00 ■
0,00 ■
0,00 0,01 0,02 0,03 0,04 0,05
-tp = 1--ср = 3-------1(1 = 5---ip = 7
4,00 3,50 3,00 2,50 2,00 1,50 1,00 0,50
0,00
^ - -
^ sJ __ -
////
7//
7/
7
0,00 0,01 0,02 0,03 0,04 0,05
-q> = 1--ф = 3......■ tp = 5---(р = 7
Рис.8 Зависимость средней плотности электронного тока от параметров ßi и % (Re3 = 500, r = 3, M = 0.6, s = 0,3)
Рис.9 Зависимость средней плотности электронного тока от параметров ßi и % (Re3 = 500, r = 30, M = 0.6, s = 0,3)
Рис.10 Зависимость средней плотности электронного тока от параметров ßi и
(Re3 = 500, r0 = 100, M = 0.6, £ = 0,3)
Из приведенных на рис.2 ^ 10 кривых следует:
1. Обнаруженная в численных экспериментах нелинейность в поведении средней плотности электронного тока в зависимости от чисел Re3 и ß (в ограниченном интервале
их изменения) имеет место при различных значениях характерных параметров задачи
M, ^ ^
2. Зависимость j от числа Маха в интервале 0,05 <M < 0,6 относительно невелика. С ростом числа M j ср также растет, что связано с зависимостью Ев от и0.
3. Зависимость jcp от r0 качественно повторяется. С уменьшением r0 jecp
увеличивается, что совпадает с данными Лафрамбуаза [6] и других авторов [7,8]. Это вызвано увеличением пристеночной напряженности электрического поля при уменьшении радиуса кривизны зонда.
4. Зависимость jecp от обратна зависимости от r0. С ростом jecp растет,
поскольку возрастает пристеночная напряженность электрического поля.
T
5. С ростом электронной температуры (параметр s = —- уменьшается) максимум jecp
Te
смещается в сторону более высоких значений ß. Это объясняется тем, что более высокотемпературные электроны замагничиваются при более высоких значениях индукции магнитного поля B.
В заключение дадим физическое объяснение обнаруженной нелинейности в поведении кривых ] = /(д). С этой целью рассмотрим распределение плотности
электронного тока ^ по обводу цилиндра. При относительно больших числах Рейнольдса
(Яеэ > 100) эффекты, связанные с конвекцией, преобладают над диффузией и
подвижностью. Вследствие этого плотность тока на боковую поверхность цилиндра мала (заряженные частицы, участвующие в конвективном движении, пролетают мимо цилиндра). Имеет место только ток электронов на лобовую часть и на теневую, если там имеются вихри. При включении осевого магнитного поля возникающие пондеромоторные силы приводят к существенному росту тока ^ на одну из боковых поверхностей цилиндра, которое
Д = 0
Д = 10 -
Д = 2,5 -10-
Д = 10 1
Рис. 11. Влияние магнитного поля на поле скоростей электронов. (Яеэ = 104; М = 0,6; го = 10, Фс = 5)
2
2
превышает уменьшение тока на лобовую и теневую области и противоположный бок цилиндра. При дальнейшем увеличении индукции магнитного поля В (д > 10 2) растет степень замагниченности электронов и это ведет к уменьшению у также, как это имеет место при (Яеэ < 100). В качестве подтверждения приведенного физического анализа на рис.11 приведены (по результатам численных экспериментов) линии тока электронов вблизи поверхности цилиндра при Яеэ = 104 и четырех значениях параметра Холла
Д = 0; 10~2; 2,5 -10-2; 10 . При Д = 0 линии тока почти параллельны боковым поверхностям цилиндра. При Д = 10 2 под действием пондеромоторных сил с одной стороны цилиндра появляется составляющая вектора скорости электронов, направленная к поверхности цилиндра. При Д = 2,5 -10 2 эта составляющая достигает максимума, а при Д = 10 1 она снова уменьшается за счет замагничивания электронов.
Библиографический список
1. В.А. Котельников, М.В. Котельников, В.Ю. Гидаспов. Математическое моделирования обтекания тел потоками столкновительной и бесстолкновительной плазмы -М.:Физматлит, 2010, 266с.
2. О.М. Белоцерковский, Ю.М. Давыдов. Метод крупных частиц в газовой динамике. Вычислительный эксперимент. - М.: Наука, 1982, 392с.
3. Ю.М. Давыдов. Метод крупных частиц // Математическая энциклопедия. - М.: Сов. экциклопедия, 1985 - т.З, с.125-129.
4. Н.С. Бахвалев, Н.П. Жидков, Г.Н. Кобельков. Численные методы - М.: Наука, 1989, 608с.
5. Д.А. Франк-Каменецкий. Лекции по физике плазмы. - М.: Атомиздат, 1968, 286с.
6. J.G. Laframboise. Theory of cylindrical and spherical Langmuir probes in collisionless plasma at rest // Rarefied Gas Dunamics - 1966 - V.II, №4, Р.22.
7. Чан П., Телбот Л., Турян К. Электрические зонды в неподвижной и движущейся плазме (теория и приложения) - М.: Мир, 1976 - 202с.
8. Котельников В.А., Ульданов С.В., Котельников М.В. Процессы переноса в пристеночных слоях плазмы. - М.: Наука, 2004. - 422с.
Сведения об авторах
Котельников Вадим Алексеевич, профессор Московского авиационного института (государственного технического университета), д.т.н., МАИ, Волоколамское ш., 4, Москва, А-80, ГСП-3, 125993; тел.: (499) 158-19-70; e-mail: mvk_home@mail.ru
Котельников Михаил Вадимович, профессор Московского авиационного института (государственного технического университета), д.т.н., МАИ, Волоколамское ш., 4, Москва, А-80, ГСП-3, 125993; тел.: (499) 158-19-70; e-mail: mvk_home@mail.ru
Мальцев-Горский Дмитрий Александрович, студент Московского авиационного института (государственного технического университета). Авиационная ул., д. 70 к.3, кв. 67, Москва 123182; тел.: 8-926-330-6393; e-mail: netdmg@msn.com
