Механика и электродинамика пристеночной плазмы

Плазмой называется ионизованный газ, в котором выполняется условие квазинейтральности. Это условие может соблюдаться лишь в объемах, характерные размеры которых велики в сравнении с некоторым масштабом длины М^ и на временах, больших по сравнению с масштабом времени Как показано в классических монографиях по физике плазмы [^[З], [44], в качестве Мь выбирается радиус дебаевского экранирования гр, а в качестве Гу^ - период плазменных (или ленгмюровских) колебаний. В пристеночной области условие квазинейтральности может не соблюдаться. Причин для этого много: влияние потенциала тела, теплового и направленного движение заряженных частиц, различных видов эмиссии электронов или инжекции других частиц с поверхности, химических реакций и т. д. Между тем роль пристеночных слоев, как правило, оказывается определяющей при создании различных приборов и устройств, использующих плазму в качестве рабочего тела (плазменные движители, технологические плазмотроны, плазмохимические реакторы, МГД-генераторы, плазменная электроника), или работающие в плазме (космические летательные аппараты (ЮТА) и космические станции (КС), гиперзвуковые летательыне аппараты (ГЛА), системы типа «Буран», «Шатл», «Колумбия» и др.). Особое место занимают методы зондовой диагностики, основанные на детальном исследовании процессов переноса заряда в пристеночной области.

Многочисленные работы в области механики и электродинамики пристеночной плазмы, обзор которых будет дан ниже по главам, не полностью охватывают важные для практики режимы течения, часто исследуют асимптотические и предельные случаи, не носят систематического характера. В диссертации исследованы с единых позиций все три режима течения: молекулярный, переходный и континуальный. Обнаружены и физически обоснованы новые нелинейные эффекты, возникающие при взаимодействии заряженных тел с плазмой. Особое внимание уделено приложениям результатов исследований в области авиационно-космической техники и зондовой диагностики. Основной метод исследования — математическое моделирование. Приведены обширные вычислительные эксперименты в различных режимах течения с учетом разнообразных физических процессов. Среди использованных вычислительных методов особое место занимает алгоритм метода крупных частиц Давыдова Ю. М. и его модификации применительно к задачам механики и электродинамики пристеночной плазмы. Результаты работы опубликованы в трех монографиях, десяти статьях в журналах, рекомендованных ВАК, доложены на международных и всесоюзных конференциях (более 25 опубликованных докладов и тезисов докладов). Ряд результатов используется в заинтересованных организациях и в учебном процессе МАИ.

Диссертация состоит из введения, четырех глав, списка использованных источников и приложения.

Основные результаты, полученные в диссертации.

1. С единых позиций в достаточно общей постановке сформулированы физические, математические и вычислительные модели механики и электродинамики пристеночной плазмы в молекулярном, столкновительном и переходном режимах. Показано, что в столкновительной плазме граничное условие на теле для концентрации притягивающихся частиц не может записываться в виде nj e = 0, как это делается практически во всех работах по данной тематике. Если выполняется неравенство rp > X > А (гр — размер телаX — средний пробегА — толщина слоя объемного заряда), то граничное условие имеет вид п, е = пь где П) > 0,2 а".

Разработаны пакеты прикладных программ (ППП) для решения задач механики и электродинамики пристеночной плазмы в молекулярном, континуальном и переходном режимах.

Все ППП оптимизированы и приспособлены для использования на ЭВМ средней мощности, в частности, показано, что условие Куранта-Фридрихса-Леви для шага по времени At в задачах электродинамики пристеночной плазмы оказывается слишком жестким. Шаг по времени в алгоритме метода крупных частиц может быть увеличен в (3-Н-) раза по сравнению с его значением, вытекающим из условия Куранта без потери устойчивости решения. ППП устойчиво работают в широком диапазоне изменения направленной скорости, индукции магнитного поля, характерного размера и потенциала тела и позволяют достичь установления возмущенной зоны за приемлемое машинное время (2-^-8 часов на компьютерах с процессором Pentium 4).

ППП содержит оригинальную графическую оболочку, позволяющую визуально изучать и контролировать эволюцию любых параметров плазмы, которые подлежат определению.

2. В молекулярном режиме получены функции распределения (ФР) заряженных частиц вблизи цилиндрического тела, обтекаемого поперечным потоком плазмы в лобовой, боковой и теневой областях. Если в лобовой и боковой областях с учетом сдвига центров тяжести ФР на V", они похожи на ФР в покоящейся плазме, то в теневой области ФР ионов разделяются на две части, причем каждая из частей соответствует потокам, обтекающим цилиндр с одной стороны и с другой. Этот эффект обнаружен впервые.

Наложение осевого магнитного поля на цилиндр еще более осложняет вид ФР ионов. В некотором интервале изменения индукции магнитного поля двугорбый характер ФР ионов смещается под действием сил Лоренца в боковую область.

Исследовано также влияние на ФР ионов и электронов столкновений типа «ион-нейтрал» и «электрон-нейтрал», которые играют роль при промежуточных числах Кнудсена. Если в покоящейся плазме при отсутствии столкновений ФР имели в проекции на плоскость (Уг и Уе) вид подковы, то столкновения приводят к частичному заполнению внутренней области подковы.

