автореферат диссертации по авиационной и ракетно-космической технике, 05.07.05, диссертация на тему:Математическое моделирование распространения струи стационарного плазменного двигателя в объеме вакуумной камеры

4851493

ТОРОПОВ Григорий Петрович

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РАСПРОСТРАНЕНИЯ СТРУИ СТАЦИОНАРНОГО ПЛАЗМЕННОГО ДВИГАТЕЛЯ В ОБЪЕМЕ ВАКУУМНОЙ

КАМЕРЫ

Сп ец и ал ьно сть 0 5.07.05 «Тепловые, электроракетные двигатели и энергоустановки летательных аппаратов»

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

1 4 ИЮп 2011

Москва-2011

4851493

Работа выполнена на кафедре «Электроракетные двигатели, энергетические и энергофизические установки» Московского авиационного института (государственного технического университета)

Научный руководитель:

д.т.н., доцент Хартов Сергей Анатольевич

Официальные оппоненты:

д.т.н., профессор Рыжов Юрий Алексеевич, K.T.H., с.н.с. Бишаев Андрей Михайлович

Ведущая организация:

Федеральное государственное унитарное предприятие ОКБ "Факел'

г"

Защита диссертации состоится « 19 » сентября 2011 года в 13 : 00 на засед нии диссертационного совета Д 212.125.08 при Московском авиационном институ (государственном техническом университете) «МАИ» по адресу: 125993, ГСП-3, / 80, Волоколамское ш. д. 4.

С диссертацией можно ознакомиться в библиотеке МАИ.

Отзыв на автореферат в одном экземпляре, заверенный печатью просим в сылать по адресу: 125993, ГСП-3, А-80, Волоколамское ш. д. 4, Ученый совет MAI

Автореферат разослан « сн

2011 года.

Ученый секретарь диссертационного совета Д 212.125.08, д.т.н., профессор

--- Ю.В. Зу

Общая характеристика работы

Актуальность темы:

В космической технике все более широкое применение в качестве тяговых модулей находят электроракетные двигатели (ЭРД) и в частности одна из их разновидностей - стационарный плазменный двигатель (СПД). Наиважнейшим вопросом разработки и совершенствования СПД, как и любого ЭРД, является проблема наземных испытаний. Одним из основных недостатков этого этапа отработки двигателей является невозможность в полном объеме имитировать условия окружающей среды, характерной для работы ЭРД на борту космического аппарата. Так, например, концентрация частиц остаточной атмосферы в вакуумных камерах на несколько порядков превосходит реальную концентрацию частиц в космическом пространстве. Отличны также от наблюдаемых при эксплуатации двигателей: пространственное распределение частиц и состав остаточной атмосферы испытательных стендов. Причиной этого является процесс взаимодействия ускоренных ионов со стенками вакуумных камер, в результате которого образуются обратные потоки атомов рабочего тела, нейтрализовавшегося на стенках, и распыленного с них материала. Эти потоки зависят от геометрии и размеров камеры. Верхний предел распыленных потоков во многом формирует требования для минимальных габаритов вакуумных камер при заданном уровне мощности испытуемого двигателя. Осаждение атомов распыленного материала стенок вакуумных камер может негативно влиять на функционирование изоляторных и электродных узлов ЭРД. В частности, для СПД это оказывает существенное влияние на электропроводящие и эмиссионные свойства поверхности керамической разрядной камеры.

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

Комплексные экспериментальные исследования указанных явлений при испытаниях СПД достаточно трудоемки и не всегда дают однозначный результат. Вместе с тем отдельные процессы относительно хорошо изучены (свободномолекулярное движение атомов; распыление материалов; перенос и осаждение распыленного материала; резонансная перезарядка и т.п.) и поддаются математической формализации. В последнее время были выполнены работы по моделированию поведения струи СПД в вакуумных камерах, которые провели группы исследователей: под руководством И. Бойда (LDJBoyd) в Мичиганском у университете (США), сотрудниками фирмы Альта (Alta, Италия) и ряд др. В данных V /

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

В качестве цели работы была выбрана разработка методики моделирования процесса распространения струи СПД в объеме вакуумной камеры в трехмерном приближении, с учетом взаимодействия ионов с остаточным газом и стенками вакуумной камеры.

Для достижения поставленной цели в работе решались следующие задачи:

1. Разработка трехмерной математической модели распространения струи стационарного плазменного двигателя в вакуумной камере.

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

3. Создание математической модели для расчета процесса осаждения распыленных частиц материала стенок вакуумной камеры на поверхностях двигателя и вакуумной камеры.

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

5. Реализация модели на примере расчета вакуумных камер различной геометрии.

Научная новизна работы состоит в том, что разработана и реализована в виде вычислительного комплекса трёхмерная инженерная модель распространения струи стационарного плазменного двигателя в вакуумной камере. Предложен виртуальный датчик - ПМА, позволяющий проводить расчеты интегральных характеристик СП, (тяга, суммарный ионный ток), а также локальных значений параметров в стру двигателя (давление, плотность ионного тока и распределение частиц по скоростям). Практическая значимость работы состоит в следующем:

- разработана методика математического моделирования процесса распространен« струи СПД, позволяющая изучать влияние на динамику атомов и ионов таки факторов, как геометрия стенок и распределение нейтрального газа в объем вакуумной камеры, а также определять характеристики потоков ионов нейтральных атомов внутри вакуумной камеры в процессе работы двигателя;

- проведено моделирование для конкретных геометрий вакуумных каме динамического распределения нейтральных частиц до и после включения двигателя;

- разработана методика численного определения локальных характеристик

заданных точках струи СПД;

- разработана методика расчета распределения перезарядившихся частиц в объеме вакуумной камеры;

- разработана методика математического моделирования процессов распыления стенок вакуумной камеры и перепыления распыленного материала;

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

Достоверность полученных результатов подтверждается соответствием между модельными и экспериментальными данными, а также совпадением полученных автором результатов с расчетами других исследователей. На защиту выносятся:

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

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

3. Результаты моделирования распределения нейтральных атомов и ионов в вакуумных камерах различной конфигурации.

4. Трехмерная инженерная модель оценки процессов распыления и перепыления материала стенок вакуумной камеры.

