автореферат диссертации по авиационной и ракетно-космической технике, 05.07.05, диссертация на тему:Исследование характеристик струй стационарных плазменных двигателей (СПД) при повышенных разрядных напряжениях

У046

5584

Архипов Алексей Сергеевич

Исследование характеристик струй стационарных плазменных двигателей (СПД) при повышенных разрядных напряжениях

Специальность 05.07.05 - «Тепловые, электроракетные двигатели и

энергоустановки летательных аппаратов»

Автореферат диссертации на соискание ученой степени кандидата технических наук

Москва-2010

004615584

Работа выполнена в Научно-исследовательском институте прикладной механики и электродинамики Государственного образовательного учреждения высшего профессионального образования «Московский авиационный институт (государственный технический университет)» (НИИ ПМЭ)

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

доктор технических наук, старший научный сотрудник Ким Владимир

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

доктор технических наук, профессор

Ходненко Владимир Павлович

кандидат технических наук Семенкин Александр Вениаминович

Ведущая организация — Федеральное государственное унитарное предприятие «Опытное конструкторское бюро «Факел» (ФГУП «ОКБ «Факел»), г. Калининград

Защита состоится « -о^^^оО^/Р 20^ года, в /•$ часов на заседании диссертационного советаЭД 212.12^.08 при Московском авиационном институте (государственном техническом университете) по адресу: 125993, г. Москва, А-80, ГСП-3, Волоколамское шоссе, д. 4.

С диссертацией можно ознакомиться в библиотеке Московского авиационного института (государственного технического университета).

Автореферат разослан » 20 г.

Ученый секретарь диссертационного совета

Д 212.125.08 д.т.н., профессор = Ю.В.Зуев

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

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

- межорбитальная транспортировка;

-поддержание длительной работоспособности КА на геостационарной орбите (ГСО) со сроками активного существования (САС) 12-15 лет;

- межпланетные перелеты и полеты в дальний космос.

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

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

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

С учетом изложенного в качестве целей диссертационной работы были выбраны следующие:

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

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

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

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

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

1. Провести анализ и исследование характеристик струи моделей современных двигателей типа СПД-100 и СПД-140 при их работе с повышенными разрядными напряжениями.

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

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

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

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

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

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

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

- результаты экспериментального определения характеристик струй двигателей, работающих при повышенных разрядных напряжениях, используются при оценке механического, эрозионного и загрязняющего воздействий двигателей типа СПД-ЮОД и СПД-140Д в ОКР «Двина ТМ» и «Луч»;

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

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

Основные положения, выносимые на защиту

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

2. Положение о решающем влиянии температуры электронов на расходимость струи СПД.

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

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

Апробация работы и научные публикации

Результаты диссертационной работы докладывались и обсуждались на научно-технических семинарах в Научно-исследовательском институте прикладной механики и электродинамики (НИИ ПМЭ), 5-ой, 6-ой, 7-ой и 8-ой Международных конференциях по неравновесным процессам в соплах и струях, 14-ой, 15-ой и 16-ой Международных конференциях по вычислительной механике и современным прикладным программным системам, Международной конференции «Фундаментальные проблемы высокоскоростных течений», 3-ей, 6-ой и 7-ой Международных конференциях «Авиация и космонавтика», 2-ой Европейской конференции по аэрокосмическим наукам, 31-ой Международной конференции по электроракетным двигателям. Основное содержание и результаты диссертационной работы изложены в 5-ти публикациях, список которых приведен в конце автореферата:

Вклад автора

Автор принимал непосредственное участие в экспериментах, в ходе которых определялись угловые распределения плотности ионного тока и энергетические спектры для моделей двигателей СПД-МО-О и СПД-140М,

5

занимался обработкой и физическим анализом полученных данных. Автором была разработана методика расчета интегральных параметров струи СПД с учетом доли двухзарядных ионов и с ее использованием был проведен анализ зависимости коэффициентов, отражающих влияние на потери тяги и тягового к.п.д. (неполноты использования рабочего тела в ускорительном канале и неполноты использования приложенного разрядного напряжения на ускорение ионов, расходимости струи и разброса ионов по скоростям) от разрядного напряжения. Автором была разработана методика и экспериментально определены параметры квазирадиальных потоков ионов вблизи среза ускорительного канала опытного образца двигателя СПД-ЮОД, работающего на режиме с высоким удельным импульсом тяги, а совместно с в.н.с. A.M. Бишаевым разработана методика расчетной оценки указанных потоков и проведены расчеты, подтверждающие механизм их формирования.

Структура и объем диссертации

Работа представляет собой рукопись объемом 150 страниц печатного текста, включая 85 рисунков, 7 таблиц, а также список цитируемой литературы, насчитывающий 90 наименований. Работа состоит из введения, четырех глав, заключения и списка цитируемой литературы.

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

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

Во второй главе описаны методика и результаты экспериментального исследования характеристик струй моделей двигателей СПД-100В, СПД-ЮО-О и СПД-140М, разработанных в НИИ ПМЭ и работающих при повышенных разрядных напряжениях, приводится анализ полученных результатов. Исследования проводились на стенде НИИ ПМЭ У-2В-1. Давление в камере при работе двигателя с расходами ксенона (2-5) мг/с составляло (2.5-5)х10~5 мм.рт.ст. В качестве основного средства диагностики струи СПД использовался плоский электростатический зонд-энергоанализатор, который устанавливался внутри вакуумной камеры на штанге, перемещающей его по

6

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

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

-60 -40 -20 О 20 40 60

Рис. 1 Относительные распределения плотности тока ускоренных ионов в струе модели двигателя СПД-ЮО-О для различных режимов работы

Рис. 2 Относительные распределения плотности тока ускоренных ионов в струе модели двигателя СПД-140М для различных режимов работы.

