Теплофизика высоких температур, 2022, T. 60, № 1, стр. 3-10
Электронный ток насыщения на цилиндрический зонд в потоке разреженной замагниченной плазмы
В. А. Шувалов 1, *, Ю. П. Кучугурный 1, Г. С. Кочубей 1, С. В. Носиков 1
Институт технической механики Национальной академии наук Украины (ИТМ)
г. Днепр, Украина
* E-mail: vashuvalov@ukr.net
Поступила в редакцию 22.04.2021
После доработки 04.08.2021
Принята к публикации 28.09.2021
- EDN: USJYNV
- DOI: 10.31857/S0040364422010033
Аннотация
Разработана процедура диагностики замагниченной разреженной плазмы с использованием электронного тока насыщения на цилиндрический электрический зонд. Получены приближенные формулы и зависимости электронного тока насыщения от угла между осью зонда и вектором индукции внешнего магнитного поля, а также от масштабных параметров, характеризующих собирание зондового тока: радиусов зонда и приэлектродного слоя, ларморовского радиуса электронов плазмы. Показано, что для углов более 65° электронный ток на цилиндрический зонд в замагниченной плазме равен электронному току насыщения в отсутствие внешнего магнитного поля. Представленные формулы и зависимости позволяют определить параметры электронов в потоке разреженной замагниченной плазмы с использованием двух выходных сигналов – зондового тока и потенциала зонда.
ВВЕДЕНИЕ
Электрические цилиндрические зонды широко используются для диагностики потоков разреженной замагниченной плазмы: на летательных аппаратах в верхней атмосфере и ионосфере Земли; в плазменных аэродинамических трубах; в струях, генерируемых электрореактивными двигателями; при исследовании магнитогидродинамических процессов и течений [1–4]. Интерпретацию зондовых измерений затрудняют достаточно сложная теория, описывающая взаимодействие зондов с плазмой, и процессы, сопровождающие собирание тока в системе “зонд‒плазма”. На вольт-амперной характеристике (ВАХ) зонда (зависимости собираемого тока от потенциала зонда) условно выделяются три участка: ветвь ионного тока насыщения, переходный участок и область насыщения электронного тока на зонд.
Собирание тока цилиндрическими электрическими зондами и электродами в потоках разреженной замагниченной плазмы широко обсуждалось в литературе, например [5–8] и др. Установлено, что при отрицательных потенциалах зонда ${{\varphi }_{\rho }}$ относительно потенциала плазмы ${{\varphi }_{\infty }}$ (${{\varphi }_{W}} = {{\varphi }_{p}} - {{\varphi }_{\infty }} < 0$) влиянием внешнего магнитного поля с индукцией ${{B}_{\infty }} \leqslant {{10}^{{ - 2}}}$ Тл на ионную ветвь ВАХ можно пренебречь [9].
Математически ионный ток на зонд в потоке разреженной замагниченной плазмы представляет многопараметрическую функцию. Интерпретация ионного тока на цилиндрический зонд, определение концентрации заряженных частиц (ионов ${{N}_{i}}$) по ионной ветви ВАХ сопряжены с необходимостью учета ряда параметров: зависимости ионного тока от ориентации оси симметрии зонда (${{{\mathbf{l}}}_{p}}$) относительно векторов скорости потока плазмы ${{U}_{\infty }}$ и индукции внешнего магнитного поля ${{{\mathbf{B}}}_{\infty }}$, числа Маха (температуры, химического состава и молекулярной массы частиц), степени неизотермичности плазмы и ряда характеристических длин, таких как радиус ${{r}_{p}}$ и длина ${{l}_{p}}$ зонда, ларморовский радиус ионов ${{r}_{i}}$, дебаевский радиус ${{\lambda }_{d}}$ плазмы, радиус приэлектродного слоя ${{r}_{s}}$.
Собирание ионного тока сопровождается процессами фотоэмиссии (в ионосфере на освещенном участке орбиты летательного аппарата и вторичной ионно-электронной эмиссии. При плотности ионного тока насыщения на цилиндрический зонд ${{j}_{{i{\text{sat}}}}}\sim 10 \times {{10}^{{ - 9}}}$ А/см2 в ионосфере на высотах 500−800 км плотность тока фотоэмиссии для материалов зонда (W, Mо, Au, Pt) лежит в пределах (2–8) × 10–9 А/см2, что составляет 50–70% собираемого ионного тока насыщения [2, 10]. Значения коэффициентов вторичной ионно-электронной эмиссии ионов ${\text{H}}_{{\text{2}}}^{{\text{ + }}}{\text{,}}\,\,{\text{H}}{{{\text{e}}}^{{\text{ + }}}}{\text{,}}$ ${\text{N}}{{{\text{e}}}^{{\text{ + }}}}{\text{,}}\,\,{\text{A}}{{{\text{r}}}^{{\text{ + }}}}{\text{,}}$ ${\text{N}}_{{\text{2}}}^{{\text{ + }}}{\text{,}}\,\,{\text{O}}_{{\text{2}}}^{{\text{ + }}}{\text{,}}\,\,{\text{K}}{{{\text{r}}}^{{\text{ + }}}}$ и ${\text{X}}{{{\text{e}}}^{{\text{ + }}}}$ на поверхностях материалов электрических зондов близки к аппроксимации ${{\gamma }_{i}} \approx 1.6 \times {{10}^{{ - 2}}}({{h}_{i}} - 2{{\chi }_{W}})$, где ${{h}_{i}}$ – потенциал ионизации иона, ${{\chi }_{W}}$ – работа выхода материала зонда [11]. Токи вторичной ионно-электронной эмиссии для большинства ионов с энергией ${{E}_{i}} \leqslant 100$ эВ составляют 5–25% от собираемого зондом ионного тока. Учет перечисленных факторов существенно затрудняет интерпретацию ионной ветви ВАХ, снижает точность определения концентрации ${{N}_{i}}$ заряженных частиц в потоке разреженной замагниченной плазмы.
Собирание электронного тока насыщения ${{I}_{{e{\text{sat}}}}}$ при положительных потенциалах цилиндрического зонда (${{\varphi }_{W}} > 0$) в потоке разреженной замагниченной плазмы зависит от ориентации оси симметрии зонда ${{{\mathbf{l}}}_{p}}$ относительно вектора индукции внешнего магнитного поля ${{{\mathbf{B}}}_{\infty }}$ и от характеристических длин ${{r}_{p}},\,{{l}_{p}},\,{{r}_{s}}$, ${{r}_{e}}$ (${{r}_{e}}$ – ларморовский радиус электронов). Эмиссионные процессы практически не влияют на собирание электронного тока насыщения: при ${{\varphi }_{W}} > 0$ фото-, вторичные и отраженные электроны движутся в тормозящем поле зонда и большая их часть возвращается на поверхность зонда [12, 13]. Плотность электронного тока насыщения ${{j}_{{e{\text{sat}}}}}$ на цилиндрический зонд в потоке разреженной замагниченной плазмы не зависит от угла между осью симметрии ${{{\mathbf{l}}}_{p}}$ зонда и вектором скорости потока ${{{\mathbf{U}}}_{\infty }}$ и практически на 1.5–2 порядка превышает плотность ${{j}_{{i{\text{sat}}}}}$ ионного тока насыщения. При этом количество параметров, характеризующих собирание электронного тока насыщения ${{j}_{{e{\text{sat}}}}}$, значительно меньше и, соответственно, процедура интерпретации электронной ветви ВАХ существенно проще, чем для ионной составляющей ${{j}_{{i{\text{sat}}}}}$ зондового тока.
Цель работы:
– получить приближенные зависимости электронного тока насыщения на цилиндрический зонд в потоке разреженной замагниченной плазмы с использованием элементов теории стока электронов на зонд [14], данных спутниковых и ракетных измерений ВАХ в ионосферной разреженной плазме, а также результатов физического (стендового) эксперимента;
– разработать процедуру определения параметров электронов в потоке разреженной замагниченной плазмы с использованием электронной ветви ВАХ, выходных сигналов ${{j}_{{e{\text{sat}}}}}$ и ${{\varphi }_{p}}$ цилиндрического зонда.
ПАРАМЕТРЫ ЭЛЕКТРОННОГО ТОКА НАСЫЩЕНИЯ
Несмотря на многочисленные публикации, приближенные и численные решения задачи о собирании ионного и электронного тока цилиндрическим зондом в замагниченной бесстолкновительной плазме, данные о непосредственном сравнении расчетных и измеренных значений ВАХ скудны. Расчеты и интерпретация ВАХ зонда в замагниченной разреженной плазме в конкретных условиях измерений затруднительны.
В [15] для электронного тока ${{I}_{{eB}}}$ на зонд произвольной формы при положительных потенциалах, близких к потенциалу плазмы ${{\Phi }_{W}} = {{e{{\varphi }_{W}}} \mathord{\left/ {\vphantom {{e{{\varphi }_{W}}} {k{{T}_{e}} \geqslant 0}}} \right. \kern-0em} {k{{T}_{e}} \geqslant 0}}$ и ${{{{T}_{i}}} \mathord{\left/ {\vphantom {{{{T}_{i}}} {{{T}_{e}}}}} \right. \kern-0em} {{{T}_{e}}}} \ll 1$ получена зависимость
В [5, 6, 14] ток электронов на слабо заряженный положительный (${{\Phi }_{W}} \geqslant 0$) цилиндрический зонд представлен в виде
При высоких положительных потенциалах по результатам измерений электронного тока насыщения IeB(θB, ΦW) на цилиндрический зонд (радиус ${{r}_{p}} = 3.8 \times {{10}^{{ - 2}}}$ см, длина ${{l}_{p}} = 20.3$ см) научного модуля NASA 18.70 [8] в ионосфере на высотах h = 250−340 км для rp/λd = 0.11, 0.021, 0.012; rp/re = = 0.014; re/Δrs = 0.3, 0.45, 2.5 установлено, что при θB ≥ 65°
Электронный ток насыщения на цилиндрический зонд для ${{\theta }_{B}} \geqslant 65^\circ $ и ${{\Phi }_{W}} = 5.0$ практически равен электронному току насыщения на зонд в отсутствие внешнего магнитного поля (${{B}_{\infty }} = 0$). Этот результат согласуется с выводом и расчетами работы [16]: для больших $\left| {{{\varphi }_{p}}} \right|$ ток, ограниченный орбитальным движением электронов, является верхним пределом, не зависит от ${{B}_{\infty }}$ и изменяется как ${{\left| {{{\varphi }_{p}}} \right|}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}$.
В общем виде электронный ток насыщения на цилиндрический зонд, произвольно ориентированный относительно вектора индукции внешнего магнитного поля ${{{\mathbf{B}}}_{\infty }}$, можно представить в виде зависимости
В [8] для оценки отношения ${{\Delta {{r}_{s}}} \mathord{\left/ {\vphantom {{\Delta {{r}_{s}}} {{{\lambda }_{d}}}}} \right. \kern-0em} {{{\lambda }_{d}}}}$ используется соотношение
Для ${{{{r}_{p}}} \mathord{\left/ {\vphantom {{{{r}_{p}}} {{{\lambda }_{d}}}}} \right. \kern-0em} {{{\lambda }_{d}}}} \leqslant 0.17$ с погрешностью менее 4% множитель F ≈ 1.0 и
(1)
${{\Delta {{r}_{s}}} \mathord{\left/ {\vphantom {{\Delta {{r}_{s}}} {{{\lambda }_{d}}}}} \right. \kern-0em} {{{\lambda }_{d}}}} = \Phi _{W}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}.$Из (1) следует
(2)
${{\left( {{{\Delta {{r}_{s}}} \mathord{\left/ {\vphantom {{\Delta {{r}_{s}}} {{{r}_{p}}}}} \right. \kern-0em} {{{r}_{p}}}}} \right)}^{2}} = 1.05 \times {{10}^{{ - 5}}}\frac{{{{l}_{p}}}}{{{{r}_{p}}}}\frac{{\varphi _{W}^{{{3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2}}}}}{{{{I}_{{e{\text{sat}}}}}}}.$При ${{{{r}_{s}}} \mathord{\left/ {\vphantom {{{{r}_{s}}} {{{r}_{p}}}}} \right. \kern-0em} {{{r}_{p}}}} \gg 1$ с погрешностью $ \leqslant {\kern 1pt} 7\% $ ${{\left( {{{\Delta {{r}_{s}}} \mathord{\left/ {\vphantom {{\Delta {{r}_{s}}} {{{r}_{p}}}}} \right. \kern-0em} {{{r}_{p}}}}} \right)}^{2}} \to $ $ \to {{\left( {{{{{r}_{s}}} \mathord{\left/ {\vphantom {{{{r}_{s}}} {{{r}_{p}}}}} \right. \kern-0em} {{{r}_{p}}}}} \right)}^{2}}$. Соотношение (2) подобно закону “3/2” для цилиндрического электрода [17]
При этом в [18–20] для параметров ${{\Delta {{r}_{s}}} \mathord{\left/ {\vphantom {{\Delta {{r}_{s}}} {{{\lambda }_{d}}}}} \right. \kern-0em} {{{\lambda }_{d}}}}$ и ${{\Phi }_{W}}$ получены соотношения, которые могут быть представлены в виде
(3)
${{\Delta {{r}_{s}}} \mathord{\left/ {\vphantom {{\Delta {{r}_{s}}} {{{\lambda }_{d}}}}} \right. \kern-0em} {{{\lambda }_{d}}}} = \gamma \Phi _{W}^{{{3 \mathord{\left/ {\vphantom {3 4}} \right. \kern-0em} 4}}}.$Из (3) следует
(4)
${{\left( {{{\Delta {{r}_{s}}} \mathord{\left/ {\vphantom {{\Delta {{r}_{s}}} {{{r}_{p}}}}} \right. \kern-0em} {{{r}_{p}}}}} \right)}^{2}} = 6.22{{\gamma }^{2}}\frac{{{{l}_{p}}}}{{{{r}_{p}}}}\frac{{\varphi _{W}^{2}}}{{{{V}_{e}}{{I}_{{e{\text{sat}}}}}}}.$При ${{\varphi }_{W}} \gg 1$ в (2), (4) φw ≃ φp. В [18, 21] для цилиндрических зондов приведена зависимость γ = = ${{{\text{1}}{\text{.278}}} \mathord{\left/ {\vphantom {{{\text{1}}{\text{.278}}} {{{{\left[ {\ln \left( {{{\Delta {{r}_{s}}} \mathord{\left/ {\vphantom {{\Delta {{r}_{s}}} {{{r}_{p}}}}} \right. \kern-0em} {{{r}_{p}}}}} \right)} \right]}}^{{1{\text{/}}2}}}}}} \right. \kern-0em} {{{{\left[ {\ln \left( {{{\Delta {{r}_{s}}} \mathord{\left/ {\vphantom {{\Delta {{r}_{s}}} {{{r}_{p}}}}} \right. \kern-0em} {{{r}_{p}}}}} \right)} \right]}}^{{1{\text{/}}2}}}}}$, из которой при 10 < Δrs/rp ≤ ≤ 300 следует ${{\gamma }_{{{\text{mid}}}}} \approx 0.60$. По результатам измерений ВАХ цилиндрических зондов в ионосфере и в лабораторной плазме при $0.17 \leqslant {{T}_{e}} \leqslant 6.15$ эВ и $5 \times {{10}^{4}} \leqslant {{N}_{e}} \leqslant 4 \times {{10}^{9}}$ см–3 [4, 21–23] имеем $0.49 \leqslant $ ≤ $\gamma \leqslant 0.58$, среднее значение ${{\gamma }_{{{\text{mid}}}}} = 0.53$. В то же время из равенства соотношений (1) и (3) $\gamma \approx \Phi _{W}^{{ - {1 \mathord{\left/ {\vphantom {1 4}} \right. \kern-0em} 4}}}$ и для $5 < {{\Phi }_{W}} \leqslant 25$ получим ${{\gamma }_{{{\text{mid}}}}} = 0.51$.
ТЕХНИКА ЭКСПЕРИМЕНТА
Экспериментальные исследования проводились на плазмодинамическом стенде Института технической механики НАН Украины (ИТМ). Стенд относится к классу плазменных аэродинамических труб. Безмасляная откачивающая система производительностью ~50 м3/с, наличие криопанелей, охлаждаемых жидким азотом, обеспечивают в вакуумной камере стенда (цилиндр диаметром 1.2 м и длиной 3.5 м) остаточное давление ~1 × × 10–5 H/м2, а при натекании газа – рабочее давление 10–4−10–3 H/м2. Источником потоков разреженной плазмы служит газоразрядный ускоритель с ионизацией рабочего газа электронным ударом, осцилляцией электронов во внешнем магнитном поле и “саморазгоном” плазмы [3].
Для диагностики потока разреженной плазмы на стенде использовались: микроволновой интерферометр, работающий на частоте 5.45 ГГц, система электрических зондов (цилиндрические, плоский и многоэлектродный зонд-анализатор) и цилиндр Фарадея. Зонды установлены на подвижных платформах (верхней и нижней) с четырьмя степенями свободы каждая. Платформы обеспечивают угловые и поперечные перемещения зондов в горизонтальной и вертикальных плоскостях и вращение вокруг вертикальной оси. Точность отсчета для линейных перемещений ~0.5×10–3 м, угловой ~0.5°. Определение параметров плазмы осуществлялось с использованием вольт-амперных характеристик зондов и сигналов микроволнового интерферометра [24].
В экспериментах использовались три цилиндрических зонда, изготовленных из молибдена, со следующими геометрическими характеристиками: ${{r}_{{p1}}} = 5 \times {{10}^{{ - 2}}}$ см, ${{l}_{{p1}}} = 8 \times {{10}^{{ - 1}}}$ см; ${{r}_{{p2}}} = 4.5 \times {{10}^{{ - 3}}}$ см, ${{l}_{{p2}}} = 9 \times {{10}^{{ - 1}}}$ см; ${{r}_{{p3}}} = 4 \times {{10}^{{ - 3}}}$ см, ${{l}_{{p3}}} = 4 \times {{10}^{{ - 1}}}$ см. При измерениях ВАХ для всех зондов реализован режим бесстолкновительного обтекания потоком плазмы.
Сигнал микроволнового интерферометра не зависит от наличия (отсутствия) внешнего магнитного поля. Концентрация электронов ${{N}_{e}}$ пропорциональна частоте зондирующей волны и фазовому сдвигу, обусловленному присутствием плазмы между антеннами [24].
В [8] показано, что концентрация электронов ${{N}_{e}}$ в потоке разреженной замагниченной плазмы может быть определена по электронному ${{I}_{{e{\text{sat}}}}}$ току насыщения цилиндрического зонда, ось которого ортогональна векторам скорости потока плазмы ${{U}_{\infty }}$ и индукции внешнего магнитного поля ${{{\mathbf{B}}}_{\infty }}$. При ${{\varphi }_{p}} \gg {{\varphi }_{\infty }}$ концентрация может быть определена из соотношения
где $\chi = \frac{{{{A}_{p}}}}{\pi }e{{\left( {{{2e} \mathord{\left/ {\vphantom {{2e} {{{m}_{e}}}}} \right. \kern-0em} {{{m}_{e}}}}} \right)}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}} \approx {\text{const}}{\text{.}}$Участок электронного тока насыщения ${{I}_{{e{\text{sat}}}}}$, соответствующий условию ${{\varphi }_{p}} \gg {{\varphi }_{\infty }}$, не сложно определить на электронной ветви ВАХ, построенной в линейном или полулогарифмическом масштабе.
В [25, 26] экспериментально установлено, что корректное определение температуры электронов ${{T}_{e}}$ в потоке разреженной замагниченной плазмы может быть осуществлено с помощью тонких цилиндрических зондов, собирающая поверхность которых перпендикулярна векторам ${{U}_{\infty }}$ и ${{{\mathbf{B}}}_{\infty }}$.
Для оценки потенциала плазмы ${{\varphi }_{\infty }}$ и температуры электронов ${{T}_{e}}$ могут быть использованы линейный участок электронной ветви ВАХ (φ∞ − φp) и соотношения [22]
Измерения зондового тока проводились в потоке разреженной плазмы азота для двух режимов работы плазменного ускорителя:
I) при температуре ${{T}_{{e1}}} = 1.2 \times {{10}^{4}}$ К, концентрации электронов ${{N}_{{e1}}} = 4.6 \times {{10}^{6}}$ см–3, скорости ионов ${{U}_{{i1}}} \approx 8.3$ км/с и двух значениях индукции внешнего магнитного поля ${{B}_{{11}}} = 150$ Гс, ${{B}_{{12}}} = 15$ Гс;
II) при ${{T}_{{e2}}} = 3.5 \times {{10}^{4}}$ К, ${{N}_{{e2}}} = 2.5 \times {{10}^{9}}$ см–3, ${{U}_{{i2}}} = 10.6$ км/с и ${{B}_{{21}}} = 150$ Гс, ${{B}_{{22}}} = 15$ Гс.
Масштабные коэффициенты, характеризующие собирание электронного тока зондами, приведены в таблице. Там же представлены масштабные коэффициенты для условий измерения электронного тока цилиндрическим зондом в ионосфере на научном модуле NASA 18.70 [8] и на космическом аппарате (КА) “Explorer-31”. При оценке масштабных коэффициентов для КА “Explorer-31” использовались данные из [27–30].
СОБИРАНИЕ ЭЛЕКТРОННОГО ТОКА НАСЫЩЕНИЯ
На рис. 1 показана нормированная на величину ${{I}_{{0e}}} = {{I}_{{e{\text{sat}}}}}\left( {{{\Phi }_{W}} = 0} \right)$ зависимость электронного тока насыщения на цилиндрический зонд, ось симметрии которого параллельна вектору индукции внешнего магнитного поля ${{{\mathbf{B}}}_{\infty }}$, от параметра ${{{{r}_{p}}} \mathord{\left/ {\vphantom {{{{r}_{p}}} {{{r}_{e}}}}} \right. \kern-0em} {{{r}_{e}}}}$:
При высоких положительных потенциалах зонда собирающей поверхностью служит поверхность приэлектродного слоя радиусом ${{r}_{s}}$. На рис. 2 приведены зависимость от ${{{{r}_{e}}} \mathord{\left/ {\vphantom {{{{r}_{e}}} {{{r}_{s}}}}} \right. \kern-0em} {{{r}_{s}}}}$ электронного тока насыщения $i_{{eB}}^{ - } = {{{{I}_{{eB}}}(0,{{\Phi }_{W}} \geqslant 5)} \mathord{\left/ {\vphantom {{{{I}_{{eB}}}(0,{{\Phi }_{W}} \geqslant 5)} {{{I}_{{eB}}}({\pi \mathord{\left/ {\vphantom {\pi 2}} \right. \kern-0em} 2},{{\Phi }_{W}} \geqslant 5)}}} \right. \kern-0em} {{{I}_{{eB}}}({\pi \mathord{\left/ {\vphantom {\pi 2}} \right. \kern-0em} 2},{{\Phi }_{W}} \geqslant 5)}}$, нормированного на его максимальное значение, когда его ось ортогональна вектору индукции внешнего магнитного поля (${{\theta }_{B}} = {\pi \mathord{\left/ {\vphantom {\pi 2}} \right. \kern-0em} 2}$), и аппроксимация авторов
(5)
$i_{{eB}}^{ - }\left( {{{\theta }_{B}} = 0,\,\,{{\Phi }_{W}} \geqslant 5} \right) = {{\left( {1 + \frac{{3\pi }}{{16}}\frac{{{{r}_{s}}}}{{{{r}_{e}}}}} \right)}^{{ - 1}}}.$В слабых магнитных полях, когда ${{{{r}_{e}}} \mathord{\left/ {\vphantom {{{{r}_{e}}} {{{r}_{s}}}}} \right. \kern-0em} {{{r}_{s}}}} \gg 1$ $\left( {{{{{r}_{e}}} \mathord{\left/ {\vphantom {{{{r}_{e}}} {{{r}_{s}}}}} \right. \kern-0em} {{{r}_{s}}}} \geqslant 50} \right)$, нормированный ток на зонд $i_{{eB}}^{ - }(0,{{\Phi }_{W}} \geqslant $ ≥ 5) ≃ 1, т.е. ${{I}_{{eB}}}({{\theta }_{B}} = 0,{{\Phi }_{W}} \geqslant 5) = $ ${{I}_{{eB}}}({{\theta }_{{B{\text{ }}}}} = 0,$ ${{\Phi }_{W}} \geqslant 5)$ $ = {{I}_{{e{\text{sat}}}}}\left( {{{\Phi }_{W}} \geqslant 5} \right)$, где ${{I}_{{e{\text{sat}}}}} = \frac{{{{A}_{p}}}}{4}e{{N}_{e}}{{\bar {V}}_{e}}\frac{2}{{\sqrt \pi }} \times $ $ \times {{\left( {1 + {{\Phi }_{W}}} \right)}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}$. Для ${{\Phi }_{W}} \geqslant 5$ и ${{\varphi }_{p}} \gg {{\varphi }_{\infty }}$ электронный ток насыщения на цилиндрический зонд составляет ${{I}_{{e{\text{sat}}}}} = \frac{{{{A}_{p}}}}{\pi }e{{N}_{e}}{{\left( {2e{{{{\varphi }_{p}}} \mathord{\left/ {\vphantom {{{{\varphi }_{p}}} {{{m}_{e}}}}} \right. \kern-0em} {{{m}_{e}}}}} \right)}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}$ и концентрация электронов ${{N}_{e}}$ определяется по измеренным значениям ${{I}_{{e{\text{sat}}}}}$ и ${{\varphi }_{p}}$.
На рис. 3 представлены зависимости нормированного $i_{{eB}}^{ - }$ электронного тока насыщения на цилиндрический зонд от угла ${{\theta }_{B}}$ между осью ${{{\mathbf{l}}}_{p}}$ и вектором индукции внешнего магнитного поля ${{{\mathbf{B}}}_{\infty }}$ для фиксированных значений безразмерного потенциала ${{\Phi }_{W}} = 5$ (работа [8] и ${{\Phi }_{W}} \approx 12.6$ при измерениях тока $i_{{eB}}^{ - }\left( {{{\theta }_{B}}} \right)$ на стенде). Электронный ток насыщения ${{I}_{{eB}}}\left( {{{\theta }_{B}},\,\,{{\Phi }_{W}}} \right)$ нормирован на его максимальное значение ${{I}_{{eB}}}({\pi \mathord{\left/ {\vphantom {\pi 2}} \right. \kern-0em} 2},{{\Phi }_{W}})$, собираемое зондом, ось которого ортогональна вектору ${{{\mathbf{B}}}_{\infty }}$.
Значения 5 на рис. 3 для ${{\Phi }_{W}} \geqslant 5$ соответствуют аппроксимации авторов
(6)
$\begin{gathered} i_{{eB}}^{ - }\left( {{{\theta }_{B}},{{\Phi }_{W}}} \right) = \frac{{{{I}_{{eB}}}\left( {{{\theta }_{B}},{{\Phi }_{W}}} \right)}}{{{{I}_{{eB}}}\left( {{\pi \mathord{\left/ {\vphantom {\pi 2}} \right. \kern-0em} 2},{{\Phi }_{W}}} \right)}} = \frac{{{{I}_{{eB}}}\left( {{{\theta }_{B}},{{\Phi }_{W}}} \right)}}{{{{I}_{{e{\text{sat}}}}}\left( {{{\Phi }_{W}}} \right)}} = \\ = \sin {\kern 1pt} {{\theta }_{B}} + i_{{eB}}^{ - }\left( {0,{{\Phi }_{W}}} \right){{\left( {1 - \sin {\kern 1pt} {{\theta }_{B}}} \right)}^{{0.7}}}, \\ \end{gathered} $При ${{\varphi }_{p}} \gg {{\varphi }_{\infty }}$, ${{\theta }_{B}} = {\pi \mathord{\left/ {\vphantom {\pi 2}} \right. \kern-0em} 2}$ имеем ${{I}_{{e{\text{sat}}}}} = \frac{{{{A}_{p}}}}{\pi }e{{N}_{e}} \times $ $ \times \,\,{{\left( {2e{{{{\varphi }_{p}}} \mathord{\left/ {\vphantom {{{{\varphi }_{p}}} {{{m}_{e}}}}} \right. \kern-0em} {{{m}_{e}}}}} \right)}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}$, откуда
(7)
${{N}_{e}} = \pi {{{{I}_{{e{\text{sat}}}}}} \mathord{\left/ {\vphantom {{{{I}_{{e{\text{sat}}}}}} {{{A}_{p}}e{{{\left( {2e{{{{\varphi }_{p}}} \mathord{\left/ {\vphantom {{{{\varphi }_{p}}} {{{m}_{e}}}}} \right. \kern-0em} {{{m}_{e}}}}} \right)}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}} = {{\chi }^{{ - 1}}}}}} \right. \kern-0em} {{{A}_{p}}e{{{\left( {2e{{{{\varphi }_{p}}} \mathord{\left/ {\vphantom {{{{\varphi }_{p}}} {{{m}_{e}}}}} \right. \kern-0em} {{{m}_{e}}}}} \right)}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}} = {{\chi }^{{ - 1}}}}}{{I}_{{e{\text{sat}}}}}\varphi _{p}^{{{{ - 1} \mathord{\left/ {\vphantom {{ - 1} 2}} \right. \kern-0em} 2}}}.$Из соотношений (6) и (7) при ${{\Phi }_{W}} \geqslant 5$ или ${{\varphi }_{p}} \gg {{\varphi }_{\infty }}$ следует
Для ${{r}_{e}} = 3.37{{\sqrt {{{T}_{e}}} } \mathord{\left/ {\vphantom {{\sqrt {{{T}_{e}}} } {{{B}_{\infty }}}}} \right. \kern-0em} {{{B}_{\infty }}}}$ получим
(8)
$\frac{{{{B}_{\infty }}}}{{\sqrt {{{T}_{e}}} }} = \frac{{5.7}}{{{\kern 1pt} {{r}_{s}}}}\left( {i_{{eB}}^{{ - 1}}\left( {0,{{\Phi }_{W}} \geqslant 5} \right) - 1} \right).$Здесь ${{T}_{e}}$ в эВ, ${{B}_{\infty }}$ в Гс, re и rs в см.
На орбите космического аппарата “Explorer-31” при витке 546, высоте h = 618 км, наклонении орбиты ${{\beta }_{H}} \approx $ 80° индукция дипольного магнитного поля Земли ${{B}_{\infty }}$ ≈ 0.47 Гс.
Для значений масштабных коэффициентов, характеризующих собирание электронного тока насыщения (см. таблицу) при $i_{{eB}}^{ - }\left( {{{\theta }_{B}} = 0,\,\,{{\Phi }_{W}} \gg 1} \right) \approx 0.49$ (рис. 8 б и 9б в [27]), из (8) получим ${{T}_{e}}$ ≈ 0.23 эВ, что в пределах погрешности ≤5% согласуется со значением ${{T}_{e}}$ ≈ 0.24 эВ (день, средняя солнечная активность [28, 31]).
В то же время для режимов А и В собирания электронного тока насыщения цилиндрическим зондом (эксперимент NASA 18.70, таблица) на высотах h = 250−340 км с наклонением орбиты ${{\beta }_{H}} \approx $ 32° при ${{T}_{e}}$ = (8.6−8.9) × 10–2 эВ и ${{i}_{{eB}}}$ ≈ 0.81 (режим А), ${{i}_{{eB}}}$ ≈ 0.46 (режим В) (рис. 4) значение индукции внешнего магнитного поля Земли ${{B}_{\infty }}$ ≈ 0.36 Гс по (8) с погрешностью ≈3% согласуется с оценкой ${{B}_{\infty }}$ ≈ 0.37 Гс.
С использованием электронных ветвей ВАХ двух взаимно ортогональных цилиндрических зондов (или двух положений цилиндрического зонда относительно вектора индукции внешнего магнитного поля ${{{\mathbf{B}}}_{\infty }}$ при ${{\theta }_{B}} = 0$ и ${{\theta }_{B}} = {\pi \mathord{\left/ {\vphantom {\pi 2}} \right. \kern-0em} 2}$) соотношение (8) устанавливает связь между индукцией магнитного поля ${{B}_{\infty }}$ и температурой электронов плазмы ${{T}_{e}}$. Процедура определения концентрации ${{N}_{e}}$ и температуры ${{T}_{e}}$ электронов в потоке замагниченной плазмы по электронному току насыщения цилиндрического зонда, где ${{\varphi }_{p}} \gg {{\varphi }_{\infty }}$, включает использование соотношений (2), (7), (8) и выходных сигналов ${{I}_{{e{\text{sat}}}}}$ и $\,{\text{ }}{{\varphi }_{p}}$.
В соответствии с данными рис. 3 и выводами работы [8] для углов ${{\theta }_{B}} \geqslant 65^\circ $ при ${{\Phi }_{W}} \geqslant 5$ и ${{r}_{p}} < {{r}_{e}} \leqslant {{r}_{s}}$ электронный ток насыщения на цилиндрический зонд в потоке разреженной замагниченной плазмы не зависит от угла ${{\theta }_{B}}$: $i_{{eB}}^{ - }$ ≃ 1.0 и IeB ≥ Iesat.
Эти выводы подтверждают и результаты измерения электронного тока насыщения на цилиндрический зонд КА “Explorer-31” [27]. Измеренные в [27] (рис. 8 б, 9б) на высоте $h = 618$ км значения $i_{{eB}}^{ - }({{\theta }_{B}},{{\Phi }_{W}})$ при ΦW = const показаны на рис. 4. Значения масштабных коэффициентов для режимов измерения приведены в таблице. Расчетные значения $i_{{eB}}^{ - }({{\theta }_{B}},{{\Phi }_{W}})$ при ΦW = const по аппроксимации (8) коррелируют с результатами измерений электронного тока насыщения на цилиндрический зонд КА "Explorer-31” на высотах h = 2200−2420 км в [27]. Зависимости $i_{{eB}}^{ - }({{\theta }_{B}},{{\Phi }_{W}})$ показаны на рис. 5.
Условия измерений | Режим | Масштабные коэффициенты | ||||
---|---|---|---|---|---|---|
${{{{l}_{p}}} \mathord{\left/ {\vphantom {{{{l}_{p}}} {{{r}_{p}}}}} \right. \kern-0em} {{{r}_{p}}}}$ | ${{{{r}_{p}}} \mathord{\left/ {\vphantom {{{{r}_{p}}} {{{\lambda }_{d}}}}} \right. \kern-0em} {{{\lambda }_{d}}}}$ | ${{{{l}_{p}}} \mathord{\left/ {\vphantom {{{{l}_{p}}} {{{\lambda }_{d}}}}} \right. \kern-0em} {{{\lambda }_{d}}}}$ | ${{{{r}_{p}}} \mathord{\left/ {\vphantom {{{{r}_{p}}} {{{r}_{e}}}}} \right. \kern-0em} {{{r}_{e}}}}$ | ${{{{r}_{e}}} \mathord{\left/ {\vphantom {{{{r}_{e}}} {{{r}_{s}}}}} \right. \kern-0em} {{{r}_{s}}}}$ | ||
Стенд ИТМ | I | 16 | 0.14 | 2.3 | 0.018−2.2 | 0.15−8 |
200 | 0.013 | 2.6 | ||||
II | 100 | 0.15 | 15.4 | 0.01−1.3 | ||
Ионосфера, модуль NASA 18.70, h = 250−340 км |
A | 534 | 0.1 | 58.8 | 0.014 | 2.5 |
B | 0.021 | 11.2 | 0.47 | |||
C | 0.012 | 6.4 | 0.3 | |||
Ионосфера, КА “Explorer-31” |
h = 618 км | 0.02 | 32.1 | 0.0088 | 0.57 | |
h = 2200−2400 км | 1533 | 0.0061 | 9.8 | 0.0043 | 0.72 |
ЗАКЛЮЧЕНИЕ
По результатам стендовых (лабораторных), спутниковых и ракетных экспериментов выявлено влияние взаимной ориентации оси цилиндрического зонда и вектора индукции внешнего магнитного поля на собирание электронного тока насыщения в потоке разреженной замагниченной плазмы. Установлено, что для углов ${{\theta }_{B}} \geqslant 65^\circ $ между осью цилиндрического зонда ${{{\mathbf{l}}}_{p}}$ и вектором индукции внешнего магнитного поля ${{{\mathbf{B}}}_{\infty }}$ электронный ток на зонд равен току насыщения в незамагниченной плазме: влияние внешнего магнитного поля на собирание зондового тока отсутствует. Получены приближенные формулы, зависимости электронного тока насыщения от отношения масштабных коэффициентов, таких как радиус зонда, ларморовский радиус электронов и радиус приэлектродного слоя.
Показано, что полученные формулы и зависимости позволяют определять кинетические параметры электронов в потоке замагниченной разреженной плазмы с использованием электронного тока насыщения, двух выходных сигналов – собираемого зондового тока ${{I}_{{e{\text{sat}}}}}$ и потенциала ${{\varphi }_{p}}$ поверхности цилиндрического зонда.
Список литературы
Котельников В.А., Котельников М.В., Кассин О.В. Зондовые измерения на борту гиперзвукового летательного аппарата // ТВТ. 2020. Т. 58. № 2. С. 175.
Whipple E.S. Potential of Surfaces in Space // Rep. Prog. Phys. 1981. V. 4. № 11. P. 1197.
Шувалов В.А., Токмак Н.А., Кучугурный Ю.П., Резниченко Н.П. Торможение намагниченного тела при взаимодействии его магнитного поля с потоком разреженной плазмы // ТВТ. 2020. Т. 58. № 2. С. 163.
Губский В.Ф. Влияние магнитного поля на измерения концентрации и температуры электронов цилиндрическими зондами в ионосфере Земли // Солнечно-земная физика. 2008. Т. 2. Вып. 12. С. 261.
Каган Ю.М., Перель В.И. О диагностике плазмы в магнитном поле с помощью тонкого цилиндрического зонда // ЖТФ. 1968. Т. 38. № 10. С. 1663.
Бакшт Ф.Ю., Дюжев Г.А., Циркель Б.И., Школьник С.М., Юрьев В.Г., Антонов С.В., Вайнберг Л.И., Казанец Г.И. Зондовая диагностика низкотемпературной плазмы в магнитном поле. Ч. III // ЖТФ. 1977. Т. 47. № 11. С. 2269.
Laframboise J.G., Sonmor L.J. Current Collection by Probes and Electrodes in Space Magnetoplasmas: A Review // J. Geophys. Res.: Space Phys. 1993. V. 98. № A1. P. 337.
Szuszczewicz E.P., Takas P.Z. Magnetosheath Effects on Cylindrical Langmuir Probes // Phys. Fluids. 1979. V. 22. № 12. P. 2424.
Алексеев Б.В., Котельников В.А. Зондовый метод диагностики плазмы. М.: Электроатомиздат, 1988. 240 с.
Смирнова В.В. К теории горячего и фотозонда // Геомагнетизм и аэрономия. 1966. Т. 6. № 2. С. 275.
Шувалов В.А. Об аккомодации энергии газовых ионов на поверхности поликристаллов // ПМТФ. 1983. № 6. С. 17.
Грановский В.А. Электрический ток в газах. М.‒Л.: Гостехиздат, 1952. Т. 1. 432 с.
Шульман А.Р., Фридрихов Б.С. Вторично-эмиссионные методы исследования твердого тела. М.: Наука, 1977. 552 с.
Мальков М.А. Сток электронов на зонд в плазме средних давлений и в плазме в магнитном поле // ТВТ. 1991. Т. 29. № 3. С. 429.
Bohm D., Burhop E.H.S., Massey H.S.W. The Use of Probes for Plasma Exploration in Strong Magnetic Field. In: The Characteristics of Electrical Discharge in Magnetic Field. Ch. 2 / Ed. by Cuthric A., Macerling R.K. N.Y.: McGrow Hill, 1949. P. 13.
Rubinstein J., Laframboise J.G. Upper-bound Current to a Cylindrical Probe in a Collisionless Magnetoplasma // Phys. Fluids. 1978. V. 21. № 9. P. 1655.
Langmuir J., Blodgett K. Current Collection by Space Charge Coaxial Cylinders // Phys. Rev. 1923. V. 22. № 4. P. 374.
Bettinger R., Walker E.H. Relationship for Plasma Sheaths about Langmuir Probe // Phys. Fluids. 1965. V. 8. № 4. P. 748.
Bettinger R.T., Chen A.A. An End Effect Associated with Cylindrical Langmuir Probes Moving at Satellite Velocities // J. Geophys. Res. 1968. V. 73. № 7. P. 2513.
Hester S.D., Sonin A.A. Ion Temperature Sensitive End Effect in Cylindrical Langmuir Probe Response at Ionosphere Satellite Conditions // Phys. Fluids. 1970. V. 11. № 5. P. 1265.
Bettinger R.T. An in situ System for Measurement of Ionospheric Parameters. In: Interaction of Space Vehicles with an Ionized Atmosphere / Ed. by Singer S.F. London: Pergamon Press, 1965. P. 163.
Шувалов В.А., Письменный Н.И., Лазученков Д.Н., Кочубей Г.С. Зондовая диагностика потоков лабораторной и ионосферной разреженной плазмы // ПТЭ. 2013. № 4. С. 98.
Хазен А.М., Шувалов В.А. Применение термоанемометра – зонда Ленгмюра для диагностики разреженной плазмы // ТВТ. 1969. Т. 7. № 5. С. 866.
Шувалов В.А., Чурилов А.Е., Турчин В.В. О диагностике струи разреженной плазмы с применением зондового и СВЧ-методов // ТВТ. 1978. Т. 16. № 1. С. 9.
Носачев Л.В., Скворцов В.В. Исследование параметров замагниченного потока синтезированной плазмы // ЖТФ. 1978. Т. 48. № 1. С. 49.
Носачев Л.В., Скворцов В.В. Характеристики зондов замагниченного потока синтезированной плазмы // ЖТФ. 1978. Т. 48. № 11. С. 2319.
Miller N.J. Some Implications of Satellite Spin Effects in Cylindrical Probe Measurements // J. Geophys. Res. 1972. V. 77. № 16. P. 2851.
Wrenn G.L., Smith P.A. Results Derived from Simulations Measurements Using the Langmur Plate and Spherical Ion Probe on Explorer-XXXI and the Ionosonde on Alouette // Proc. IEEE. Spec. Iss. 1969. V. 57. № 6. P. 1085.
Findlay J.S., Brace L.A. The Cylindrical Electrosta-tic Probes Employed on Alouette-II and Explorer-XXXI Satellites // Proc. IEEE. Spec. Iss. 1969. V. 57. № 6. P. 1054.
Donley J.L., Brace L.N., Findlay J.A., Hoffman J.H., Wrenn G.L. Comparison of Results on Explorer-XXXI Direct Measurement Probes // Proc. IEEE. Spec. Iss. 1969. V. 57. № 6. P. 1078.
Гуревич А.В., Шварцбург А.Б. Нелинейная теория распространения радиоволн в ионосфере. М.: Наука, 1973. 272 с.
Дополнительные материалы отсутствуют.
Инструменты
Теплофизика высоких температур