Апробация работы и научные публикации. Результаты работы прошли апробацию на 9 и 11 Международных научно-технических конференциях студентов и аспирантов «Радиоэлектроника электротехника и энергетика» в 2003 и 2005гг.; на «Всероссийском конкурсе на лучшие научные работы студентов по естественным, техническим наукам (проекты в области высоких технологий) и инновационным научно-образовательным проектам» в 2004г.; на 3 и 6 Международных конференциях «Авиация и космонавтика» в 2004 и 2007гг.; на Всероссийской конференции молодых ученых и студентов «Информационные технологии в авиационной и космической технике» в 2008г.; на 29, 30 и 31 Международных конференциях по электроракетным двигателям в 2005, 2007 и 2009гг. Основное содержание и результаты диссертационной работы опубликованы в трех статьях (в журналах, рекомендованных ВАК РФ) и 5-ти научно-технических отчетах.

Структура и объем работы. Диссертация изложена на 99 страницах, включая размещенные в тексте 77 рисунка и 2 таблицы, а также списка использованных источников из 55 наименований; и состоит из введения, четырех шав, заключения (выводов) и приложений, объемом 20 стр.

Основное содержание работы

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

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

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

Рассмотрены созданные другими авторами модели струй СПД. Показано, что их основным недостатком является двумерная реализация.

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

Во второй главе представлено описание стационарной трехмерной математической модели распространения струи плазменного двигателя в вакуумной камере. Особенностями реализации ее составных частей являются:

Моделирование давления в вакуумной камере. Распределение нейтральных частиц имеет ключевое значение для анализа основных процессов, происходящих в вакуумной камере: для определения концентрации ионов, образующихся в результате перезарядки, электрического потенциала пространства, плотности тока первичных ионов в струе и др. Динамический вакуум в системе образуют атомы, двигающиеся в объеме камеры, при условии равенства потоков частиц из работающего двигателя и частиц, выходящих в средства откачки. Наиболее эффективным методом для оценки распределения нейтральных частиц, является трассировка их траекторий с использованием метода Монте-Карло. При этом вводятся следующие ограничения:

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

- траектория движения частицы строится до места ее пересечения со стенками вакуумной камеры;

- при попадании на стенку ион рабочего тела рекомбинирует и уходит от нее в виде атома со скоростью, определяемой с учетом максвелловского распределения для температуры стенки;

- т.к. поверхность стенок камеры имеет шероховатость, то при попадании на неё частицы точно неизвестно, в каком направлении она отразится. Поэтому в алгоритме расчета используется случайное отражение по «закону косинуса» относительно нормали к точке касания;

- атомы двигаются строго по прямолинейным траекториям, для учета соударений вводится коэффициент ослабления потока;

- при прохождении атомом расстояния, большего длины свободного пробега, генерируется его соударение с изменением скорости и направления (для данного расчета предварительно производится предрасчет концентрации нейтральных частиц без учета соударений). Сечение упругого взаимодействия атомов ксенона определяется с помощью модели взаимодействия двух шаров разных скоростей методом «твердых сфер», описанной в главе 3 диссертации. Скорость и направление движения частиц задается из базы данных заранее сделанных предрасчетов;

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

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

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

^) = 4гг-(^)3/2-ехр , (1)

где V - хаотическая скорость частицы, Т - средняя температура отраженных частиц. Для корректного использования максвелловского распределения в рассматриваемых моделях подбирается для него степень свободы «Хи-распределения». В данной работе был использован метод наименьших квадратов (см. рис.1).

Рисунок 1. Результат подбора Хи-распределения к Максвелловскому по скоростям для заданной температуры газа 300К

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

Моделирование распределения концентрации ионов перезарядки. Количество перезаряженных атомов в расчетной ячейке определяется при прохождении моделируемого первичного иона через неё. Зная точную траекторию движения моделируемого иона в ячейке и значение сечения перезарядки для заданной энергии иона, можно определить среднюю длину пробега частиц в среде с учетом того, что скорость движения атомов относительно ионов пренебрежимо мала:

где Я - средняя длина пробега иона, п — концентрация нейтральных частиц в ячейке, Опер - сечение перезарядки:

амр = (А-В1о8(Е))-10Г20, (3)

Где £-энергия первичного иона, А = 87,3±0,9, В = 13,6±0,6.

С использованием соотношения (2) можно определить, какое число ионов в

8

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

Моделирование электрического поля. Основной задачей данного расчета является определение распределения потенциала в объеме вакуумной камеры. Граничные условия задачи: начальный потенциал пространства и температура электронов в определенной точке. Для их задания используются данные экспериментальных измерений в фиксированной точке вакуумной камеры. Принимается, что плазма квазинейтральна. Для описания ее параметров используются трл основных уравнения магнитогидродинамики (уравнение неразрывности для электронного газа, уравнение сохранения импульса, уравнение сохранения энергии), которые решаются с помощью метода конечных разностей:

1. Уравнение неразрывности:

|-пе + Ч(пеУе) = пеп„ < <х^е >, (4)

где Уе - направленная скорость движения электрона; - тепловая скорость движения электрона.

2. Уравнение сохранения импульса:

| (тепеК)+тЛ(К*)К =-еп.Е-Уре +ЯС, (5)

с учетом: ре=пекТе, Яс=—, = ' де = < ЪМ > пе+< оеаре > пе

СУ "1е1ае

преобразуется в обобщенный закон Ома:

-Ч<р+—ЩпекТе) еп„

Из соотношения (6) с учетом заданных пе, Ve и Те и уравнения неразрывности тока электронов:

У/ = 0 (7)

определяется распределение потенциала плазмы, применяя конечно-разностную схему для вычисления градиентов и лапласианов, при этом расчет параметров на новой итерации (слое) производится с помощью значений на предыдущей итерации (слое).

3.Уравнение сохранения энергии:

а /3 \ 3

-(-nekTe) + -пе(уеЧ)кТе + peWe =

VAeVTe +;Е -3^ûenek(Te - Т„)-пеп„ < atve > et, (8)

- 2.4 к2П.Те

где ар = —з— --—.

е .. «et méS

Ше

В расчетах использовалось распределение Больцмана для определения потенциала плазмы:

<Р = <Р*+-еТе'\пф) (9)

С учетом этого определяется влияние напряженности электрического поля на траекторию движения частиц. Например, для оси X движение описывается соотношениями:

2