Анализ угловых и относительных угловых распределений плотности ионного тока, полученных для моделей двигателей СПД-ЮО-О и СПД-140М, работающих на различных режимах, выявил следующие тенденции:

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

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

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

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

О 200 400 600

-и^кю в —а— 1м=5оо в —а

(а)

и>,в

еоо юоо 1200

№700 В Ш=500 В

?ПП ЛПП ЛОО япл -1№300В -в—№500 В —I

иа.В

100(1 1?00 -№«00 В

(В)

(Нс/сШа,

800 1000 1200 ■и(1=7оо в о иа-эоо в

йПП япо

-Ш=600В -

и«, в

1ППО 15П0

к— 1М=»00В

(б) (Г)

Рис. 3 Функции распределения ионов по энергии для моделей двигателей СПД-ЮО-О (а,б) и СПД-140М (в,г), работающих при расходах через анод та = 2.4 мг/с и та = 4.5 мг/с соответственно и различных разрядных напряжениях (а,в - 0 градусов от оси СПД, б,г - 10 градусов от оси СПД)

В результате исследования функций распределения ионов по энергии для моделей двигателей СПД-ЮО-О (см. Рис. 3 а,б) и СПД-140М (см. Рис. 3 в,г), было показано, что:

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

-отмеченные закономерности сохраняются при увеличении разрядного напряжения.

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

Еср, ЭВ

900

-40 -20 0 =300 В, т=2.4 мг/с -1=700 В, т=2.4 мг/с -1=300 В, т=4.5 мг/с

20 40 60 -Ш=500 В, т=2.4 мг/с -Ш=900В, т=2.4 мг/с

•60 -40 -20 О ~в-~ ш=300 В, тМ.5 мг/с --А— Ш=800 В, т=4.5 мг/с --*— Ш=600 В, т=7.02 мг/с

20 40 60 —иЧ=500 В, т=4.5 мг/с -Ш=800 В, т=5.51 мг/с

(Ю (б)

Рис. 4 Распределения средней энергии ионов в струе моделей двигателей СПД-ЮО-О (а) и СПД-140М (б) по углу отклонения направлений их вылета от оси двигателя

Была выявлена также разница в зависимости средней энергии ионов от угла р при одном и том же разрядном напряжении, но при разных расходах рабочего газа (см. Рис. 4 а), а именно: при разрядном напряжении 300 В и расходе рабочего газа 2.4 мг/с средняя энергия ионов на периферии струи значительно выше, а в приосевой области - заметно ниже, чем при расходе 4.5 мг/с. Отмеченная разница определяется изменением функции распределения ионов по энергии в соответствующих частях струи.

Отмечено, что в приосевой части потока при -10° </?<+10° с увеличением разрядного напряжения ширина распределения мало изменяется по сравнению с увеличением средней энергии или разрядного напряжения. Сказанное подтверждается зависимостью полуширины распределения от разрядного напряжения и отношения этой полуширины к разрядному напряжению (см. Рис. 5 а,б).

700 900

Диа, р*10*

йиа/Ш, р=10*

300 400 500 -*-Д1)а, р=0' * ¿иа/иа. р-0'

-Диа, р=10' -ЛиаЛМ р»10*

(а) (б)

Рис. 5 Зависимости полуширины распределения ионов по энергии в приосевой части струи для моделей двигателей СПД-ЮО-О (а) и СПД-140М (б)

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

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

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

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

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

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

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

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

движущихся по одному и тому же направлению, одинаковы.

При этих допущениях система основных уравнений для определения доли двухзарядных ионов сводится к следующей:

1,М

1)./,= Ыя1]\{К,Р)ьхъ№Р\ 2).Л+//2=1; 3). -г;

о тае{\ + цг)

м

4). <К>=у ]<К> (/?)/,№№/?<//?; 5). 81 !р = цм < У„ >.

А о

о

где /( - полный ионный ток ускоренных ионов в струе, /л, и - доля однозарядных и двухзарядных ионов соответственно, г]т - коэффициент использования рабочего газа в ускорительном канале, та, М и е -соответственно массовый расход через анод, масса иона и заряд электрона, < У1г > - осевая составляющая скорости ионов, g и / - соответственно ускорение свободного падения тел на Земле и измеренное значение удельного импульса тяги.

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

