Моделирование ионной эрозии стенок канала разрядной камеры стационарного плазменного двигателя тема диссертации и автореферата по ВАК РФ 05.07.05, кандидат технических наук Абгарян, Вартан Карленович

• Абгарян, Вартан Карленович
• кандидат технических наук
• 2009

Оглавление диссертации кандидат технических наук Абгарян, Вартан Карленович

ВВЕДЕНИЕ.

ГЛАВА 1. СОСТОЯНИЕ РАЗРАБОТКИ СПД И ОСНОВНЫЕ ЗАДАЧИ

РАБОТЫ.

1.1. Введение.

1.2. Схема стационарного плазменного двигателя.

1.3. Эрозия стенок канала разрядной камеры в СПД.

1.4. Механизмы распыления поверхностей керамических материалов.

1.4.1. Электронное распыление.

1.4.2. Ионное распыление поверхностей твердых тел.

1.5. Зависимости коэффициента распыления.

1.5.1. Энергетическая зависимость.

1.5.2. Зависимости коэффициента распыления от угла падения ионов.

1.5.3. Пространственные распределения распыленного вещества.

1.5.4.Зависимость коэффициентов распыления от температуры мишени.

1.6. Моделирование эрозии стенок канала разрядной камеры СПД.

1.6.1. Поэтапная методика.

1.6.2. Модели расчета эрозии стенок РК.

1.7. Изменение микрорельефа облучаемых ионами поверхностей.

1.7.1. Механизмы возникновения рельефа на гладкой поверхности.

1.7.2. Регулярные структуры на поверхности.

1.7.3. Модели возникновения волнообразного рельефа.

1.7.4. Перенапыление распыленного материала облучаемой ионами протяженной поверхности. Граничные условия.

1.8. Потоки распыляемого материала стенок канала РК СПД.

1.8.1. Потоки распыляемого материала стенок вне канала.

1.8.2. Потоки распыляемого материала стенок внутри канала.

1.9. Постановка задачи.

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

КАМЕРЫ СПД.

2.1. Введение.

2.2. Теоретическое описание модели.

2.2.1. Допущения. Обозначения.

2.2.2. Уравнение эрозии поверхностей стенок канала.

2.2.3. Численное интегрирование уравнения эрозии.

2.2. 4. Пределы интегрирования по азимутальному углу.

2.2. 5. Затенение потоков перенапыления на стенки канала РК СПД.

2.2. 6. Распыляющая способность.

2.2.7. Углы падения ионов на стенки канала РК. 2.2.8. Угловая зависимость относительного коэффициента распыления.

2.2.9. Математическая постановка задачи.

2.3. Оптимизационный метод и параметры оптимизации.

2.3.1. Параметры оптимизации.

2.3.2. Процедура оптимизации параметров модели.

2.4. Программа и различные варианты расчета профилей эрозии.

2.5. Результаты расчетов.

2.5.1. Эрозия стенок канала РК СПД. i' 2.5.2. Поток распыляемого материала вне канала.

2.6. Выводы.Ю<

ГЛАВА. 3. МОДЕЛЬ РАСЧЁТА ИЗМЕНЕНИЯ МИКРОГЕОМЕТРИИ

ОБЛУЧАЕМОЙ ИОНАМИ ШЕРОХОВАТОЙ ПОВЕРХНОСТИ. 10:

3.1. Уравнение эрозии облучаемой поверхности.10!

3.1.1. Вывод уравнения эрозии.10!

3.1.2. Потоки перенапыляемых частиц.10:

3.2. Граничные условия. Переход к криволинейным координатам.10'

3.3. Особенности модели в криволинейной геометрии.

3.3.1. Геометрия модели и траектории частиц.

3.3.2. Расстояния между элементами на облучаемой поверхности.11!

3.3.3. ГТоправка на кривизну поверхности.11:

3.4. Построение шероховатой поверхности статистическим методом.1 Ь

3.5. Потоки перенапыления при возможном затенении промежуточными элементами поверхности.11!