^=<7Ех; Рх=та/, Ях = У0х(+Ц~ (10)

Результаты расчета электрического поля учитываются при перерасчете траектории движения перезарядившихся частиц.

В третьей главе представлена программная реализация математической модели.

В первом разделе дано краткое описание приложений и библиотек, созданных для реализации проекта.

Во втором разделе приводятся математические выкладки, которые использовались при разработке вспомогательной численной сущности -виртуального датчика - поверхностного математического анализатора (ПМА), используемого при расчетах различных интегральных характеристик (тяга, суммарный ионный ток), а также локальных значений параметров (давление, плотность ионного тока, распределение частиц по скоростям):

Определение давления. В процессе моделирования на поверхность ПМА попадает N условных моделируемых частиц. Масса, которая переносится каждой условной частицей за секунду времени - т/И(общ), где Щобщ) - общее число моделируемых частиц. Формула расчета давления принимает следующий вид:

где Vtfzj - проекция вектора скорости иона на нормаль к площадке ПМА, S -площадь поверхности анализатора.

Определение тяги. Для расчета тяги двигателя использовалась следующая методика: назначалось измерительное сечение (ПМА) в компьютерной модели, на расстоянии нескольких сантиметров от двигателя, так, что оно перекрывало все сечение камеры. При попадании на площадку ПМА ускоренного иона, фиксировалось значение его скорости и косинус угла между направлением движения иона к нормали площадки по абсолютной величине. Тяга в таком случае будет равна:

T = ^^Vi'Cos(aùl (12)

Определение плотности тока первичных ионов. Плотность тока ионов на ПМА определяется из формулы:

1 W(o6uÙ-M-S

где F - постоянная Фарадея; ц - молярная масса ксенона; S - площадь ПМА; Z -валентность иона. Поскольку ПМА устанавливаются на определенном расстоянии от двигателя, пучок с моделируемыми ионами, проходя через среду нейтральных атомов в вакуумной камере, перезаряжается, что учитывается с помощью коэффициента ослабления.

Определение плотности тока ионов перезарядки. По результатам моделирования распределения первичных ионов в вакуумной камере известно общее количество ионов, образующихся в результате перезарядки в отдельных ячейках. Их траектории строятся с учетом уравнения движения. В процессе моделирования на поверхность ПМА попадает N условных частиц, Nnep - общее число ионов, образующихся в результате перезарядки в объеме камеры, Л^ -общее число моделируемых траекторий. Плотность тока на анализаторе определяется из формулы:

где е - элементарный заряд, S - площадь ПМА.

В третьем разделе рассмотрен процесс упругого взаимодействия атомов ксенона. Представлены расчеты сечений упругого взаимодействия с использованием различных теорий (модели твердых сфер; модели твердых сфер переменного диаметра; модели с потенциалом взаимодействия по Леннард-Джонсу; модели с потенциалом отталкивания по Борну-Майеру). Проведенный анализ показал, что при температурах менее 1000К целесообразно использовать модель твердых сфер, а при больших температурах, модель с потенциалом отталкивания по Борну-Майеру.

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

В пятом разделе приведен тест модели расчета распыления на полусфере, с использованием датчиков ПМА. Полученные результаты позволили убедиться в правильности построения модели процесса перепыления материала стенок вакуумной камеры в трехмерном пространстве.

В шестом разделе рассмотрены оценки погрешностей модели. Точность расчета давления в вакуумной камере зависит от погрешности датчика определения массового расхода газа, погрешности определения коэффициента выхода для насосов, оценки температуры стенки камеры и косвенно зависит от измеренного экспериментально углового распределения ионов, выбора модели отражения атома от стенки камеры. Оценки показывают, что суммарная погрешность входных данных может быть ~10%, при этом погрешность самого расчета не превосходит 30%. Это сопоставимо с экспериментальными данными - у современных вакуумметров погрешность измерения составляет 10...30%.

Погрешность измерения ионных токов на периферии струи составляет 10...20%. Согласно оценкам суммарная погрешность при расчете концентрации или тока струи в определенном участке камеры, не превышает 10...20%.

Погрешность вычисления скорости распыления атомов и скорости перепыления атомов стенок вакуумной камеры порядка 10...30%, а плотность ионного тока перезарядившихся ионов рассчитывается с точностью 30...50%.

В четвертой главе приведены результаты моделирования распространения струи нескольких моделей СПД в различных вакуумных камерах таблицы 1 и 2, а также тестирования отдельных элементов программного продукта. В приложении А дано подробное описание начальных и граничных условий. Представлено пять различных двигателей (среди них стандартные модели СПД-70, СПД-100 и их аналоги),

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

Таблица 1.

Наименование, изготовитель т., мг/с Up, В Np, Вт Вак. камера,

СПД-70, ОКБ «Факел» 2,35 300 690 1

СПД-70, ОКБ «Факел» 2,35 300 690 6

КМ-5, «Центр Келдыша» 3,5 500 1350 2

Аналог СПД-100, Италия 5,4 300 1500 4

ВРТ-4000, США 14,1 300 3000 5

СПД-100, ОКБ «Факел» 5 300 1350 3

Таблица 2.

Обозн. Название (тип средств откачки) Габариты [диаметр/ длина], м Объем, м3 Площадь насосов, м2

1 Стенд 2ИУВ, кафедра 208 МАИ (масл.) 0,9 х 3,5 3 1,28

2 Стенд КВУ-90, «Центр Келдыша» (крио.) 3,7x10 89 2,51

3 Стенд 71-3-90, ОКБ «Факел» (крио.) 2x7 50 ~1

4 Стенд Centrospazio (Alta, Италия) (крио.) 1,9x4 7,7 1,17

5 Вакуумная камера 12 V, Arnold Engineering Development Center (США) (крио.) 5х 12 164 41

6 Стенд в У2В, кафедра 208 МАИ (масл.) 2x6 20 1,28

7 Стенд в Мичиганском университете (США) (крио.) 6x9 232 7,87

Распределение нейтральных аггомов при работающем двигателе проанализировано с использованием данных экспериментов в четырех установках: 1, 2, 3 и 7, а также сравнено с расчетами других исследовательских групп: фирмы Alta и Мичиганского университета. Достигнута удовлетворительная сходимость с экспериментом (в пределах погрешности измерения датчиков вакуумметра).

На рисунке 2 представлен результат расчета концентрации нейтральных атомов в вакуумной камере 7 (использовалось уравнение Менделеева-Клайперона для пересчета концентрации частиц в давление). Расчет для данной камеры, выполненный в Мичиганском университете, совпадает с экспериментом по форме кривой распределения, однако абсолютные значения в 8...10 раз меньше экспериментальных. Результаты моделирования настоящей работы демонстрируют хорошую сходимость с экспериментом (рис. 2).

1,5Е-03

1.3Е-03

1ДЕ-03

I Э,0Е-04

7,0Е-04

5,0Е-04

3,0Е-О4

|\

11 11 11

11 11 \4 1 и г чч \

\ 1 А \ --и / 1 1

✓ / у' * г ¡¿а ■■ *ч ; V

✓ г *

-3 -2 -1 0 1 2 3 4 5 6 7

удаление от двигателя, м

Рисунок 2. Сравнение результатов расчетов с экспериментом (пунктир -расчет; сплошная - давление, измеренное экспериментально)

Второй раздел посвящен результатам расчета распространения струи плазменного двигателя в объеме вакуумной камеры. Для анализа модели и сравнения с экспериментом был создан модуль, позволяющий имитировать зондовую систему измерений с помощью ПМА. Моделировались две дуговые арки с датчиками: одна с размерами датчиков такими же, как диаметр зонда в эксперименте (7мм), расположенными через 2,5 градуса. Вторая с ПМА диаметром 35мм, закрепленными через 6 градусов. Первая использовалась для расчета ионного тока, вторая для определения энергетического распределения первичных ионов.

На рис.3 и 4 приведены результаты расчета потока ускоренных ионов в вакуумной камере 2. Полученные данные расчета ионного тока имеют хорошую сходимость с экспериментом (рис. 3), а представленные на рис. 4 распределения подтверждают существующие представления о постепенной перестройке кольцевой структуры струи СПД и соответствуют измеренным данным.

угол от оси а, радиан

Рисунок 3. Сравнение плотности моделируемого тока на 91 ПМА диаметром 18 мм (ромбы), находящихся на расстоянии 1м, с экспериментальными измерениями (сплошная линия)

Рисунок 4. Результаты расчета пространственного распределения плотности

ионного тока на расстояниях 0,1; 0,5; 1 м (сверху вниз) Интегральное значение тяги для двигателя КМ-5, рассчитанное для 3000000 моделируемых ионов на ПМА-тягомере равнялось 81,07мН с учетом 100% ионизации массового расхода 3,5мг/с. Измерения тяги двигателя в данном эксперименте показали значение ~75мН.

В третьем разделе представлены расчеты распыления струей двигателя

стенок вакуумной камеры 1. Для данной камеры характерна следующая особенность: около 60% ускоренных ионов попадают в днище камеры, остальные распределяются по цилиндрическим стенкам. Наиболее «распыляемой» частью камеры, за исключением днища, является область цилиндрических стенок на удалении 0,4...1,3 м от среза двигателя. В расчетах было получено, что скорость осаждения перепыленного материала со стенок камеры на поверхность I разрядной камеры двигателя может доходить до 4-Ю"1'м/с. Ранее (в гл.1) было приведены оценки предельных скоростей осаждения загрязнений, характерных для длительных испытаний двигателя СПД-100 в ОКБ «Факел», которые составляли -10"''м/с. Исходя из этого, для рассматриваемой камеры сделан вывод о необходимости осуществления специальных мероприятий для снижения загрязняющих потоков. В качестве таких действий было предложено использовать мишень-ионоприемник потока в торце камеры с призматическими выступами на поверхности и защиту стенок вставками из материала с низким коэффициентом распыления. Результаты моделирования показали, что наибольший эффект (снижение до 5 раз) достигается при использовании в зонах интенсивного распыления стенок вставок из графита (рис.5).

Рисунок 5. Плотность потока перепыления материала при различных покрытиях стенок камеры

В четвертом разделе главы помещены результаты моделирования потенциала в вакуумной камере с использованием магнитогидродинамической модели. Проанализировано влияние электрического поля на характер движения перезарядившихся ионов. Проведено сравнение плотности тока перезаряженных ионов в задней части камеры, рассчитанное с помощью модели, с результатами, которые были получены в камере 6 для двигателя СПД-70. Для данной камеры было подсчитано давление нейтральных частиц при коэффициенте выхода из насосов 0,55 и статическом давлении 7,23-10"3Па (рис.6).

Г

а) нержавеющая сталь