Поиск решения приведенной системы осуществлялся методом последовательных приближений. При практической реализации решения дополнительно рассчитывались полуугол расходимости потока ускоренных ионов /?095, в пределах которого распространяется 95 % ионов, средняя энергия ионов <е(>, а также коэффициенты отражающие влияние на величину тягового к.п.д. двигателя угловой расходимости струи гц, разброса ионов по скоростям Цу и неполноты использования разрядного напряжения на ускорение ионов Г!с.

Расчет проводился по разработанному автором алгоритму, на основе которого в среде МаСаЬ была написана специальная расчетная программа для решения системы уравнений, по которым определялись все вышеуказанные параметры. Результаты расчетов представлены в таблицах 1 и 2.

Таблица 1

Интегральные параметры моделей двигателей СПД-ЮО-О и СПД-140М (расход в таблице пересчитан с учетом поправки на давление в вакуумной камере)_

мг/с В Л. А "о Вт т, мН СТ> Вт/мН Пт /„А Д>.95 > град

СПД-ЮО-О

2.53 307 2.12 651 37.0 17.6 1492 0.42 1.71 0.81 45.90

2.53 492 2.24 1102 51.6 21.4 2082 0.48 1.85 0.83 44.50

2.53 695 2.44 1696 65.0 26.1 2625 0.49 1.96 0.80 44.90

2.53 886 2.74 2428 75.8 32.0 3060 0.47 2.18 0.79 46.40

СПД-140-М

4.95 299 4.45 1328 78.4 16.9 1613 0.47 3.66 0.82 48.50

4.95 510 4.51 2296 110.3 20.8 2269 0.53 4.05 0.90 45.30

4.95 802 4.72 3788 148.6 25.5 3057 0.59 4.32 0.92 46.20

Таблица 2

Характеристики струи моделей двигателей СПД-ЮО-О и СПД-140М (расход в таблице пересчитан с учетом поправки на давление в вакуумной камере)

мг/с В <Б(> , эВ км/с <У1:>, км/с Иг 1* П, Л у п. Чт 1тР

СПД-ЮО-О

2.53 307 218 14.92 16.40 0.022 0.91 0.86 0.96 0.71 0.42 0.43

2.53 492 381 20.82 22.06 0.065 0.94 0.86 0.96 0.77 0.48 0.50

2.53 695 548 26.25 27.24 0.105 0.96 0.88 0.95 0.79 0.49 0.51

2.53 886 711 30.60 30.60 0.184 1.00 0.81 0.94 0.80 0.47 0.49

СПД-140-М

4.95 299 231 16.13 17.34 0.092 0.93 0.87 0.94 0.77 0.47 0.48

4.95 510 382 22.69 23.39 0.159 0.97 0.90 0.95 0.75 0.53 0.55

4.95 802 641 30.57 30.58 0.201 1.00 0.89 0.93 0.80 0.59 0.61

Расчет доли двухзарядных ионов в ускорительном канале дал достаточно реалистичные значения и соответствует ожидаемым тенденциям, а именно:

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

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

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

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

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

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

-коэффициент г]р, учитывающий потери энергии из-за расходимости

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

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

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

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

11(1, в

(а)

(б)

Рис. 6 Зависимости тягового к.п.д., рассчитанного по результатам измерения тяги щ и характеристик струи Т1Тр от разрядного напряжения для моделей двигателей СПД-ЮО-О (а) и СПД-140М (б)

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

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

V

А*-77-. (3)

<ауш>п,

где Уа - скорость атомов, < <т,Ке >= к, - коэффициент скорости ионизации, осредненный по функции распределения ионов, пе - концентрация электронов.

При невысоком и постоянном разрядном напряжении и понижении расхода главным фактором увеличения очевидно, является уменьшение концентрации ионов и электронов в ускорительном канапе. Согласно результатам измерений в модели двигателя СПД-ЮОВ температура электронов возрастает с увеличением разрядного напряжения, вместе с этим коэффициент скорости ионизации к) возрастает быстрее, чем падает концентрация плазмы. Таким образом, с увеличением разрядного напряжения должно происходить уменьшение длины свободного пробега атомов до их ионизации и увеличение вероятности их ионизации. И это может служить объяснением увеличения суммарного ионного тока в струе двигателя, наблюдавшегося при испытаниях моделей двигателей СПД-ЮО-О и СПД-140М при увеличении разрядного напряжения с 300 В до 700 В и выше.

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

кТ

в радиальном направлении у выхода из канала можно оценить как Д<рг « —-.

е

Тогда характерное значение радиальной скорости ионов в струе можно представить в виде:

. 2е 2кТ

А/р'м=)11Г>

(4)

а характерное значение продольной составляющаеи скорости ионов можно определить как:

М

(5)

Таким образом, в первом приближении расходимость ионного потока

может быть оценена с использованием tgQ =

К enPd

, т.е.

в ж arctg

кТ,

(6)

С учетом того, как введены значения параметров в соотношениях (4)-(6), расчетное значение угла в сравнивалось с усредненным по потоку значением угла отклонения от оси, в качестве которого был выбран полуугол Д, 5 конуса с вершиной, расположенной в выходной плоскости и на оси двигателя, в пределах которого, распространяется 50 % потока ускоренных ионов. Для расчета угла в бралось максимальное значение температуры электронов Г,, измеренное с помощью пристеночных зондов в канале двигателя. Расчеты были проведены для модели двигателя СПД-ЮОВ.

Результаты расчета 9 и Д,5 представлены на Рис. 7. Видно, что полученные зависимости близки по характеру и по величине. Учитывая это, можно объяснить характер изменения расходимости струи с изменением разрядного напряжения характером изменения максимальной величины температуры электронов в ускорительном канале от разрядного напряжения.

300 350 400 450 500 550 600 650 700

Рис. 7 Зависимости угла в и полуугла ßos от разрядного напряжения для модели двигателя СПД-ЮОВ, работающей при расходе через анод 2.5 мг/с

В четвертой главе приводятся результаты расчетного и экспериментального исследований распределений параметров плазмы в ближней зоне струи СПД. Для расчетов была использована физико-математическая модель разработанная при лидирующей роли A.M. Бишаева, в которой для описания поведения струи вводятся функции распределения ионов и нейтралов /(*,#), g(x,w)

15

соответственно. Подробное описание используемой модели представлено в работе [1] из списка публикаций, приведенного в конце автореферата.

С помощью указанной физико-математической модели была проведена серия расчетов с целью исследования влияния увеличения разрядного напряжения на распределение макропараметров струи. Расчеты проводились для двигателя типа СПД-100. Его габариты были заданы в виде куба с размером ребра 150 мм. Внутренний и внешний радиусы выходного отверстия ускорительного канала - 35 мм и 50 мм, соответственно. Размеры счетной области в абсолютных величинах составляли 350x350x900 мм. Распределения осевой и продольной составляющих плотности ионного тока в отверстии брались из полученных ранее экспериментальных данных. Температура ионов на выходе из отверстия задавалась равной максимальной температуре электронов в канале, измеренной пристеночными зондами, для сопоставимых по расходу и напряжению режимов. Аналогично задавались параметры и для нейтральной компоненты.

В результате проведенного численного исследования были получены распределения основных параметров потока ионов и нейтралов соответствующие режимам работы двигателя типа СПД-100 с расходом газа через анод 2.4 мг/с и разрядными напряжениями 300 В, 500 В, 700 В и 900 В. Изоконцентрали ионов в продольном сечении для разрядных напряжений 300 В и 900 В показаны на Рис 8 а,б.

(а) (б)

Рис. 8 Линии уровня равной плотности ионов в продольном сечении и структура радиальных потоков для двигателя типа СПД-100, работающего на режимах с расходом через анод 2.4 мг/с и разрядными напряжениями 300 В (а) и 900 В (б)

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

- происходит сужение потока ионов к оси симметрии двигателя;

- сжимаются распределения концентрации ионов;

-значения концентрации, а, следовательно, и плотности тока, в ядре потока возрастают примерно пропорционально разрядному напряжению.

В целом указанные закономерности согласуются с обнаруженными и проанализированными в Главе 2 тенденциями изменения параметров струи СПД с ростом разрядного напряжения. Следует также отметить, что во всех расчетных сериях были получены радиально направленные и «обратные» ионные потоки, структура которых приведена на Рис 8 а,б. Как видно из

приведенных рисунков, указанные

Ji, мА/см1

плоскость среза ускорительного канала СПД

.....О.....рсчет, Ud=300 В, т=2.4 ыг/с

—в— рсчег, Ud=500 В, т=2.4 мг/с А рсчет, Ud=700 В, т=2.4 мг/с —в—рсчвг, Ud=800 В, т=2.4 мг/с —М—рсчег, Ud=900 В, 1Ш=2.4 мг/с —эксперимент, Ud=800 В, ld=2-6A

потоки могут оказывать распыляющее воздействие на катод-нейтрализатор, который

располагается в зоне

распространения этих потоков.

На Рис. 9 приведены зависимости распределения

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

экспериментальные данные,

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

следующем:

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

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

Были также рассчитаны обратные ионные потоки на область лежащую в плоскости среза ускорительного канала. Полученные в результате численного эксперимента данные удалось сравнить с измерениями обратных ионных потоков, которые проводились для модели двигателя СПД-100В во время 100-часовых ресурсных испытаний на режиме работы с разрядным напряжением 750 В и расходом газа через анод 2.4 мг/с. Было достигнуто

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

удовлетворительное качественное соответствие расчетных и экспериментальных данных.

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

В Главе 1 было показано, что во время работы двигателя квазирадиальные потоки ионов могут вызывать эрозию элементов двигателя и КА, расположенных вблизи среза ускорительного канала двигателя. Для оценки значимости этих потоков автором помимо приведенных выше расчетов была разработана методика их экспериментального определения с помощью зондов Ленгмюра. В силу того, что исследуемые потоки ионов имеют относительно невысокие значения, необходимо было предпринимать специальные меры с тем, чтобы увеличить сигнал, поступающий с коллектора зонда. С этой целью были разработаны и созданы зонды с увеличенной собирающей поверхностью.

Указанные потоки определялись для опытного образца двигателя СПД-100Д, производства ОКБ «Факел». Зонды с помощью крепежей, выполненных из немагнитных материалов, устанавливались на наружный магнитный полюс двигателя по азимуту через каждые 90 градусов таким образом, чтобы их ось пересекала ось симметрии двигателя и была бы ей перпендикулярна, а коллекторы были бы ориентированы к струе (см. Рис. 10).

Визуальный осмотр отработавшего двигателя показал, что зона с наблюдаемой эрозией поджигающих электродов катодов начинается на уровне плоскости среза ускорительного канала и заканчивается в 4-5 мм за плоскостью

-зонды 0 мм от среза -зонды 6.3 мм от среза

135 180

Рис. 11 Распределение радиальных ионных потоков по азимуту

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

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

- для зондов №№ 2,4,5 располагающихся на расстоянии 6.3 мм от среза ускорительного канала, значения потоков ионов на коллекторы по азимуту изменяются не значительно (см. Рис. 11);

- плотность ионного тока на зонд № 1, расположенный между катодами в плоскости среза

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

- с увеличением расстояния от среза ускорительного канала в продольном направлении (вниз по потоку) плотность ионного тока пропорционально возрастает примерно до расстояния 30 мм, а далее практически не изменяется (см. Рис. 9).

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

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

максвелловская функция распределения выходящих из катода нейтралов и соответственно все ее параметры.

* ■

¡-ИВ

в

А В

о

(а)

&

!' ©

§ 9

а м

ВП

ЯК Шт кШшшН

ЕВ л...

Двнии

(б)

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

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

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

1. Проведено экспериментальное исследование характеристик струй отечественных двигателей типа СПД-100 и СПД-140 и выявлены следующие закономерности:

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

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

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

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

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

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

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

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

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

-полуугол расходимости струи несколько уменьшается с увеличением разрядного напряжения до (500-700) В, а затем слабо возрастает, при этом характер его зависимости от разрядного напряжения может быть объяснен зависимостью радиальной составляющей скорости ионов от температуры электронов и характером ее изменения в ускорительном канале при увеличении разрядного напряжения;

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

-потери тяговой эффективности из-за расходимости струи составляют 10-15 %, а разброс ионов по скоростям дает наименьшие (~5 %) потери в тяговой эффективности двигателя.

3. Совместно с другими сотрудниками НИИ ПМЭ проведены измерения характеристик струй опытных образцов двигателей типа СПД-100Д и СПД-140Д на режимах работы с высоким удельным импульсом тяги. Результаты

этих измерений используются в НИИ ПМЭ и ОКБ «Факел» для оценки воздействия струй указанных двигателей на элементы конструкции транспортного модуля, разрабатываемого в рамках ОКР «Двина ТМ» в НПО имени СЛ. Лавочкина, и на элементы конструкции КА «ЛучМ», разрабатываемого в ОАО «Информационные системы связи».

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

В частности показано, что с увеличением разрядного напряжения:

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

-в ядре потока возрастают значения концентрации и плотности тока

ионов;

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

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

Архипов A.C., Бишаев A.M. Численное моделирование в трехмерной постановке струи плазмы, выходящей в окружающее пространство из стационарного плазменного двигателя // Журнал вычислительной математики и математической физики, 2007, Т.47, №3, С.491-506.

2. Arkhipov A.S., Bishaev A.M., Kim V. Simulation of rarefied plasma dynamics in the plume of the stationary plasma thruster // Paper 5-SP-56 presented at the 2nd European Conference for Aerospace Sciences, Brussels Belgium, July 1-7,2007.

3. Arkhipov A.S., Kim V, Sidorenko E.K Analysis of Energy Balance in the Discharge of SPT Using Results of Its Integral Parameters and Plume Characteristics Measurements // Paper IEPC-2009-097 presented at the 31st International Electric Propulsion Conference, Ann Arbor, Michigan, USA, September 20-24,2009.

А.Архипов A.C., Бишаев A.M., Ким В. Влияние различных факторов на распределение параметров в струе стационарного плазменного двигателя // Математическое моделирование, 2010, Т.22, №2, С.124-138. 5. Архипов А. С., Ким В.П., Сидоренко Е.К. Анализ баланса энергии в разряде СПД по результатам измерений его интегральных параметров и характеристик струи // Вестник Московского авиационного института, 2010, Т. 17, №5, С.121-129.

Подписано в печать 01.11.201 Ог. Печать цифровая. Усл.п.л.1,5 Тираж 100 экз. Заказ № 777 Отпечатано в ООО «Реглет» 125315 г. Москва, Ленинградский проспект, д.74 к.1 Тел: 790-47-77; 661-60-89

Введение.

ГЛАВА 1. Анализ состояния разработки СПД и основные задачи диссертационной работы.

1.1 Состояние разработки СПД.

1.2 Основные направления дальнейшего развития СПД.

1.3 Результаты исследования струй СПД.

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

1.5 Цели и основные задачи диссертационной работы.

ГЛАВА 2. Экспериментальное исследование характеристик струй СПД, работающих при повышенных разрядных напряжениях.

2.1 Общая характеристика струй СПД.

2.2 Краткое описание стенда, методики и средств измерения параметров потока ускоренных ионов в струе СПД.

2.2.1 Вакуумный стенд У-2В-1 и его системы.

2.2.2 Электростатический зонд-энергоанализатор.

2.2.2.1 Принцип работы датчика.

2.2.2.2 Особенности функционирования датчика.

2.2.2.3 Конструкция датчика.

2.2.2.4 Недостатки датчика.

2.2.3 Основные измеряемые характеристики и определяемые параметры

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

2.3.1 Краткое описание моделей СПД, для которых проведены исследования характеристик струи.

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

2.3.3 Анализ закономерностей изменения характеристик струи двигателя при изменении режима его работы.

2.4 Выводы по главе.

ГЛАВА 3. Анализ зависимости интегральных параметров струи от разрядного напряжения и их взаимосвязи с тяговой эффективностью двигателя.

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

3.2 Зависимости интегральных параметров струи от разрядного напряжения

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

3.4 Анализ факторов, снижающих тяговую эффективность двигателя.

3.4.1 Неполнота использования рабочего тела в ускорительном канале двигателя.

3.4.2 Неполнота использования разрядного напряжения на ускорение ионов

3.4.3 Расходимость струи двигателя.

3.5 Выводы по главе.

ГЛАВА 4. Расчетное и экспериментальное исследование распределений параметров плазмы в ближней зоне струи СПД.

4.1 Расчетное исследование влияния разрядного напряжения на распределения параметров плазмы в струе СПД.

4.2 Экспериментальное исследование квазирадиальных ионных потоков.

4.3 Расчетное исследование влияния катода и закрутки ионов на распределения параметров плазмы в струе СПД.

4.4 Выводы по главе.

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

• межорбитальная транспортировка;

• поддержание длительной работоспособности КА на геостационарной орбите (ГСО) со сроками активного существования (САС) 12-15 лет;

• межпланетные перелеты и полеты в дальний космос.

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

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

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

С учетом изложенного в качестве целей диссертационной работы были выбраны следующие:

• уточнение закономерностей формирования распределений плотности тока ускоренных ионов по углу отклонения направлений их вылета из двигателей типа СПД-100 и СПД-140 отечественной разработки и распределений по энергии ионов, движущихся по различным направлениям при работе таких двигателей при повышенных разрядных напряжениях;

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

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

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

1. Провести анализ и исследование характеристик струи моделей современных двигателей типа СПД-100 и СПД-140 при их работе с повышенными разрядными напряжениями.

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

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

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

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

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

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

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

• результаты экспериментального определения характеристик струй двигателей, работающих при повышенных разрядных напряжениях, используются при оценке механического, эрозионного и загрязняющего воздействий двигателей типа СПД-100Д и СПД-140Д в ОКР «Двина ТМ» и «Луч М»;

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

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

Основные положения, выносимые на защиту

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

2. Положение о решающем влиянии температуры электронов на расходимость струи СПД.

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

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

Апробация работы и научные публикации

Результаты диссертационной работы докладывались и обсуждались на научно-технических семинарах в Научно-исследовательском институте прикладной механики и электродинамики (НИИ ПМЭ), 5-ой, 6-ой, 7-ой и 8-ой Международных конференциях по неравновесным процессам в соплах и струях, 14-ой, 15-ой и 16-ой Международных конференциях по вычислительной механике и современным прикладным программным системам, Международной конференции «Фундаментальные проблемы высокоскоростных течений», 3-ей, 6-ой и 7-ой Международных конференциях «Авиация и космонавтика», 2-ой Европейской конференции по аэрокосмическим наукам, 31-ой Международной конференции по электроракетным двигателям. Основное содержание и результаты диссертационной работы изложены в следующих публикациях:

1. Архипов А. С., Бишаее A.M. Численное моделирование в трехмерной постановке струи плазмы, выходящей в окружающее пространство из стационарного плазменного двигателя // Журнал вычислительной математики и математической физики, 2007, Т.47, №3, С.491-506.

2. Arkhipov A.S., Bishaev A.M., Kim V. Simulation of rarefied plasma dynamics in the plume of the stationary plasma thruster // Paper 5-SP-56 presented at the 2nd European Conference for Aerospace Sciences, Brussels, Belgium, July 1-7, 2007.

3. Arkhipov A.S., Kim V, Sidorenko E.K. Analysis of Energy Balance in the Discharge of SPT Using Results of Its Integral Parameters and Plume Characteristics Measurements // Paper IEPC-2009-097 presented at the 31st International Electric Propulsion Conference, Ann Arbor, Michigan, USA, September 20-24, 2009.

4. Архипов А. С., Бишаев A.M., Ким В. Влияние различных факторов на распределение параметров в струе стационарного плазменного двигателя // Математическое моделирование, 2010, Т.22, №2, С.124-138.

5. Архипов А.С., Ким В.П., Сидоренко Е.К. Анализ баланса энергии в разряде СПД по результатам измерений его интегральных параметров и характеристик струи // Вестник Московского авиационного института, 2010, Т.17, №5, С.121-129.

Вклад автора

Автор принимал непосредственное участие в экспериментах, в ходе которых определялись угловые распределения плотности ионного тока и энергетические спектры для моделей двигателей СПД-ЮО-О и СПД-140М, занимался обработкой и физическим анализом полученных данных. Автором была разработана методика расчета интегральных параметров струи СПД с учетом доли двухзарядных ионов и с ее использованием был проведен анализ зависимости коэффициентов, отражающих влияние на потери тяги и тягового к.п.д. (неполноты использования рабочего тела в ускорительном канале и неполноты использования приложенного разрядного напряжения на ускорение ионов, расходимости струи и разброса ионов по скоростям) от разрядного напряжения. Автором была разработана методика и экспериментально определены параметры квазирадиальных потоков ионов вблизи среза ускорительного канала опытного образца двигателя СПД-100Д, работающего на режиме с высоким удельным импульсом тяги, а совместно с в.н.с. A.M. Бишаевым разработана методика расчетной оценки указанных потоков и проведены расчеты, подтверждающие механизм их формирования.

Структура и объем диссертации

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

4.4 Выводы по главе

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

1. С увеличением разрядного напряжения: сужаются распределения концентрации и потоков ионов, что объясняет полученный экспериментально факт сужения струи при увеличении разрядного напряжения; в ядре потока возрастают значения концентрации и плотности потока, что объясняет отмеченное в Главе 2 увеличение плотности тока в приосевой части струи с увеличением напряжения.

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

3. Показана необходимость учета потока газа из катода и закрутки ионов в магнитном поле при будущих работах по моделированию процессов в ближней зоне струи.

ЗАКЛЮЧЕНИЕ

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

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

2. Разработана методика расчета коэффициентов, отражающих влияние неполноты использования рабочего тела в ускорительном канале и неполноты использования разрядного напряжения на ускорение ионов, наличия двухзарядных ионов, расходимости струи, разброса ионов по скоростям и общих «непроизводительных» затрат энергии в разряде на потери тяги и снижение тяговой эффективности. Получены новые данные о зависимости названных коэффициентов и доли двухзарядных ионов от разрядного напряжения. В результате показано, что: значения перечисленных коэффициентов, рассчитанных по результатам измерений характеристик струи дают значения тягового к.п.д., согласующиеся со значениями, полученными в результате прямых измерений интегральных параметров; коэффициент использования рабочего вещества в ускорительном канале СПД и среднее значение доли двухзарядных ионов монотонно возрастают с увеличением разрядного напряжения; полуугол расходимости струи несколько уменьшается с увеличением разрядного напряжения до (500-700) В, а затем слабо возрастает, при этом характер его зависимости от разрядного напряжения может быть объяснен зависимостью радиальной составляющей скорости ионов от температуры электронов и характером ее изменения в ускорительном канале при увеличении разрядного напряжения; средняя по струе энергия ионов несколько возрастает с увеличением разрядного напряжения, достигая значений порядка 80 % от энергии, соответствующей этому напряжению для однозарядных ионов; потери тяговой эффективности из-за расходимости струи составляют 1015 %, а разброс ионов по скоростям дает наименьшие (~5 %) потери в тяговой эффективности двигателя.

3. Совместно с другими сотрудниками НИИ ПМЭ проведены измерения характеристик струй опытных образцов двигателей типа СПД-100Д и СПД

140Д на режимах работы с высоким удельным импульсом тяги. Результаты этих измерений используются в НИИ ПМЭ и ОКБ «Факел» для оценки воздействия струй указанных двигателей на элементы конструкции транспортного модуля, разрабатываемого в рамках ОКР «Двина ТМ» в НПО имени С.А. Лавочкина, и на элементы конструкции КА «Луч М», разрабатываемого в ОАО «Информационные системы связи».

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

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

1. Морозов A.M. Исследование стационарного электромагнитного ускорения плазмы // Диссертация на соискание ученой степени доктора физико-математических наук, Москва, Институт атомной энергии имени И.В. Курчатова, 1965.

2. Козубскш КН., Мурашко В.М., Рылов Ю.П. и др. СПД работают в космосе // Физика плазмы, 2003, Т.29, №3, С.277-292.

3. Kim V., Kozubsky K.N., Murashko V.M., Semenkin A. V. History of the Hall Thrusters Development in USSR // Paper IEPC-2007-142 presented at the 30th International Electric Propulsion Conference, Florence, Italy, September 17-20, 2007.

4. Arhipov B.A., Vinogradov V.N., Kozubsky K.N. et al. Development and Application of Electric Thrusters at EDB «Fakel» // Paper IEPC-1997-004iLpresented at the 25 International Electric Propulsion Conference, Cleveland, Ohio, USA, August 24-28, 1997.

5. Arkhipov В., Kim V., Koryakin A. et al. Small SPT Development and Tests // Proceedings of the 3rd International Conference on Spacecraft Propulsion, Cannes, France, 2000, pp.399-401.

6. Manzella D., Jacobson R., Hofer R. Laboratory Model 50 kW Hall Thruster // Paper AIAA-2002-3676 presented at the 38th Joint Propulsion Conference, Indianapolis, Indiana, USA, July 7-10, 2002.

7. Горшков О.А., Муравлев В.А., Шагайда А.А. Холловские и ионные плазменные двигатели для космических аппаратов // Москва, «Машиностроение», 2008, 280с.

8. TacSat-2 // Материалы электронной энциклопедии космонавтики (Encyclopedia Astronáutica), http://www.astronautix.com/craft/tacsat2.htm.

9. Myers R.M. Overview of Major U.S. Industrial Electric Propulsion Programs // Paper AIAA-2004-3331 presented at the 40,h Joint Propulsion Conference, Fort Lauderdale, Florida, USA, July 11-14, 2004.

10. McLean C.H., McVey J.B., Lazarovici C.E. et al. Innovative Low Power Hall Effect Thruster // Paper AIAA-2000-3247 presented at the 36th Joint Propulsion Conference, Huntsville, Alabama, USA, July 16-19, 2000.

11. Kim V., Kozlov V., Lazurenko A. et al. Development and Characterization of Small SPT // Paper AIAA-1998-3335 presented at the 34th Joint Propulsion Conference, Cleveland, Ohio, USA, July 13-15, 1998.

12. Belikov M.B., Gorshkov O.A., Jakupov A.B., Khartov S.A. Experimental Research of SPT Low-Power Perspective Model // Paper AIAA-1998-3786 presented at the 34th Joint Propulsion Conference, Cleveland, Ohio, USA, July 13-15, 1998.

13. Biagioni L., Cesari U., Saverdi M. et al. Development Status of the HT-100 Miniaturized Hall Effect Thruster System // Paper AIAA-2005-3682 presented at the 41st Joint Propulsion Conference, Tucson, Arizona, USA, July 10-13,2005.

14. Loyan A. V., Maksymenko T.A. Performance Investigation of SPT-20M Low Power Hall Effect Thruster // Paper IEPC-2007-100 presented at the 30lh International Electric Propulsion Conference, Florence, Italy, September 1720, 2007.

15. Fujioka Т., Kitano Т., Tahara H. et al. Optimization on Magnetic Field and Acceleration Channel for Low Power Hall Thrusters // Paper IEPC-2003-015 presented at the 28th International Electric Propulsion Conference, Toulouse, France, March 17-21, 2003.

16. Choe W., Lee J., Seo M. et al. Development of a 10 mN Class Hall Thruster for STSAT-3 // Paper IAC-09-C4.4.1 presented at the 60th International Astronautical Congress, Daejeon, Republic of Korea, 12-16 October, 2009.

17. Albarede L., Bouchoule A., Lazurenko A. et al. Characterization ofPPS-1350 Type Thrusters Under Increased Discharge Voltages and Comparison with Hybrid Codes Simulation Results // Paper IEPC-2005-136 presented at theL

18. International Electric Propulsion Conference, Princeton, New Jersey, USA, October 31 November 4, 2005.

19. Kim V. et al. Modern Trends of Electric Propulsion Activity in Russia // Paper IEPC-1999-004 presented at the 26th International Electric Propulsion Conference, Kita-Kyushu, Japan, October 17-21, 1999.

20. Bouchoule A., Kim V., Kozlov V. et al. Investigation of the SPT Operation Particularities Under High Discharge Voltages // Paper IEPC-2003-211 presented at the 28th International Electric Propulsion Conference, Toulouse, France, March 17-21, 2003.

21. Kim- V., Kozlov V., Skrylnikov A. et al. Investigation of the Local Plasma Parameter Distributions in the SPT Accelerating Channel Under Increased

22. Discharge Voltages // Paper IEPC-2005-004 presented at the 29th International Electric Propulsion Conference, Princeton, New Jersey, USA, October 31 November 4, 2005.

23. Kim V., Kozlov V., Shylnikov A. et al. Investigation of the SPT Operation Under High Discharge Voltages // Paper IEPC-2005-004 presented at the 29 International Electric Propulsion Conference, Princeton, New Jersey, USA, October 31 November 4, 2005.

24. Bouchoule A., Kim V., Kozlov V., Lazurenko A., Skrylnikov A., Vial V. DualMode Operation of Stationary Plasma Thrusters // Journal of Propulsion and Power, January-February 2006, Vol.22, No.l, pp.3 8-47.

25. Kim V., Kozlov V., Popov G., Skrylnikov A., Sidorenko E., Umnitsin L. High Voltage SPT Studies // Paper SP2008-127 presented at the 5th International Spacecraft Conference, Heraclion, Crete, Greece, May 5-9, 2008.

26. Kim V., Kozlov V., Popov G., Skrylnikov A. Plasma Parameter Distribution in SPT-70 Plume // Paper IEPC-2003-107 presented at the 28th International Electric Propulsion Conference, Toulouse, France, March 17-21, 2003

27. Manzella D.H. Stationary Plasma Thruster Ion Velocity Distribution // Paper AIAA-1994-3141 presented at the 30th Joint Propulsion Conference, Indianapolis, Indiana, USA, June 27-29, 1994.

28. Gallimore A.D., King L.B., Marrese C.M. Transport-Property Measurements in the Plume of an SPT-100 Hall Thruster // AIAA Journal of Propulsion and Power, May-Junel998, Vol.14, No.3, pp.327-335.

29. Gallimore A.D., King L.B. Ion Energy Diagnostics in the Plasma Exhaust Plume of a Hall Thruster // AIAA Journal of Propulsion and Power, September-October 2000, Vol.16, No.5, pp.916-922.

30. Gallimore A.D., King L.B. Mass Spectral Measurements in the Plume of an SPT-100 Hall Thruster // AIAA Journal of Propulsion and Power, November-December 2000, Vol.16, No.6, pp. 1086-1092.

31. Gallimore A.D., Gilchrist B.E., Ohler S. Microwave Plume Measurements of a Closed Drift Hall Thruster // AIAA Journal of Propulsion and Power, November-December 1998, Vol.14, No.6, pp. 1016-1021.

32. Boyd I. Hall Thruster far Field Plume Modeling and Comparison to Express Flight Data // Paper AIAA-2002-0487 presented at the 40th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, USA, January 14-17, 2002.

33. Arkhipov В., Day M., Gnizdor R., Kim V. et al. SPT-100 Module Lifetime Test Results // Paper AIAA-1994-2854 presented at the 30th Joint Propulsion Conference, Indianapolis, Indiana, USA, June 27-29, 1994.

34. Day M., Kim V., Kozlov V., Popov G. et al. Investigation of the Nonoperating Cathode Erosion Reasons // Paper AIAA-1996-2710 presented at the 32nd Joint Propulsion Conference, Lake Buena Vista, Florida, USA, July 1-3, 1996.

35. Борисов Б.С., Корсун А.Г. II Тезисы докладов 5-ой Всесоюзной конференции по плазменным ускорителям и ионным инжекторам, Москва, «Наука», 1982.

36. Gabdullin F.F., Korsun A.G., Tverdokhlebova E.M. Mathematical Model of Hypersonic Plasma Flows Expanding in Vacuum II Computer Physics Communications, 2004, Vol.164, No. 1-3, pp.434-441.

37. Жевандров П.Н., Морозов A.M., Якунин С.А. Динамика плазмы, образующейся при ионизации разреженного газа // Физика плазмы, 1984, Т.10, Вып.2, С.353-360.

38. Khartov S., Metois P., Nadiradze A., Perrin V. Simulation Tools for the Plasma Propulsion and Satellite Environment // Paper presented at the 52nd International Astronautical Congress, Toulouse, France, October 1-5, 2001.

39. Bird G.A. Molecular Gas Dynamics and the Direct Simulation of Gas Flows // Oxford University Press, Oxford, 1994.

40. Haas J.M., Hastings D.E., Gallimore A.D., Marrese C.E., Oh D.Y. Modeling of Stationary Plasma Thruster-100 Thruster Plumes and Implication for Satellite Design // AIAA Journal of Propulsion and Power, March-April 1999, Vol.15, No.2, pp.345-357.

41. Boyd I., Carrigues L., Koo J., Keidar M. Progress in Development of a Combined Device/Plume Model for Hall Thrusters // Paper AIAA-2000-3520 presented at the 36 Joint Propulsion Conference, Huntsville, Alabama, USA, July 16-19, 2000.

42. Cai C., Boyd I., Sunz Q. Three-Dimensional Particle Simulation of Plume Flows From Hall Thrusters // Paper IEPC-2005-048 presented at the 29th International Electric Propulsion Conference, Princeton, New Jersey, USA, October 31 November 4, 2005.

43. Бишаев A.M., Калашников В.К., Ким В. Численное исследование струи разреженной плазмы стационарного ускорителя с замкнутым дрейфом электронов (УЗДП) // Физика плазмы, 1992, Т. 18, Вып.6, С.698.

44. Бишаев A.M., Калашников В.К, Ким В., Шавыкина А.В. Численное моделирование плазменной струи стационарного плазменного двигателя, распространяющейся в среде низкого давления // Физика плазмы, 1998, Т.24, №11, С.989-995.

45. Архипов А.С., Бишаев A.M. Численное моделирование в трехмерной постановке струи плазмы, выходящей в окружающее пространство из стационарного плазменного двигателя // Журнал вычислительной математики и математической физики, 2007, Т.47, №3, С.491-506;

46. Коган М.Н. Динамика разреженного газа // Москва, Главная редакция физико-математической литературы, 1967.

47. А. Бишаев A.M. Применение методов кинетической теории для решения задач разреженных газов и плазмы // Диссертация на соискание ученойстепени доктора физико-математических наук, Москва, ФГНУ «НИИ ПМЭ», 2005.

48. Волков Б.И., Морозов А.И., Свешников А.Г., Якунин С.А. Численное моделирование ионов в системе с замкнутым дрейфом // Физика плазмы, 1981, Т.7, Вып.2, С.245-253.

49. Лившиц Е.М., Питаевский Л.П. Физическая кинетика // Москва, «Наука», 1979.

50. Vahrenkamp R.P. Measurement of Double Charged Ions in the Beam of a 30-cm Mercury Bombardment Thruster // Paper AIAA-1973-1057 presented at the AIAA 10th Electric Propulsion Conference, Lake Tahoe, Nevada, USA, October 31 November 2, 1973.

51. Gallimore A., King L. Propellant Ionization and Mass Spectral

52. Bareilles J. Modélisation 2D Hybride d'un Propulseur a Effet Hall pour Satellites // Ph.D. Thesis, Université Paul Sabatier, Toulouse, 2003, 139p.

53. Леонард С., Хаддлстоун Р. Диагностика плазмы // Москва, «Мир», 1967.

54. Архипов А.С., Бишаев А.М., Ким В. Влияние различных факторов на распределение параметров в струе стационарного плазменного двигателя // Математическое моделирование, 2010, Т.22, №2, С. 124138;