3.5.1. Методика определения затенения потока перенапыления промежуточными элементами поверхности.Ill

3.5.2. Разбиение поверхности на элементы и особенности расчета на ЭВМ 12(

3.6. Результаты расчетов. Выводы.12!

Введение диссертации (часть автореферата) на тему «Моделирование ионной эрозии стенок канала разрядной камеры стационарного плазменного двигателя»

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

Стационарные плазменные двигатели (СПД) успешно используются в системах коррекции орбит космических аппаратов с длительным сроком активного существования. Ресурс современных СПД находится на уровне 10000 часов. При сертификации двигателя для подтверждения ресурса двигателя необходимо проведение полномасштабных испытаний с названной длительностью. Стоимость таких испытаний составляет величину порядка 10б $ США и, таким образом, проблема обеспечения ресурса является в процессе разработки современных СПД одной из наиболее сложных и дорогостоящих.

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

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

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

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

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

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

Цели и основные задачи

Основная цель работы состояла в

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

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

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

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

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

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

Научная новизна работы

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

1. Построена расчетная модель, позволяющая на основе оптимизационного подхода с использованием экспериментальных профилей проводить расчёты' эрозии стенок РК с высокой точностью для двигателей параметрического ряда современных СПД. Что продемонстрировано на примере двигателя СПД-100.

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

3. Произведены расчёты эрозии стенок разрядной камеры с учетом напыления распыленного материала в канале РК СПД и показано, что толщина"напыленной в прианодной области пленки двигателя типа СПД-100 за время работы, равной ресурсу, может достигать 0.5 мм.

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

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

Научная и практическая ценность

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

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

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

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

Применение нового метода замыкания границ расчетной области (РО) в разработанной методике уменьшает компьютерные затраты почти на порядок для 2х-мерной и примерно на полтора порядка для Зх-мерной РО по сравнению с затратами при применении периодических граничных условий, обычно используемых в такого рода задачах. Такой метод позволяет повысить физическую адекватность расчетной модели, и при равной погрешности обеспечивает значительный выигрыш в счете.

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

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

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

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

Апробация работы

Материалы диссертационной работы докладывались и обсуждались на Международных конференциях по ионным двигателям в Индианаполисе в 1994 и в 2003гг., XXVIII-й Международной конференции по физике взаимодействия заряженных частиц с кристаллами, Москва, 1998 г, XIV Международной конференции по взаимодействию ионов с поверхностью (Звенигород, 1999), X Международном совещании «Радиационная физика твердого тела» (Севастополь, 2000г.), XV Международной конференции по взаимодействию ионов с поверхностью (Звенигород, 2001), XVI Международной конференции по взаимодействию ионов с поверхностью (Звенигород, 2003), XVII Международной конференции по взаимодействию ионов с поверхностью (Звенигород, 2005), Международной конференции по ионным двигателям в г. Кальяри, Сардиния, 2-4 июня, 2004, 4-й международной конференции «Авиация и космонавтика», (Москва, 10-13 октября 2005г.), Международной конференции STAR-2006 «Компьютерные технологии решения прикладных задач тепломассопереноса и прочности», (г.Саров,4-6 апреля 2006г.), XVIII Международной конференции по взаимодействию ионов с поверхностью (Звенигород, 2007)

Публикации

По материалам диссертации опубликовано 18 работ, список основных из которых приведён в конце автореферата. и

Вклад автбра

Личный вклад автора в работы, написанные в соавторстве и вошедшие в диссертацию, в большинстве работ является определяющим.

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

Диссертация состоит из введения, трех глав и заключения. Она содержит 137 страниц текста, в том числе 35 рисунок и список литературы из 109 наименований.

Заключение диссертации по теме «Тепловые, электроракетные двигатели и энергоустановки летательных аппаратов», Абгарян, Вартан Карленович

ЗАКЛЮЧЕНИЕ

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

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

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

3. Показано, что за время работы СПД в прианодной области канала РК напыляется пленка из распыляемого материала толщиной до 0.5мм.

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

5. Для адекватного численного моделирования распыления протяженной (в плоскости) поверхности предложены новые граничные условия. Разработана процедура замыкания границ расчетной области поверхности, благодаря чему, граничные условия исключены из рассмотрения. Такой прием позволяет на порядок сократить объем счета. Его можно применять в задачах, где размеры расчетной области много больше характерных размеров задачи (здесь размеры шероховатости), и присутствует массо- или обмен излучением между различными областями рассматриваемых поверхностей или объемов.

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

Автор считает своим приятным долгом выразить благодарность профессору В.Киму за содействие, оказанное на всех этапах выполнения работы от постановки задачи до программной реализации модели расчета эрозии в СПД; профессору И.И. Шкарбану за многолетнее сотрудничество и помощь при выполнении представленной работы. Также автор благодарит А.А. Семенова, Т.В. Забровскую за техническую помощь, а также Е.В.Лукьянову и А.Ю. Фельдмана за всестороннее содействие, постоянно оказываемого при выполнении работы.

Список литературы диссертационного исследования кандидат технических наук Абгарян, Вартан Карленович, 2009 год

1. Козубский КН., Мурашко В.М., Рылое Ю.П., Трифонов Ю.В., Ходненко

2. B.П., Ким В., Попов Г.А., Обухов В.А. СПД работают в космосе // Физика плазмы, том 29, №3, 2003, С. 277 292.

3. Плазменные ускорители/ Под ред. Л.А.Арцимовича. М.: Машиностроение, С.5, 1973г.

4. Л.А.Арцилювич, JI.A., Морозов А.И., Снарский Р.К.// Космические исследования, 1974г.,т.12, №3, С.451

5. А.Морозов А.И. Физические основы космических электрореактивных двигателей// в кн. «Элементы динамики потоков в ЭРД, С.252, М., Атомиздат, 1978.

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

7. А.И.Морозов, Ю.В.Есипчук, А.М.Капулкин. Влияние конфигурации магнитного поля на режим работы ускорителя с замкнутым дрейфом электронов (УЗДП) //ЖТФ, 1972, т. XLII, вып.З, С.612-619.

8. С.Д.Гришин, Л.В.Лесков, Н.П.Козлов. Электрические ракетные двигатели// М.: Машиностроение, 1975, С.167-181.

9. М.В. Belikov, A.I. Vasin, О.A. Gorshkov, V.A. Muravlev, R.N. Rizcikhanov. LOW-POWER HALL-EFFECT THRUSTER for small spacecrafts/ZProceedings of International Symposium on Space Propulsion, P. 391, Shanghai, China, 2004.

10. В. Ким, А.И. Скрылъников, E.K Сидоренко. Оценка мощности, выделяющейся на стенках разрядной камеры СПД// Труды конференции

11. Авиационно-космическая техника и технология, №10/36, ХАИ, 2006, с.112-116

12. S.A.Khartov, D.M.Napolov, A.A.Perfiliev, J.V.Zikeeva. Experimental Investigation of the Possibility Argon and Oxygen Using as a Propellant for the SPT// Proceedings of 3rd International Conference "Space Propulsion", Cannes, France, 2000.

13. Г.Н. Тилинин. Экспериментальное исследование высокочастотных плазменных колебаний в УЗДП// Журнал технической физики, 1976, т.47,В.8, С.1684-1691.

14. Кирдяшов К.П., Бугрова А.К, Десятсков А.В., Морозов А.И. СВЧ -колебания в ускорительном канале СПД-Атон// Письма в ЖЭТФ. 2005, Т.31,В14, С.7-15.

15. Бигиаев A.M., Гаврюшин В.М. Герасимов В.Ф. Исследование влияния колебаний в разрядной цепи и проводимости стенок канала на процессы в УЗДП// в кн.:1У Всесоюзн. конф. по плазменным ускорителям и ионнымvинжекторам. Тезисы докладов ВНТИЦД978, с. 68-69.

16. В.И Козлов. Динамика электронов в разряде стационарного плазменного двигателя//Вестник МАИ, 2007, т. 14, №3, С.75-94.

17. Морозов A.M., Есипчук Ю.В., Капулкин, Невровский В.А., Смирнов В.А. Азимутально-несимметричные колебания и аномальная проводимость в ускорителе с замкнутым дрейфом электронов (УЗДП)// Журнал технической физики, 1973, T.XL1II, В.5, С.972-982.

18. Архипов Б А. Исследование катодов нового поколения для стационарных плазменных двигателей (СПД)// Диссертация на соискание ученой степени доктора технических наук, Калининград, 1998г, С.302.

19. Бишаев А.М, Ким В. Исследование локальных параметров плазмы в ускорителе с замкнутым дрейфом электронов и протяженной зоной ускорения// ЖТФ, 1978,48,№9, С. 1853-1857.

20. Л.М.Бйшаев, В.Ким. Исследование влияния конфигурации стенок канала на распределение параметров плазмы в УЗДП// в сб.: Источники и ускорители плазмы, Харьков: ХАИ, 1981, вып.5, стр. 3-8.

21. А.И.Бугрова, В.Кгш. Современное состояние физических исследований в ускорителях с замкнутым дрейфом электронов и протяженной зоной ускорения// в сб. Плазменные ускорители и ионные инжекторы, М.; Наука, 1984, С.107-129.

22. Kim В. Анализ закономерностей износа изолятора, ограничивающего ускорительный канал УЗДП// в кн. Источники и ускорители плазмы. Харьков: ХАИ, 1982г., вып.5, С.7

23. Кгш В. Физико-технические основы проектирования стационарных плазменных двигателей (СПД)// Диссертация на соискание ученой степени доктора технических наук, Москва, 1984, С.411.

24. Баранов В.К, Васин А.И., Петросов В.А. Закономерности износа изолятора степени канала СПД// в науч. техн. сборн. «Ракетно-космическая техника» вып. 3(131), С. 54-70, НИИТП, Москва 1991 г.

25. В.Ким, В.И.Козлов, В.М.Мурашко. Исследование особенностей работы стационарных плазменных двигателей при повышенных разрядных напряжениях // Космонавтика и ракетостроение, 33(52). С.21-27, 2008г.

26. А.А.Семенов, С.А.Ситников, А.Ю. Фельдман, И.И. Шкарбан. Влияние концентрации нитрида бора на коэффициент распыления композиционныхкерамических материалов// Материалы 16-ой Международной конференции

27. Взаимодействие ионов с поверхностью» (ВИП-2003), 25-29 августа 2003 г., г. Звенигород, М.: 2003. т. 1, С. 180.

28. С.С.Еловиков, Е.Ю.Зыкова, С.А. Постников, Ю.А.Рыжов, И.И.Шкарбан, В.Е.Юрасова. Радиационная стойкость нитрида бора и керамик на его основе к низкоэнергетическому электронному облучению // Известия РАН, Сер. Физическая, Т.71 (2007), №5, С.761-764.

29. В.К.Абгарян, С.Ю. Михеев, М.В. Прокофьев, Ю.А. Рыжов, И.И.Шкарбан. Массовые спектры частиц, эмитируемых с поверхностей керамик, облучаемыхплазменными потоками // Известия РАН.Серия физическая. Т. 70. № 6, С.879-882, 2006г.

30. Н.В. Плегиивцев. Катодное распыление// М.: Атомиздат, 1968,

31. Плегиивцев Н.В., Бажин A.M. Физика воздействия ионных пучков на материалы //М. Вузовская книга 1998.

32. В.Е. Юрасова. Взаимодействие плазмы с поверхностями твердых тел. Избранные труды//М.: 1999г.

33. Е.С.Машкова. Современные тенденции в исследовании распыления твердых тел// в кн. "Фундаментальные и прикладные аспекты распыления твердых тел",- под ред. Е.С. Машковой, М.: Мир, 1989

34. Хофер В. Распределения распыленных частиц по углам, энергиям и массам// в кн "Распыление под действием бомбардировки частицами", под. Ред. Р. Бер.иша и К. Виттмаака., вып. Ill, М. Мир, 1998, с. 26

35. П.Зигмунд.в кн. «Распыление твердых тел ионной бомбардировкой», т.1. Под ред. Бериша Р.М.-М., Мир, 1984

36. П.Зигмунд, в кн. «Распыление под действием бомбардировки частицами»,-вып. III под ред. Р. Бериша, К. Виттмака М. Мир, 1998 г.

37. М.Т.Робинсон, в кн. «Распыление твердых тел ионной бомбардировкой». т.1. Под ред. Бериша Р.М.-М., Мир, 1984

38. В. Экштайн. Компьютерное моделирование взаимодействия частиц с поверхностью твердого тела// М.Мир, 1995г.

39. С.С. Еловиков, Е.Ю. Зыкова, А.С. Мосунов, А.А. Семёнов, И.И. Шкарбан, В.Е. Юрасова. Зависимость распыления нитрида бора от энергии, массы и угла падения ионов. // Известия РАН. Серия Физическая, 2002, том 66, № 4, с. 558-561 v

40. Фирсов О.Б. О зависимости распыления мишени от угла падения бомбардирующих ионов// ДАН СССР, 1969, т. 189, С. 302.

41. B.Rubin, J.L.Topper, A.P.Yalin. Total and Differential Yields of Boron Nitride Measured by Quartz Crystal Microbalance and Weight Loss // Proceedings of International Electric Propulsion Conference, p.74, Florence, Italy, September 1720, 2007.

42. А.С.Мосунов, Ю.А.Рыэ1Сов, А.А.Семёнов, И.И.Шкарбан, Д.С.Коллигон, В.Е.Юрасова. Угловое распределение частиц, распыленных медленными ионами. //Известия РАН, Сер. Физическая, Т.68 (2004) №11, 1665-1672.

43. B.K. Абгарян, И.И. Шкарбан. Компьютерное моделирование ионного облучения монокристаллических поверхностей произвольной ориентации. // Известия АН. Сер. Физическая, 1998, т.62, №4, С.716-722.

44. Рыжов Ю.А., Шкарбан И.И. Тепло- и массообмен между потоками газов и твердой поверхностью// М. Изд. МАИ, 1975, с. 23.

45. Рыжов Ю.А. В книге: «Проблемы механики и теплообмена в космической технике М. Машиностроение, 1982, с. 99.

46. В.К. Абгарян, В.М. Арбатский, А.Б. Надирадзе, Е.Б. Паршина, В.В. Шапошников. II Труды XVII-й Международной конференции «Взаимодействие ионов с поверхностью», Звенигород, 2005.

47. Перроте А.И., Сторчак М.А. Ускоренные испытания изделий на надежность//М.: Машиностроение 1975 г.,45 с.

48. Пешее Л.Я., Степанова М.Д. Основные теории ускоренных испытаний на надежность // Минск: Наука и техника, 1972 г., 166 с.

49. Бугрова А.И., Версог^кий B.C., Ниекин В.Т. Способ локального форсирования скорости износа стенок канала ускорителей плазмы // в сб. Ракетно-космическая техника НИИТП, 1991 г. Вып. 3 (31), С. 89.

50. Латышев Л.А., Соколоверов А.П., Хартов С.А., Чуян Р.К Поэтапная методика испытаний ускорителей с замкнутым дрейфом электронов // в сб. Ракетно-космическая техника НИИТП. 1991 г. Вып. 3 (31), С. 71.

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

52. R.Hofer, /. Micalidis, I. Kats, D. Gobel. BPT-4000 Hall Truster Discharge Chamber Erosion Model Comparison with Qualification Life Time Test Data,

53. Proceedings of International Electric Propulsion Conference, p.267, Florence, Italy, September 17-20, 2007.

54. S. Chen, M. Matrines-Sanchez. Modeling of Hall Thruster Lifetime and Erosion Mechanisms // Proceedings of International Electric Propulsion Conference, p.250, Florence, Italy, September 17-20, 2007.

55. M. Gamero-Castafio, I.Kaiz. Estimation of Hall Thruster Erosion Using HPHall // Proceedings of 29th International Electric Propulsion Conference, p.303, Princeton University, USA, October 31-November 4, 2005.

56. J.T.Yim, I,D.Boyd, M. Keidar. Hall Thruster Erosion Prediction Using A Hydrodinamic Plazma Model and Sputtering Simulation// Proceedings of International Electric Propulsion Conference, p.34, Florence, Italy, September 1720,2007.

57. E. Sornmer, M.K.Allis, M.A.Cappelli. Wall Erosion in 2D Thruster Simulation // Proceedings of 29 International Electric Propulsion Conference, p. 189, Princeton University, USA, October 31-November 4, 2005.

58. Калмыков Б.М., Рыжов Ю.А. О деформации поверхностей, уносимых под действием ионной бомбардировки // ПМТФ, с. 132, 1967.

59. Гришин С.Д., Косогоров А.В., Марахтанов М.К, Суслов Н.И. Обработка поверхности низкоэнергетичными ионами плазменных ускорителей // Оптико-механическая промышленность, с. 30, № 11, 1979 г.

60. А.М.Борисов, Е.С.Машкова, А.С. Немое, С.А. Камнева, В.А. Курнаев, Н.Н.Трифонов. Влияние ионно-индуцированного рельефа на высокодозное распыление графита// Вопросы атомной науки и техники. Серия термоядерный синтез, вып.2,2004, с.65-72.г*

61. Wittmaack. Effect of surface roughening on secondary ion yield and erosion rates of silicon to oblique oxigen bombardment // S. Vac. Sci. Technol. 1990, v. A8, pp. 2246 2250.

62. В.П. Ким, А.А. Семёнов, А.Ю. Фельдман, И.И. Шкарбан. Распыление керамик ионами Хе и Кг при энергиях 100-400 эВ. // Сб. трудов 15 Международной конференции «Взаимодействие ионов с поверхностью» -Москва, 2001, т.1, с.134-136.

63. Васильев М.А. Структура и динамика поверхностных переходных металлов// 1985г., Киев, Наукова думка, С.289.

64. Емельянов В.И. Дефектно-деформационная самоорганизация поверхностных упорядоченных структур при ионном и лазерном облучении // Труды XVII Международной конференции «Взаимодействие ионов с поверхностью», Звенигород, 2005, т. 1, стр. 48-51.

65. Шелякин Л.Б., Шульце Р. Д.Г., Юрасова В.Е. Некоторые особенности формирования поверхностной структуры монокристаллов при ионной бомбардировке // Физика плазмы ,1975 ,т.1, с. 488-495.

66. Begrambekov L.B., Zakharov A.M., Telkovsky V.G. Peculiarities and mechanism of cone growth under ion bombardment // Nucl. Inst. And Meth. Phys. Res. В 115, 1996, p. 456-460.

67. Carter &. Physics and applications of ion beam erosion // J. Phys. D: Appl. Phys. 2001 v.34 p. R1-R22.

68. Stevie F.A., Kahova P.M., Simons P.S., Chi P. Secondary ion yield changes in Si and Ga As due topography changes during Cb+ or Cs+ ion bombardment // J. Vac. Sci. Techn. 1988, v. A7, p. 76-80.

69. Смирнов В.К., Лепшин П.А., Кривелевич С.А., Кибалов Д.С. Зависимость процесса рельефообразования при бомбардировке кремния ионами азота от энергии и угла бомбардировки // Неорганические материалы, 1988 г., т. 34, с. 1081-1084.

70. Бачурпн В.И. Распыление твердых тел ионами инертных и химически активных газов при фазовых превращениях// Диссертация на соискание ученой степени д.ф.м.н. 2007, Ярославль.

71. Chim К., Sanyal М.К., Bhattacharyya. Energy dependent wavelength of the ion induced nanoscale ripple // Phys. Rev., 2002, v. B66, p. 153404 -1-4.

72. Carter G., Vishnyakov V.M., Nobes M.Y. Ripple topography development on ion bombarded Si // Y. Appl. Phys., 1995, v. 78, p. 3559 3565.

73. Mayer T.M., Chason E., Howard A.Y. Roughening instability and ion induce viscous relaxation on Si02 surface // Y. Appl. Phys. 1994, v. 76, p. 1633 - 1643.

74. Chason E., Chan W.L. Kinetic mechanism in ion-induce ripple formation on Cu(001) surface // NIMB, 2006, v. 242, p. 232-236.

75. Rusponi S., Boragho C., Valbusa U. Ripple structure on metal surfaces, induced by ion sputtering // Phys. Low Dim. Struct. 1998, v. 11/12, p. 52-64.

76. Bradley R.M., Haper J. M.E. Theory of ripple topography induced by ion bombardment // J. Vac. Sci. Technol. 1988, v. A6, p.p. 2390-2395.

77. Micheky Т. Nanostructures by grazing incidence ion: ripple patterns; a thermal coarsening and subsurface channeling // Труды XVIII Международной конференции «Взаимодействие ионов с поверхностью», Звенигород, 2007, т. 1, с. 64-67.

78. Carter G. Viscoelastic relaxation and sputter-depth profiling of amorphous materials // SIA., 1997 v. 25, p. 36-40.

79. Rudy A.S., Smirnov V.K. A model of wave like structures formed by ion sputtering of amorphous materials 11 Physica Status Solidi (b), 199, v. 213, p. RI -R2.

80. Beslov-Y., Wilson L.H. Theory of redeposition of sputtered flux on surface asperities//NIMB. 182/183, 1981, p. 275-281.

81. В.Н.Опекунов, С.Ю.Михеев, Ю.А.Рыжов, И.И.Шкарбан, В.Е.Юрасова. Стадии эрозии поверхности при ионном облучении. // Сб. трудов 18-й Международной конференции «Взаимодействие ионов с поверхностью». Звенигород, 2007, Т.1, С.182-185.

82. V.N.Opelcunov, S.Yu.Mikheev, V.E.Yurasova, Yu.A.Ryzhov, A.A.Khaidarov, I.I.Shkarban. Stages of surface erosion under ion irradiation. // Vacuum, 82 (2008), P.352-358.a»

83. В.К.Абгарян. Граничные условия при численном моделировании ионного распыления поверхности шероховатых тел // Поверхность. Рентгеновские, синхротронные и нейтронные исследования 2008, №6, С. 33-38

84. Egorov V., Nadiradze A., Shaposhnikov V., Arbatskiy V. Ion Beam Study the Angle Distribution of Elements Contamination's in Xenon Jet of Electric Rocket Engine // IBA-14, Drezden, 26-30 Julay, 1999.

85. В.К Абгарян, Ю.А. Рыжов. Прогнозная модель расчёта эрозии стенок канала разрядной камеры СПД и ионно-плазменных источников// Вестник Московского авиационного института, 2008, Т. 15 № 5, С.40-43.

86. Баранцев Р.Г. Взаимодействие разреженных газов с обтекаемыми поверхностями., М. Мир, 1975 г., 255с.

87. Аксенова О.А. Фрактальная модель шероховатости поверхности взаимодействующей с разреженным газом // Математическое моделирование, 2001,т. 13, номер 7, сс.99-103.

88. Топорец А.С. Оптика шероховатой поверхности // М. Машиностроение, 1988, 191с.

89. A.Nakai, T.Aoki, TSeki, J.Matsuo,G.Takaoka, I.Yamada. Modeling of surface smoothing process by cluster ion beam irradiation///Nucl. Inst. And Meth. Phys. Res., B206(2003)c.842-845

90. Звенигород, 25-29 августа 2003 г. M.: 2003, Т.2, с.373-376.