б) графит

1 184е-0:

1.114е-02 1 0"8МП

i 00Se-02 "36е-03

9 0J4e-03f

s jj.v-огЩ

"983 »-0з| " "б51е-03| ' _:80е-03

6.5"8е-0зИ

Рисунок 6. Результаты расчета давления нейтральных частиц Хе в камере 6 при работающем двигателе [1/м3]

Экспериментально была зафиксирована угловая зависимость распределения плотности ионного тока в струе рассматриваемого двигателя. По данной кривой был определен расход ксенона с использованием формулы:

m = f*/2(S(a)-i-Mion)da, (15)

где т - массовый расход газа, S(a) - площадь элемента сферы, ограниченная телесным углом da,j - плотность ионного тока на заданном угле, е - элементарный заряд, Mion- масса иона.

Расчет данного интеграла показал хорошую сходимость с экспериментальными данными (2,2мг/с - расчет и 2,3мг/с - измерение). Погрешность экспериментального распределения плотности ионного тока (по интегральному перерасчету массового расхода) составляла около -6%.

Расчет тяги с помощью ПМА-тягомера, для экспериментально измеренного распределения плотности тока, дал значение на уровне ~39,5мН. При испытаниях величина тяги данного двигателя изменялись в диапазоне 38,5...39,5мН в зависимости от давления остаточного газа в камере. Погрешность пространственного распределения массового потока в расчетах для данной камеры равняется -2%.

Пересчет массового расхода из углового распределения плотности ионного тока (рис. 7) по формуле (15) дал результат ~2,38мг/с (экспериментальные измерения - 2,3мг/с). При этом необходимо учесть, что экспериментальная зависимость (рис. 7) не была симметрична относительно центральной оси, что повлияло на точность моделирования, кроме того при измерениях некоторая часть первичных ионов «перезарядилась» и не фиксировалась зондом. Все это и послужило причиной отклонения результатов моделирования от экспериментальных данных.

-90 -80 -70 -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60 70 80 90 угол от оси, градус

Рисунок 7. Угловое распределение плотности ионного тока в струе СПД на расстоянии 0,5м (экспериментальные измерения - сплошная линия; расчет для системы из 30 ПМА - пунктир)

1.000000&4-017

5.411695е+01б

2.928<И5е+016

1.584893е+016

8.576959Н-015

4.641589е+015

' 2.5И886еИИ5

■ 1.359356е+015

7.356423еИН4 ; ■

' 3.981072(^014 Р

2.154435е+014 : 1

1.165914еМИ4

6.309573еИ)13 ■

3.414549е*013