3. В молекулярном и континуальном режимах для случая, обтекания цилиндра поперечным потоком плазмы получены поля скоростей и концентраций заряженных частиц, а также самосогласованного электрического поля в лобовой, боковой и теневой областях. Детально исследована область ближнего следа. Если в молекулярном режиме распределение концентраций притягивающихся частиц имеет монотонный характер с максимумом в лобовой части и минимумом в теневой, то в режиме сплошной среды обнаружен диапазон изменения параметров задачи, в котором концентрации и плотности токов в боковой части цилиндрического тела имеют локальный максимум и минимум, возникновение которых объясняется совместным влиянием конвекции, диффузии и подвижности. Выявлены также условия, при которых концентрация и плотность тока в теневой области имеют локальный максимум, что объясняется вихревым движением нейтральной компоненты, которая увлекает за собой заряженные частицы. Исследованы также предельные случаи, когда превалирующее влияние имеет направленная скорость или электростатическое притяжение заряженных частиц к телу.

4. Ряд новых нелинейных эффектов обнаружен при совместном учете направленной скорости потока и внешнего осевого магнитного поля. Как в молекулярном, так и в континуальном режимах интегральный ионный ток на цилиндр нелинейно изменяется в зависимости от индукции осевого магнитного поля. В определенном интервале изменения параметров задачи он может с ростом поля как увеличиваться, так и уменьшаться. При этом точка, в которой достигается максимум плотности тока, смещается с изменением индукции В по угловой координате 0 (см. рис. 18а). Локальные максимумы и минимумы, которые имели место в столкновительном режиме без магнитного поля, имеют место и в присутствии магнитного поля. При этом с одной стороны цилиндра они смещаются в сторону большей плотности тока (рис. 18а), а с другой — в обратную сторону, образуя довольно сложную картину зависимости плотности тока от угловой координаты 0. Объяснить эти эффекты удается совместным влиянием направленного движения и электромагнитных сил.

5. Получен необходимый для практики зондовых измерений набор вольт-амперных характеристик цилиндрических зондов, расположенных поперек потока и плоских пристеночных зондов, расположенных на большой обтекаемой плазмой пластине.

Характеристики получены как в бесстолкновительном, так и в столкновительном режимах.

Разработано несколько новых методов обработки зондовых характеристик, в частности, для цилиндрического зонда, расположенного в поперечном потоке бесстолкновительной плазмы;

— для цилиндрического зонда в поперечном потоке столкновительной плазмы, если слой объемного заряда не бесконечно тонкий и число Рейнольдса имеет не слишком большое значение;

— даны рекомендации по зондовым измерениям в турбулентном режиме;

— рекомендации по использованию двойных зондов с целью уменьшения их взаимного влияния;

— рекомендации по использованию плоских ориентированных зондов в бесстолкновительном потоке плазмы;

— создан вариант автоматизированной системы обработки зондовых характеристик;

— обоснована возможность корректного проведения зондовых изменений в следе за спутником или телом, находящимся в потоке столкновительной плазмы. 1.

6. Разработана теория и методика применения нестационарного зонда (включая зондовую измерительную аппаратуру) в молекулярном и столкновительном режимах, позволившая расширить возможности зондового метода диагностики плазмы: показана возможность измерения еще одного параметра — температуры ионов, что в классической зондовой теории считалось невозможным.

В вычислительных экспериментах получены необходимые для обработки зондовых характеристик теоретические значения времени релаксации в плазме т в зависимости от радиуса и потенциала тела, отношения температур ионов и электронов, величины направленной скорости потока, индукции внешнего магнитного поля, плотности тока эмиссии электронов и других параметров.

7. Сформулирована физическая, математическая и вычислительная модели пристеночного слоя плазменного якоря как основного элемента электромагнитного рельсового ускорителя. Получено распределение потенциалов и концентраций заряженных частиц в пристеночном слое в зависимости от плотности тока эмиссии с катода.

8. С помощью математических и вычислительных моделей, разработанных в гл. 1 диссертации, исследована структура возмущенной области вблизи спутника, движущегося в ионосферной плазме и параметры собственной атмосферы при наличии фотоэмиссии и вторичной эмиссии с его поверхности. Полученные данные позволяют изучать влияние спутника на другое тело, попавшее в его «след». Эти же данные позволяют построить теорию электрического зонда для проведения зондовых измерений в возмущенной зоне, в том числе и в «следе».

9. На основе математических и вычислительных моделей, изложенных в гл. II диссертации разработано два электромагнитных метода воздействия на параметры ионизованного пограничного слоя с целью его радиопросветления:

— инжекция в пограничный слой с борта гиперзвукового летательного аппарата потока отрицательных ионов;

— возбуждение в пограничном слое импульса внешнего магнитного поля, поперечного относительно скорости потока плазмы.

С помощью вычислительных экспериментов показано, что оба метода технически осуществимы и позволяют понизить концентрацию электронов на порядок и более.

10. Предложена методика оценочного расчета электродинамического способа управления вектором тяги плазменного движителя. Для отклонения потока создается поперечное магнитное поле. Методами математического моделирования рассчитана зависимость угла поворота от величины индукции магнитного поля и других параметров потока плазмы.