1.847850(^013

!1.000000(М-013

Рисунок 8. Расчет распределения концентрации первичных ионов [1/м3]

На рис. 8 представлено распределение концентрации первичных ионов, рассчитанное по экспериментальным данным (для наглядности визуализации приведен нижний предел для значения давления, получаемого в расчете), а на рис.9 - распределение плотности тока первичных и перезарядившихся ионов вдоль оси камеры.

0 1 2 3 4 5 6

Расстояние вдоль оси камеры от заднего торцевого фланца, м

Рисунок 9. Результаты расчета распределения плотности тока ионов вдоль оси камеры 6 (первичные ионы - синяя линия; перезарядившиеся ионы на ПМА, повернутых к двигателю - коричневая линия; перезарядившиеся ионы на ПМА, повернутых от двигателя - оранжевая линия)

о

0,00001

0,0001 0,001 0,01 од

Расстояние вдоль оси камеры от заднего торцевого фланца, м

Рисунок 10. Результаты расчета распределения потенциала вдоль центральной оси камеры 6

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

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

Рассчитанное распределение потенциала на оси камеры показано на рис.10. При построении данного распределения использовалось одно экспериментально измеренное значение потенциала в точке, смещенной от центральной оси на 60 градусов и на удалении от двигателя 0,5м. Расчетное поле потенциалов имеет хорошую сходимость с экспериментальными данными (порядка 1В).

На рисунке 11 представлено сравнение с экспериментальными данными результатов расчета распределения ионов в струе по радиусу камеры. Данный график можно считать обобщающим тестом, подтверждающим выполнение в разработанной модели заявленных требований к погрешностям вычисления распределения ионов перезарядки. 1,6

£ <

1,4

¡5 0,8 и

0

1 0,6

0,2

с ) _ ..

> ь : \

\ ® V

^ Ч V*

О*» "О. о..

ОД 0,2 0,3 0,4 0,5 0,6 0,7 0,8

Удаление от центральной оси к стенке камеры, м

0,9

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

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

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

В заключении диссертации сформулированы следующие выводы:

1. Разработана трехмерная математическая модель и программные продукты для численного моделирования процессов: распределения нейтрального газа в объеме вакуумной камеры при работе СПД, движения первичных ионов, их взаимодействия с остаточной атмосферой и распространения перезарядившихся частиц с учетом влияния электрического поля.

2. Проведено сравнение результатов расчета с данными измерений отдельных параметров в вакуумных камерах различной формы и размеров. Получена сходимость результатов моделирования с экспериментальными данными на уровне 10...30%.

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

4. Разработана трехмерная математическая модель процессов распыления и пере-пыления материалов конструкций вакуумной камеры под действием ионной струн СПД, в которой впервые был применен алгоритм построения пространственного распределения распыленных частиц с помощью сплайн аппроксимации по заданным измерениям в двух плоскостях. С использованием данной модели выработаны рекомендации по уменьшению потока продуктов распыления для камеры конкретной геометрии.

Основные результаты опубликованы в следующих работах:

1. Торопов Г.П. Методика численного анализа распределения частиц в струе источника ионов // Сб. тез. докл. 9 Международной научно-технической конференции студентов и аспирантов «Радиоэлектроника электротехника и энергетика». -М.: МЭИ, 2003. - стр. 78,79.

2. Торопов Г.П., Хартов С.А. Моделирование процессов распространения струи плазменного двигателя в вакуумной камере // Сб. трудов Федеральной итоговой научно-технической конференции «Всероссийского конкурса на лучшие научные работы студентов по естественным, техническим наукам (проекты в области высоких технологий) и инновационным научно-образовательным проектам». -М.: МИЭМ. - 2004. - стр.37-39.

3. Торопов Г.П. Моделирование процессов распространения струи плазменного двигателя в вакуумной камере // Сб. тез. докл. 11 Международной научно-технической конференции студентов и аспирантов «Радиоэлектроника электротехника и энергетика». -М.: МЭИ, 2005. -стр. 112,113.

4. Toropov G.P. Khartov S.A. А 3D model calculating sputtering and depositing processes under Electric Propulsion thruster testing in a vacuum chamber IEPC-2005-074 // Proc. of 29th International Electric Propulsion Conference, 31 Oct. - 5 Nov. 2005, Princeton University, Princeton, USA.

5. Chervakov V., Khartov S., Toropov G. Numerical procedure permitting calculation of how a vacuum chamber influences Electric Propulsion thruster jet expansion. IEPC-2007-2 // Proc. of 30th International Electric Propulsion Conference, 17-20 Sept. 2007, Florence, Italy.

6. Toropov G.P. Khartov S.A. A 3D model of electron dynamics near the outlet section of SPT thruster IEPC-2007-343 // Proc. of 30th International Electric Propulsion Conference, 17-20 Sept. 2007, Florence, Italy.

7. Toropov G.P. A 3D boundary sheath analysis on Faraday probe measurements in a non-equilibrium plasma. IEPC-2007-4 //Proc. of 30th International Electric Propulsion Conference, 17-20 Sept. 2007, Florence, Italy.

8. Переславцев A.A., Торопов Г.П., Хартов C.A. Исследования процесса переноса электронов в прикатодной области разряда стационарного плазменного двигателя // Космонавтика и ракетостроение. - 2008. - Т.52. -№3. -с. 84-90.

9. Торопов Г.П., Хартов С.А. Математическая модель распространения струи плазменного ускорителя в объеме вакуумной камеры // Вестник Московского авиационного института. -2008. -Т.15.- №4.- с. 25-34.

ВВЕДЕНИЕ.

1. СОВРЕМЕННОЕ СОСТОЯНИЕ ИССЛЕДОВАНИЙ СТРУЙ СПД.

1.1 Особенности испытаний электроракетных двигателей в наземных условиях.

1.2 Состав остаточной атмосферы в вакуумнойкамере при работе СПД.

1.3 Эмпирические способы оценки параметров, влияющих на работу СПД в вакуумнойкамере.

1.4 Обзор экспериментальных исследований струй СПД в вакуумных камерах и в космосе.

1.5 Математические модели для расчета струй* СПД в, вакуумных камерах. Постановка задачи исследования.

2. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССОВ РАСПРОСТРАНЕНИЯ СТРУИ СТАЦИОНАРНОГО ПЛАЗМЕННОГО ДВИГАТЕЛЯ В ВАКУУМНОЙ КАМЕРЕ.

2.1 Моделирование уровня динамического вакуума.

2.2 Моделирование распределения концентрации первичных ионов.

2.3 Моделирование распределения концентрации ионов перезарядки.

2.4 Распыление и перепыление материала стенок камеры.

2.5 Оценка электрического поля.

2.6 Рассмотрение пристеночного слоя.

3. РЕАЛИЗАЦИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ РАСЧЕТА СТРУИ ПЛАЗМЕННОГО ДВИГАТЕЛЯ В ОБЪЁМЕ ВАКУМНОЙ КАМЕРЫ.

3.1 Описание компьютерной программы для расчета распространения струи и характеристик газа в вакуумной камере.

3.2 Метод определения характеристик с помощью поверхностного математического анализатора (ПМА).

3.2.1 Определение давления на заданную поверхность ПМА.

3.2.2 Определение тяги.

3.2.3 Определение плотности тока первичных ионов.

3.2.4 Определение плотности тока ионов перезарядки.

3.3 Анализ сечения упругого взаимодействия атомов ксенона.

3.4 Расчет коэффициентов Клаузинга для труб различного сечения.

3.5 Тестирование модели перепыления на полусфере.

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

4. РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ СТРУЙ СПД В РАЗЛИЧНЫХ КАМЕРАХ.

4.1 Расчет давления и концентрации в вакуумной камере при работе СПД.

4.2 Расчет распределения первичных ионов.

4.3 Расчеты распыления и перепыления.

4.4 Расчет распределения перезарядившихся ионов с учетом электрического поля.

Выводы.

В космической технике все более широкое применение в качестве тяговых модулей находят электроракетные двигатели (ЭРД), в частности одна из их разновидностей — стационарный плазменный двигатель (СПД). Наиважнейшим вопросом разработки и совершенствования СПД, как и любого ЭРД, является проблема наземных испытаний [1].

Одним из основных недостатков наземных испытаний ЭРД является невозможность в полном объеме имитировать условия окружающей среды, характерной для работы двигателя в космическом пространстве. Так, например, концентрация частиц остаточной атмосферы в вакуумных камерах на несколько порядков превосходит реальную концентрацию частиц в космическом пространстве [1,2,3,4,5]. Отличны также от наблюдаемых при эксплуатации двигателей пространственное распределение частиц и состав остаточной атмосферы испытательных стендов. Причиной этого является процесс взаимодействия ускоренных ионов со стенками вакуумных камер, в результате которого образуются обратные потоки атомов рабочего тела, нейтрализовавшегося на стенках, и распыленного с них материала. Эти потоки зависят от геометрии и размеров камеры. Верхний предел распыленных потоков во многом формирует требования для минимальных габаритов вакуумных камер при заданном уровне мощности испытуемого двигателя.

При распространении в объеме вакуумной камеры, струя ускоренных ионов из СПД взаимодействует с хаотично двигающимися нейтральными частицами остаточной атмосферы, образуя перезарядившиеся- ионы. Экспериментально подтверждено, что при этом измеряемый ионный ток в ядре струи уменьшается, а на периферии (на углах свыше 30 градусов от оси двигателя) значительно возрастает [6].

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

Комплексные экспериментальные исследования указанных явлений при испытаниях СПД достаточно трудоемки и не всегда дают однозначный результат. Вместе с тем отдельные процессы относительно хорошо изучены [8,9,10,11] (свободномолекулярное движение атомов; распыление материалов; перенос и осаждение распыленного материала; резонансная перезарядка и т.п.) и поддаются математической формализации. В последнее время были выполнены работы по моделированию поведения струи СПД в вакуумных камерах, которые провели группы исследователей: под руководством И. Бойда (I.D.Boyd) в Мичиганском университете (США), сотрудники фирмы Альта (Alta, Италия) и ряд др.[13,14]. В данных расчетах, в отдельных случаях, удалось добиться качественного соответствия результатов моделирования и экспериментальных данных. Однако модели, использованные в* этих работах, являются двухмерными и не позволяют оценивать влияние на распространение струи СПД реальной геометрии стенок камеры и расположения откачных средств.

В связи с этим тема данной работы, посвященной моделированию распространения струи СПД в объеме вакуумной камеры в трехмерном приближении, является актуальной.

В качестве цели работы была выбрана разработка методики моделирования процесса распространения струи СПД в объеме вакуумной камеры в трехмерном приближении, с учетом взаимодействия ионов с остаточным газом и стенками вакуумной камеры.

Для достижения поставленной цели в работе решались следующие задачи:

1. Разработка трехмерной математической модели распространения струи стационарного плазменного двигателя в вакуумной камере.

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

3. Создание математической модели для расчета распределения осаждения распыленных частиц материала стенок вакуумной камеры на поверхностях двигателя и вакуумной камеры.

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

Научная новизна работы состоит в том, что разработана и реализована в виде вычислительного комплекса трёхмерная инженерная модель распространения струи стационарного ^ плазменного двигателя в вакуумной камере. Предложен виртуальный датчик — поверхностный математический анализатор (ПМА), позволяющий проводить расчеты интегральных характеристик СПД (тяга, суммарный ионный ток), а также локальных значений параметров в струе двигателя (давление, плотность ионного тока, распределение частиц по скоростям).

Практическая значимость работы состоит в следующем:

- разработана методика моделирования процесса распространения струи СПД, позволяющая изучать влияние на динамику атомов и ионов таких факторов, как геометрия стенок и распределение нейтрального газа в объеме вакуумной камеры, а также определять характеристики потоков ионов и нейтральных атомов внутри вакуумной камеры в процессе работы двигателя;

- проведено моделирование для конкретных геометрий вакуумных камер динамического распределения нейтральных частиц до и после включения двигателя;

- разработана методика численного определения локальных характеристик в заданных точках струи СПД;

- разработана методика построения распределения перезарядившихся частиц в объеме вакуумной камеры;

- разработана методика моделирования процессов распыления стенок вакуумной камеры и перепыления распыленного материала; - разработана методика расчета интегральных параметров двигателя (тяга, суммарный ионный ток).

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

На защиту выносятся:

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

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

3. Результаты моделирования динамики атомов и ионов в вакуумных камерах различной конфигурации.

4. Трехмерная инженерная модель оценки процессов распыления и перепыления материала стенок вакуумной камеры.

Выводы

1. Разработана трехмерная математическая модель и программные продукты для численного моделирования процессов: распределения нейтрального газа в объеме вакуумной камеры при работе СПД, движения первичных ионов, их взаимодействия с остаточной атмосферой и распространения перезарядившихся частиц с учетом влияния электрического поля.

2. Проведено сравнение результатов расчета с данными измерений отдельных параметров в вакуумных, камерах различной формы и размеров. Получена сходимость результатов моделирования с экспериментальными данными на уровне 10. .30%.

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

4. Разработана трехмерная математическая модель процессов распыления и перепыления материалов конструкций вакуумной камеры под действием ионной струи СПД, в которой впервые был применен алгоритм построения пространственного распределения распыленных частиц с помощью сплайн аппроксимации по заданным измерениям в двух плоскостях. С использованием данной модели выработаны рекомендации по уменьшению потока продуктов распыления для камеры конкретной геометрии.

1. Biagioni L., Kim V., Nicolini D. et.al. Basic issues in electric propulsion testing and the need for international standards. IEPC-2003-230 // Proc. of 28th International Electric Propulsion Conference, 17-21 March 2003, Toulouse, France.

2. Korsun A.G., Tverdokhlebova E.M., Gabdullin F.F. The distinction between the EP plume expansion in space and in vacuum chamber. IEPC-2005-073 // Proc. of 29th International Electric Propulsion Conference, 31 Oct. 4 Nov. 2005, Princeton University, USA.

3. Nakles M.R., Hargus Jr.W. A background pressure effects on internal and near-field ion velocity distribution of the BHT-200-X3 Hall thruster // Proc. of 55th JANNAF Propulsion Meeting, 8- 12 Dec. 2008, Orlando, Florida, USA.

4. Hofer R.R., Peterson P.Y., Gallimore A.D. Characterizing vacuum facility backpressure effects on the performance of a Hall thruster. IEPG-01-45 // Proc. of 27th International Electric Propulsion Conference, 15-19 Oct. 2001, Pasadena, .CA, USA.

5. Hofer R. R., Walker M. L.R., Gallimore A. D. A comparison of nude and colliitmated faraday probes for use with Hall thrusters. IEPC-01-20 // Proc. of 27ш International Electric Propulsion Conference, 15-19 Oct. 2001, Pasadena, CA, USA.

6. Горшков. О.А., Муравлёв B.A., Шагайда А.А. Холловские и. ионные плазменные двигатели для космических аппаратов / под ред. академика РАН А.С. Коротеева. — М.: Машиностроение, 2008.

7. Проблемы прикладной физики. Распыление твердых тел ионной бомбардировкой. Физическое распыление одноэлементных твердых тел / под ред. Р. Бе-риша; пер. с англ. под ред. В.А. Молчанова. М., 1984.

8. Кошмаров, Ю.А., Рыжов Ю.А. Прикладная динамика разреженного газа -М.: Машиностроение, 1977.

9. Шагайда, А. А. Метод численного моделирования течений разреженного газа и его применение для расчета электрофизических устройств : автореферат дис. кандидата физико-математических наук : 01.02.05/Моск. физ.-техн. ин-т (гос. ун-т), 2000.

10. Bugrova A. I., Morozov A. I. Influence of vacuum conditions on the SPT operation. IEPC-95-46 // Proc. of 24th International Electric Propulsion Conference, 19-23 Sept. 1995, Moscow, Russia.

11. Приданников С.Ю. Исследование характеристик стационарных плазменных двигателей при длительной работе: дис. канд. тех. наук: 05.07.05 / При-данников Сергей Юрьевич. Калининград (обл.), 2003. - 156с.

12. Boyd I. D., Sun Q., Cai Ch., Tatum К. E. Particle simulation of Hall thruster plumes in the 12V vacuum chamber. IEPC-2005-62 // Proc. of 29th International Electric Propulsion Conference, Oct. 31 — Nov. 4 2005, Princeton University, USA.

13. King L.B. Transport-property and mass spectral measurements in the plasma exhaust plume of a hall-effect space propulsion system: Ph.D. thesis (Aerospace Engineering), Michigan University, 1998.

14. Boerner J.J., Boyd I. D. Numerical simulation of Faraday probe measurementsiKin a multi-component non-equilibrium plasma. IEPC-2005-085 // Proc. of 29 International Electric Propulsion Conference, 31 Oct. — 4 Nov. 2005, Princeton University, USA.

15. Miller J. S., Pullins S. H., Levandier D. J., et.al. Xenon charge exchange cross sections for electrostatic thruster models // Journal of Applied Physics. — 2002. Vol 91, No. 3.-pp. 984-991.

16. Gulczinski F.S., Gallimore A.D., Williams, G.J., et.al. Laser induced fluorescence measurements of ion velocities in the plume of a Hall effect thruster AIAA99.2424// Proc. of 35th Joint Propulsion Conference, 20-23 June, 1999, Los Angeles, CA, USA.

17. Асхабов C.H., Бургасов М.П., Фишгойт B.B., и др. Исследование струи стационарного плазменного двигателя с замкнутым дрейфом электронов // Физика плазмы. 1981. - том 7, вып. 31. — с. 225-230.

18. Manzella D.H., Sankovic J.M. Hall thruster ion beam characterization. AIAA-95-2927 // Proc. of 31st Joint Propulsion Conference, 10-12 July 1995, San Diego, CA, USA.

19. Арцимович JLA., Морозов А.И., Снарский P.K., и др. Разработка стационарного плазменного двигателя (СПД) и его испытания на ИСЗ «Метеор» // Космические исследования. 1974. — том 12, № 3. — с. 451-468.

20. KimV., Kozlov V., Popov G., Skrylnikov A. Plasma parameter distribution determination in SPT-70 plume. IEPC-2003-107 // Proc. of 28th International Electric Propulsion Conference, 17-21March, 2003, Toulouse, France.

21. Randolph Т., Kim V., Kaufman H., K. et.al. Facility effects in stationary plasma thruster testing. IEPC-93-093 // Proc. of 23rd International Electric Propulsion Conference, 19-23 Sept., 1993, Seattle, USA.

22. Kozubsky K., Randolph Т., Kim V., et.al. Effect of background nitrogen and oxygen on insulator erosion in the SPT-100. IEPC-93-092 // Proc. of 23rd International Electric Propulsion Conference, 19-23 Sept., 1993, Seattle, USA.

23. Pagnon D., Touzeau M., Lasgorceix P. Control of the ceramic erosion by optical emission spectroscopy: parametric studies of the SPT100-MML. AIAA 2004-3773 // Proc. of 40th Joint Propulsion Conference, 11-14 July, 2004. Fort Lauderdale, Fl., USA.

24. Mazouffre S., Pagnon D., Bonnet J. Two ways to evaluate the Xe+ ion flow velocity in a Hall effect thruster-LIF spectroscopy and Fabry-Pdrot interferometry.

25. AIAA 2004-3949 // Proc. of 40th Joint Propulsion Conference, 11-14 July, 2004. Fort Lauderdale, Fl., USA.

26. Hargus Jr. W., Nakles M. R. A low power Hall effect thruster discharge plasma oscillation characteristics // Proc. of 55th JANNAF Propulsion Meeting, 8-12 Dec., 2008. Orlando, Fl., USA.

27. Gray H.L., Kamath U.P. INTELSAT 10 plasma propulsion, system initial flight operations. AIAA 2005-3672 // Proc. of 41st Joint Propulsion Conference, 10-13 July, 2005. Tucson, AZ, USA.

28. Gray H., Provost S., Glogowski M., Dermaire A. INMARSAT 4F1 plasma propulsion system initial flight operations. IEPC-2005-082// Proc. of 29th International Electric Propulsion Conference, 31' Oct. 4 Nov. 2005; Princeton University, USA.

29. Demiare A. Gray H. L. Plasma propulsion system functional chain first three years in orbit on Eurostar 3000. IEPC-2007-060 // Proc. of 30th International Electric Propulsion Conference, 17-20 Sept. 2007, Florence, Italy.

30. Lichtin D. A. An Overview of Electric Propulsion activities in US industry-2005. AIAA 2005-3532// Proc. of 41st Joint Propulsion Conference, 10-13 July, 2005. -Tucson, AZ, USA.

31. Semenkin A., Kim V., Gorshkov O., Jankovsky R. Development of electric propulsion standards current status and further activity. IEPC-01-070 // Proc. of 27th International Electric Propulsion Conference, 15-19 Oct. 2001 Pasadena, CA, USA.

32. WalkerM. L.R., Gallimore A.D., Cai C., Boyd I.D. Pressure map of a facility as a function of flow rate to study facility effects. AIAA 2002-3815// Proc.of 38th Joint Propulsion Conference, 7-10 July, 2002. Indianapolis, IN, USA.

33. Andrenucci M., Biagioni L., Passaro A. PIC/DSMC models for Hall effect thus-ter plumes: present status and ways forward. AIAA 2002-4254 // Proc.of 38th Joint Propulsion Conference, 7-10 July, 2002. Indianapolis, IN, USA.

34. Физико-химические процессы в газовой динамике / под ред. Г.Г. Черного, С.А. Лосева. М.: Научный мир, 2007.

35. Горелова Г.В., Кацко И.А. Теория вероятностей и математическая статистика в примерах и задачах с применением EXCEL. Ростов-на-Дону: Феникс, 2002.

36. Ивановский Г.Ф., Петров В. И. Ионно-плазменная обработка материалов М:Радио и связь. 1986.

37. Семенов А.А., Шкарбан И.И. Распыление потоками ионов поверхностей элементов конструкций ионно-плазменных источников // Ракетно-космическая техника, вып. 3 (131), М.: НИИТП, 1991. с. 42-53.

38. Бронштейн И.Н., Семендяев К.А. Справочник по математике для инженеров и учащихся втузов. — М.: Наукал 1981.

39. Морозов А.И., Есипчук Ю.В., Тилинин Г.Н., и др. Экспериментальное исследование ускорителя с замкнутым дрейфом; электронов // ЖТФ. 1972. -том 42, № 1. — с. 54-63.

40. Boyd I.D., Yim J. Т. Hall thruster plume simulation using a detailed hybrid model. AIAA 2004-3952 // Proc. of 40th Joint Propulsion Conference, 11-14 July, 2004. Fort Lauderdale, Fl., USA.

41. Kim V., Bishaev A., Lazourenko A., Auweter-Kurtz M. 3-Dimensional simula1.tion of plasma dynamics in SPT. IEPC-01-340 // Proc. of 27ш International Electric Propulsion Conference, 15-19 Oct. 2001, Pasadena, CA, USA.

42. Bishaev A.M., Kalashnikov V.K., Kim V., Shavykina A.V. Numerical modeling of the propagation of a plasma jet produced by a stationary plasma thruster in a low-pressure gas // Plasma, Physics Report. 1998. - v.24, No 11. - pp. 989-995.

43. Bishaev A.M., Kim V., Lazourenko A.V. Simulation of plasma dynamics in SPT // Proc. of 3rd International Conference on Spacecraft Propulsion, 10-13 Oct. 2000, Cannes, France / ESA SP-465, December 2000. pp. 385-391.

44. Bird, G. A. Molecular gas dynamics and the direct simulation of gas flows. Clarendon Press, 1994.

45. Abrahamson A. Born-Mayer-Type Interatomic potential for neutral ground state atoms with Z=2 to Z=105 // Physical Review. Feb. 1968.

46. Розанов JI.H. Вакуумная техника M.: Высшая школа. - 1990.1. Напряжение на зонде,В

47. Рисунок Б.5. Вольтамперная характеристика зонда для угла 30 градусов

48. Рисунок Б.6. Вольтамперная характеристика зонда для угла 40 градусов

49. Обработка волътамперной характеристики Использовалась следующая формула для получения из вольтамперной характеристики распределения по энергиям 54.:йЕ1. Ш1+11. Ш;1. М+1 1. Б.1)

50. На графиках (рис. Б.7.Б.12) представлены результаты обработки для соответствующих углов двигателя №2.

51. Рисунок Б. 11. Распределение плотности вероятности ионов по энергиям в струе на угле 30 градусов от оси двигателя

52. Рисунок Б. 12. Распределение плотности вероятности ионов по энергиям в струе на угле 40 градусов от оси двигателя

53. Рисунок Б. 13. Распределение плотности вероятности ионов по энергиям в струе от оси двигателя ( а 0 градусов; б - 7,5; в — 15; г - 22,5)

54. Г. 4. Стенд исследовательского центра Centrospazio35.

55. Рисунок Г.6. Трехмерная модель вакуумной камеры

56. Г. 5. Вакуумная камера 12V в Arnold Engineering Development Center. 3D модель, построенная по приведенным в статье 13. размерам камеры

57. Рисунок Г.7.Трехмерная модель камеры

58. Рисунок Г. 8. Криогенные панели верхнего и нижнего насосовв данной камере.1. Г.6. СтендУ2В в МАИ

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

60. Рисунок Г.9. Чертеж камеры 23.