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

Шелестов, Александр Сергеевич

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

Шелестов, Александр Сергеевич АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Петрозаводск МЕСТО ЗАЩИТЫ

2006 ГОД ЗАЩИТЫ

01.04.04 КОД ВАК РФ

Диссертация по физике на тему «Моделирование процесса зарядки пылевой частицы и установления межчастичного расстояния в плазме низкого давления»

Автореферат диссертации на тему "Моделирование процесса зарядки пылевой частицы и установления межчастичного расстояния в плазме низкого давления"

ПЕТРОЗАВОДСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

ШЕЛЕСТОВ АЛЕКСАНДР СЕРГЕЕВИЧ

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

НИЗКОГО ДАВЛЕНИЯ

Специальность 01.04.04 - физическая электроника

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

Петрозаводск 2006

Работа выполнена на кафедре электроники и энергетики физико-технического факультета Петрозаводского государственного университета

Научный руководитель: доктор физико-математических наук,

профессор, Сысун В. И.

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

доктор физико-математических наук Мольков С.И. кандидат физико-математических наук Филиппов A.B.

Ведущая организация: Институт теплофизики экстремальных состояний Объединенного института высоких температур РАН, Москва

Защита диссертации состоится «_»_ г. в_часов на

заседании диссертационного совета К 212.190.01 при Петрозаводском государственном университете по адресу: Петрозаводск, пр. Ленина, 33, отзывы направлять по адресу:

185640, Петрозаводск, пр. Ленина, 33, ПетрГУ, ученому секретарю диссертационного совета К 212.190.01.

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

Автореферат разослан «_»_2006 г.

Ученый секретарь

диссертационного совета К 212.190.01 кандидат физико-математических наук

Алешина Л. А.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

являлось бы научно-обоснованной базой для объяснения имеющихся опытных данных.

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

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

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

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

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

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

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

Апробация работы: Основные результаты диссертационной работы были доложены на: заочной электронной конференции «Современные наукоёмкие технологии», проводимой Российской академией естествознания (РАЕ), 20-25 февраля 2005г., заочной электронной конференции «Фундаментальные исследования», проводимой РАЕ, 20-25 февраля 2005г., Второй международной научно-практической конференции «Ис-

следование, разработка и применение высоких технологий в промышленности» 7-9 февраля 2006 г., Санкт-Петербург, Россия., конференции «Фундаментальные и прикладные проблемы современной физики», Демидовские чтения, Москва-Екатеринбург-Томск, Россия, 25 февраля — 6 марта 2006 г., Второй всероссийской конференции молодых ученых «Физика и химия высокоэнергетических систем», 4-6 мая 2006 г., Томск, Россия., IV конференции «Фундаментальные и прикладные исследования. Образование, экономика и право», Римини, Италия, 9—16 сентября 2006 г., заочной электронной конференции «Математическое моделирование», проводимой Российской академией естествознания (РАЕ), 15-20 сентября 2006г,

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

Вклад автора. Исследования проведены в период 2003 — 2006 г. при непосредственном участии автора, которым сформулированы и обоснованы задачи диссертации. Разработка математической модели проведена при содействии Сысун А. В. Результаты раздела 2 были получены при содействии Подопригора А. В. В коллективных работах автору принадлежат изложенные в диссертации защищаемые положения. Все численные расчеты и написание текстов программ проведены автором.

Обьем и структура диссертации. Диссертация состоит из введения, 5 разделов, заключения, списка литературы и трех приложений. Диссертация содержит 108 стр., включая 25 рисунков, 10 таблиц и 55 наименований библиографических ссылок.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

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

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

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

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

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

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

В ряде работ для этого предлагается использовать метод молекулярной динамики. В работе [I] Р1С-МС моделированием в ячейке вблизи частицы определялся потенциал и заряд пылевой частицы в аргоне с давлением 0,1 тор. И хотя в широком диапазоне изменялась концентрация плазмы и, соответственно, дебаевский радиус, фиксация межчастичного расстояния не позволила обнаружить его влияние на заряд и потенциал частицы и, тем более, определить критерий его установления, В работе [2] методом молекулярной динамики определялся заряд и потенциал пылевой частицы в плазме неона. Радиус расчетной сферической ячейки был фиксирован и выбирался достаточно большим, чтобы тепловой хаотический поток ионов извне на ячейку существенно превышал уход ионов на пылевую частицу. Ионизация в объеме ячейки не учитывалась, поэтому такое приближение более подходит для одиночной пылевой частицы в плазме.

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

Исходя из этих выводов, сформулированы следующие задачи работы.

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

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

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

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

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

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

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

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

Раздел 2 содержит описание одномерной модели в сферической системе координат и результаты численного моделирования процесса зарядки пылевой частицы в плазме низкого давления путем прямого интегрирования уравнения Пуассона с учетом ионизации в межчастичной области. Расчет производился в ячейке Зейтца-Вигнера, которая представлялась сферической с радиусом г^, определяемым концентрацией частиц га=(4/3я п^'173. При кубической решетке пылевого кристалла г0 = 0,62Л, где Ь — межчастичное расстояние. При гранецентрирован-

ной решетке, чаще всего наблюдаемой в экспериментах, г0 = 0,5 9/г.

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

В приближении нулевой температуры ионов в режиме свободного пролета ионов уравнение Пуассона (1), представленное аналогично Ленгмюру в виде уравнения "плазма-слой" [4], решалось численно с частными аналитическими решениями на границе ячейки Зейтца-Вигнера и тонких сеточных слоях в промежутке.

1 8 д<р. __ "Т

г дг дг

еп

е<!

ехр

е<р

г'гЫг

' /

ехри г ! /27, ;

(о

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

'Л

|4лгг'2 «Дг'^с/г* (2),

1. =

1Ш

4 V лт

кТ 1

(3)

Нормировкой координаты г на дебаевский электронный радиус и потенциала на электронную температуру Те получен безразмерный вид уравнения (1), что показывает что ^ и Те, совместно с атомным весом иона являются определяющими параметрами плазмы.

Безразмерные величины:

координата X ~-=

Л/

тт е<? , потенциал и —-,

кТ

концентрация электронов и ионов Пе =

где пеа- концентрация электронов на границе ячейки Зейтца-Вигнера.

Скорость ионов V1 —

ЦК V м

, частота ионизации А —

гX.

___

л]кТе / А/

где со; - ионная плазменная частота 0)( —

е п

еО

у£0кТе ;

Тогда уравнение (1) будет иметь вид (Г)

д2и 2 ди , ,

дххдх

где выражения для концентрации электронов и ионов будут следующими (4):

л = ехр «У), п = — Г =тт • (4)

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

Для большей устойчивости численного решения использовались параболическая интерполяция потенциала по трём точкам (5):

тт(лЛ_(Ь-х)(2}г-х) х(2Ь-х) х(И-х) 2А2 + Ьг и1

где Ь — шаг дискретизации, х — отсчитывается от точки хо ,а также приближённое аналитическое решение для концентрации на каждом шаге (6):

• ■ £ о I х ,х

Г¥+ьй+е

хк У**

где а=Ц-2Ц+1 +Ц+2, Ь=-3 Ц+4 Ц+,- Ц+2, с=2 Ц-2ик, причём

<1х -1 (2 ах+Ь)

- — -агс51п

I-

т/ах2 +Ьх+с (Ь2 -4ас)"2

Результаты численных расчётов для распределения потенциала в разряде паров № для ряда значений г/Л^ (отношения радиуса зонда-пылинки к электронному дебаевскому радиусу) и гУЯ,/ (отношения размера ячейки Зейтца-Вигнера к электронному дебаевскому радиусу) приведены на рисунке 1.

-л»

-в*

ог

ЙЗ

/ /у / X

1 1 /

• 11/ V -

}

11111 I I I I

ГГ.

На всех кривых ход потенциала хорошо описывается приближенным аналитическим выражением (7);

е.. л2

<р{г) = <ра~ г

чО"«/

(7)

,где а- радиус пылинки, (ра — ее потенциал, а г^-размер ячейки

ЗеЙтца-Вигнера.

Одновременно рассчитывались полный заряд электронов (8) в ячейке, приведенный к заряду дебаевской сферы, заряд пылевой частицы (9), определенный по градиенту потенциала на ее поверхности, и общая потенциальная энергия электрического поля в ячейке (10), нормированная на полную тепловую (кинетическую) энергию электронов в дебаевской сфере:

4 я$пе(г)г2дг а'=3|и/(ж)*2& = —- (8)

д<р

дх

Я.

(9)

^-Г)ЕгАпггФ

**-им?е — <10>

2 * 2 &

Показаны существенные зависимости этих величин от размеров ячейки Зейтца-Вигнера. и радиуса пылинки

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

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

Прямое интегрирование уравнения Пуассона (1) возможно не всегда. Помимо условия нулевой температуры ионов Т. — О ограничением такого решения являлось появление неустойчивостей при больших от-

ге/

ношениях у , где а — радиус частицы, вследствие неопределенности

правой части выражения для концентрации ионов (4) при близком к нулю значении знаменателя.

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

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

их отношения ч/т. При этом не изменяются уравнения движения (11) и, соответственно, их решения.

^ = МУпГ = со)Ш. (11)

<Л М дг г

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

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

р(Я) = ~ехр ло

Л у

(12)

г« Л --1- - средняя*,и„а пробега «г-сечение—ения,

О^о

У0 -концентрация атомов).

Концентрация электронов принималась распределенной по Больц-ману (4). В контрольных численных расчетах использовалось более точное выражения для концентрации электронов (13), [5], учитывающего уход электронов на частицу:

Ни, -иг) ) кТ I 'г1

г /

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

был малым, \е<Ра \ — показал на близость Больцмановского распределения даже в этом случае. Таким образом распределение Больцмана достаточно точно описывает концентрацию электронов и его использование допустимо для дальнейших численных расчетов.

Электрическое поле определялось на каждом временном шаге решением уравнения Пуассона на сетке, по радиусу состоящей из N=200 - 500 интервалов, в зависимости от радиуса расчетной области При этом компьютерная переменная времени нормировалась на ионную плазменную частоту/' = Г -¿У,. Решение начиналось от внешней границы с граничными

условиями <ры ~ 0,| | =0. Потенциал пылевой частицы <ра определялся в конце счета. Заряд частицы находился согласно (9), он оказывался равным разности суммарных зарядов ионов и электронов в расчетной области с небольшими статистическими отклонениями, что указывает на правильность решения уравнения Пуассона.

Процесс счета продолжался до стабильного состояния системы, критерием которого выступало равенство ионного (2) и электронного (3) токов на пылинку.

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

Рассмотрим критерий установления межчастичного расстояния в пылевой плазме. Для системы заряженных частиц в качестве критерия устойчивого состояния используется минимум полной энергии электрического поля (14):

С другой стороны, критерием устойчивого положения заряженной частицы яаво внешнем электрическом поле с потенциалом <рвнещ является минимум потенциальной энергии этого заряда (15):

(15)

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

где га - радиус частицы.

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

здесь Е - напряженность поля на поверхности частицы, причем С ~ £0Еп, где Е„ - нормальная составляющая напряженности. Для проводящей частицы ЕиЕп совпадают, а, ввиду малой поверхности пылевой частицы (га~10'6 м.), это условие можно применять и для непроводящей частицы, из - за поверхностного перераспределения элементарных зарядов при возникновении даже малейшей разности потенциалов. За критерий установления межчастичного расстояния можно принять максимум значения этой поверхностной силы, действующей на частицу со стороны электрического поля. Другая частица будет притягиваться к данной до установления максимума силы. Смещение пылевой частицы в

(16):

(16),

(17),

обе стороны приведет к возникновению результирующей силы, возвращающей ее обратно.

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

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

а V Т I

— = 0,01 -г-1; = 0,05 -г 5; 0,01,0.1; — = 1, 10, 1000, где

Л/ Те хл

I екТе

— I-— - электронный дебаевский радиус,—длина пробега иона

до столкновения с атомом.

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

Значения потенциала и заряда пылинки, полученные по результатам численного эксперимента при нулевой температуре ионов в зависимости от межчастичного расстояния для различных радиусов пылинки, представлены на рисунке 2. Отдельные точки справа на графиках — расчеты, произведенные по радиальной теории, подобно тому, как это было проделано в работе [6]. Также была построена зависимость плавающего потенциала пылинки от ее радиуса в рамках радиальной теории для Ые, аналогично тому, как это было проделано в работе [6] для Аг и Н. Как видно из рисунков, радиальная теория является предельным случаем предложенной модели.

Заряд пылинки хорошо описывается выражением (19):

= 3д <ра Г \ О+Я]

Потенциальная энергия пылинки в поле окружающей ее плазмы можно приближенно рассчитать по (20):

(га-а)

а сила, действующая на пылинку, описывается согласно (21):

Л*' Л

2а

(21)

ООО 1,00 1.М 3.00 400 3,00 «.00 7,СО «00

Рис.2 17

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

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

Сил*

Рис.3

Подобным образом были получены зависимости потенциала, заряда и силы от размера ячейки, рассчитанные при отличном от нуля отношении ионной и электронной температур 7] / !ГС. При близких положениях максимума силы увеличивается потенциал и заряд пылевых частиц.

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

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

В плазме с фрактальными частицами (пыль) 11асто образуются . пылевые структуры в виде нитей (цепочек) вытянутых вдоль электрического поля. Часто такие цепочки наблюдаются на окраинах трехмерных плазменных пылевых структур (Рисунок 5).

-< *

1».»

Рис. 5

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

Рис. 6

Потенциал описываем двумерным уравнением Пуассона в цилиндрической системе координат (22):

д2<р 1 д<р д2<р е . ч с граничными условиями (23)

9[ша-0: =(^)»,/=0,' = =0(23,

дг у 5/* ш &

Потенциал определялся методом итераций (100-200) на каждом временном шаге. Ионы движутся в электрическом поле согласно уравнениям движения (24). В цилиндрической системе имеем три компоненты скорости У2,

дК Э/

е дф , V* дУ2

М дг

е д(р_ М дг

(24)

Так как в силу симметрии = 0, то угловая скорость определяется условием сохранения углового момента.

Полное число макроинов в ячейке достигало десятков миллионов. Расчетное время одного режима (одна точка) в двумерном случае составляло несколько суток.

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

Эарцд

Рис. 7

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

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

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

Рис. 8

Отдельными точками на графике зависимости межчастичного расстояния от размеров пылинки показаны экспериментальные данные.

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

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

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

3. . Показано, что при заданных параметрах плазмы и размерах частиц существует установившийся размер ячейки Зейтца-Вигнера.

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

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

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

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

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

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

Основные результаты опубликованы в виде статей и материалов докладов конференций:

1. Сысун В.И., Хахаев А.Д., Олещук О.В., Шелестов A.C. Заряд и потенциал пылевой частицы в плазме низкого давления с учётом ионизации в области возмущения // «Физика плазмы» 2005, т.31, №9, с. 834-841.

2. Сысун В.И., Хахаев А.Д., Олещук О.В., Шелестов A.C. Формирование ионного потока на пылевую частицу в плазме // «Современные наукоёмкие технологии» 2005, №3, с. 30-34.

.3. Шелестов A.C., Подопригора A.B. Моделирование ионного потока на пылевую частицу в плазме методом «крупных частиц в ячейке» // «Фундаментальные исследования» 2005, №3, с. 60-62.

4. Шелестов A.C., Сысун A.B., Двумерная модель заряда пылевой частицы в плазме низкого давления // сборник трудов «Высокие технологии, фундаментальные и прикладные исследования, образование», изд -во Политехнического университета, том 5 стр. 310-312, материалы Второй международной научно - практической конференции «Исследование, разработка и применение высоких технологий в промышленности» 7-9 февраля 2006 г., Санкт-Петербург, Россия.

5. Шелестов A.C., Сысун В.И., Численное моделирование параметров низкотемпературной плазмы методом крупных частиц // Материалы конференции «Фундаментальные и прикладные проблемы современной физики», Демидовские чтения, Москва-Екатеринбург-Томск, Россия, 25 февраля — 6 марта 2006 г., стр. 21 - 22.

6. Шелестов A.C., Модель расчета параметров кристаллических структур пылевой плазмы низкого давления // Материалы Второй всероссийской конференции молодых ученых «Физика и химия высокоэнергетических систем», 4-6 мая 2006 г., Томск, Россия, стр. 338 - 341

7. Сысун A.B., Шелестов A.C., Критерий установления межчастичного расстояния в пылевой плазме // Материалы IV конференции «Фундаментальные и прикладные исследования. Образование, экономика и право», Римини, Италия, 9-16 сентября 2006 г.

8. Sysun А. V., Shelestov A.S., Criterion of intergrain distance establishing in dusty plasma // European Journal of Natural History, 2006, №5, pp. 8688.

9. Шелестов A.C., Сысун A.B., Моделирование процесса зарядки пылевой частицы в плазме низкого давления методом молекулярной динамики // Материалы заочной электронной конференции «Математическое моделирование», проводимой Российской академией естествознания (РАЕ), 15-20 сентября 2006г.

ЦИТИРОВАННАЯ ЛИТЕРАТУРА

[1] J. P. Baeuf, Physical Review A., v.46, №12, p.7910, 1992.

[2] A.B. Зобнин, А.П. Нефедов, B.A. Синелыцнков, B.E. Фортов, ЖТЭФ, 2000, т. 118, вып 3(9), стр. 554-559.

47-67.

[3] А.Ф. Паль, Д.В, Сивохин, А.Н. Старостин, A.B. Филиппов, В.Е. Фортов, Физика плазмы, 2002, т.28, №1» стр. 32-44.

[4] Tonks L., Langmnuir I. Phisycal Review. 1929. v.34. p.876-922.

[5] Lam S.H. The Langmuir Probe in a Collisionless Plasma, Physics Fluids, v.8, pp. 73-87, 1965.

[6] Kennedy R.V., Allen J.E. J. Plasma Phisics. V.67. 2002. pp.243-250.

Подписано в печать 16.11.2006 Формат 60x84 1/16 Бумага офсетная Уч-изд. л. 1. Тираж 100 экз, Изд. № 268

Государственное образовательное учреждение высшего профессионального образования ПЕТРОЗАВОДСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

Типография Издательства ПетрГУ 185910, г.Петрозаводск, пр.Ленина, 33

Содержание диссертации автор исследовательской работы: кандидата физико-математических наук, Шелестов, Александр Сергеевич

ВВЕДЕНИЕ.

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

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

1.1.1 Методы теории зондов применительно к пылевой плазме.

1.1.2 Теория радиального дрейфа.

1.1.3 Теория ограниченных орбит.

1.1.4 Гидродинамическая модель диффузионного приближения.

1.2 Обзор экспериментальных работ по плазменно-пылевым структурам.

1.2.1 Эксперименты по определению межчастичного расстояния.

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

1.2.3 Основные силы, действующие на пылевые частицы в плазме.

1.3 Обзор работ по математическому моделированию.

1.4 Выводы по литературному обзору и постановка задач.

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

2.1 Описание метода.

2.2 Результаты численных расчётов.

2.3 Условия применимости и трудности метода.

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

3.1 Описание метода молекулярной динамики.

3.2 Моделирование максвелловского распределения ионов по скоростям.

3.3 Концентрация электронов вблизи отталкивающего сферического зонда с учетом ухода электронов на зонд.

3.4 Критерий установления межчастичного расстояния.

3.5 Исходные уравнения.

3.6 Безразмерные переменные.

3.7 Алгоритм вычислений.

3.8 Результаты расчета.

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

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

4.2 Исходные уравнения.

4.3 Интегральные параметры процесса.

4.4 Безразмерные переменные.

4.5 Взвешивание и раздача заряда и полей в цилиндрических координатах.

4.6 Максвелловское распределение. Цилиндрические координаты.

4.7 Линейное взвешивание по концентрации.

4.8 Алгоритм двумерного метода молекулярной динамики.

4.9 Анализ погрешности типов взвешивания при различных радиальных изменениях концентраций (цилиндрический случай).

4.10 Результаты моделирования.

5. СРАВНЕНИЕ РЕЗУЛЬТАТОВ ЧИСЛЕННОГО ЭКСПЕРИМЕНТА С ОПЫТНЫМИ ДАННЫМИ.

5.1 Итоговые зависимости.

5.2 Экспериментальные данные.

Введение диссертация по физике, на тему "Моделирование процесса зарядки пылевой частицы и установления межчастичного расстояния в плазме низкого давления"

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

Пылевая плазма представляет собой ионизованный газ, содержащий пылинки — частицы твёрдого вещества. Такая плазма широко распространена во Вселенной: в планетных кольцах, хвостах комет, межпланетных и межзвёздных облаках, вблизи искусственных спутников Земли, в пристеночной области термоядерных установок с магнитным удержанием, а также в плазменных реакторах, дугах, разрядах. Экспериментально в лабораторных условиях пылевую плазму впервые получил Ирвинг Лэнгмюр в 20-х годах прошлого века [1]. Однако, активно изучать её начали лишь в последнее время. Повышенный интерес к свойствам пылевой плазмы возник с развитием технологий плазменного напыления и травления в микроэлектронике, а также производства тонких плёнок. Наличие твёрдых частиц, которые попадают в плазму в результате разрушения электродов и стенок разрядной камеры приводит к загрязнению поверхности полупроводниковых микросхем и возмущает плазму часто непредсказуемым образом. Чтобы предотвратить или хотя бы уменьшить эти негативные явления, необходимо разобраться в том, каким образом идут процессы образования и роста конденсированных частиц в газоразрядной плазме, как плазменные пылинки влияют на свойства разряда, что именно приводит к их структурному упорядочиванию.

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

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

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

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

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

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

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

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

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

3. Показано, что при заданных параметрах плазмы и размерах частиц существует установившийся размер ячейки Зейтца-Вигнера.

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

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

- Заочной электронной конференции «Современные наукоёмкие технологии», проводимой Российской академией естествознания (РАЕ), 20-25 февраля 2005г.

- Заочной электронной конференции «Фундаментальные исследования», проводимой РАЕ, 20-25 февраля 2005г.

- Второй международной научно-практической конференции «Исследование, разработка и применение высоких технологий в промышленности» 7-9 февраля 2006 г., Санкт-Петербург, Россия.

- Конференции «Фундаментальные и прикладные проблемы современной физики», Демидовские чтения, Москва-Екатеринбург-Томск, Россия, 25 февраля - 6 марта 2006 г.

- Второй всероссийской конференции молодых ученых «Физика и химия высокоэнергетических систем», 4-6 мая 2006 г., Томск, Россия.

- IV конференции «Фундаментальные и прикладные исследования. Образование, экономика и право», Римини, Италия, 9-16 сентября 2006 г.

- Заочной электронной конференции «Математическое моделирование», проводимой Российской академией естествознания (РАЕ), 15-20 сентября 2006г. и опубликованы в виде статей и тезисов докладов конференций;

1. Сысун В.И., Хахаев А.Д., Олещук О.В., Шелестов А.С. Заряд и потенциал пылевой частицы в плазме низкого давления с учётом ионизации в области возмущения // «Физика плазмы» 2005, т.31, №9, с. 834-841, Plasma Physics Reports, Vol. 31, № 9,2005, pp.772-779.

2. Сысун В.И., Хахаев А.Д., Олещук О.В., Шелестов А.С. Формирование ионного потока на пылевую частицу в плазме // «Современные наукоёмкие технологии» 2005, №3, с. 30-34.

3. Шелестов А.С., Подопригора А.В. Моделирование ионного потока на пылевую частицу в плазме методом «крупных частиц в ячейке» // «Фундаментальные исследования» 2005, №3, с. 60-62.

4. Шелестов А.С., Сысун А.В., Двумерная модель заряда пылевой частицы в плазме низкого давления // сборник трудов «Высокие технологии, фундаментальные и прикладные исследования, образование», изд - во Политехнического университета, том 5 стр. 310-312, материалы Второй международной научно - практической конференции «Исследование, разработка и применение высоких технологий в промышленности» 7-9 февраля 2006 г., Санкт-Петербург, Россия.

5. Шелестов А.С., Сысун В.И., Численное моделирование параметров низкотемпературной плазмы методом крупных частиц // Материалы конференции «Фундаментальные и прикладные проблемы современной физики», Демидовские чтения, Москва-Екатеринбург-Томск, Россия, 25 февраля - 6 марта 2006 г., стр. 21-22.

6. Шелестов А.С., Модель расчета параметров кристаллических структур пылевой плазмы низкого давления // Материалы Второй всероссийской конференции молодых ученых «Физика и химия высокоэнергетических систем», 4-6 мая 2006 г., Томск, Россия, стр. 338 - 341

7. Сысун А.В., Шелестов А.С., Критерий установления межчастичного расстояния в пылевой плазме // Материалы IV конференции «Фундаментальные и прикладные исследования. Образование, экономика и право», Римини, Италия, 9-16 сентября 2006 г.

8. Sysun А.V., Shelestov A.S., Criterion of intergrain distance establishing in dusty plasma // European Journal of Natural History, 2006, №5, pp.86-88.

Вклад автора. Исследования проведены в период 2003 - 2006 г. при непосредственном участии автора, которым сформулированы и обоснованы задачи диссертации. Разработка математической модели проведена при содействии Сысун А. В. Результаты раздела 2 были получены при содействии Подопригора А. В. В коллективных работах автору принадлежат изложенные в диссертации защищаемые положения. Все численные расчеты и написание текстов программ проведены автором.

Обьем и структура диссертации. Диссертация состоит из введения, 5 разделов, заключения, списка литературы и трех приложений. Диссертация содержит 108 стр., 28 рисунков, 10 таблиц и 55 наименований библиографических ссылок.

Заключение диссертации по теме "Физическая электроника"

Выводы

1) При линейном и квадратичном уменьшении концентрации с радиусом ближе всего линейное взвешивание по концентрации

2) При спаде концентрации обратно-пропорционально радиусу - точнее всего линейное взвешивание по заряду.

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

4.10 Результаты моделирования

Результаты моделирования двумерного цилиндрического случая для Te/Tj=0, при различных размерах пылинки и ячейки Зейтца-Вигнера приведены в таблице 9.

ЗАКЛЮЧЕНИЕ

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

Из проведенных численных экспериментов и теоретического анализа вытекают следующие основные результаты диссертационной работы:

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

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

Список источников диссертации и автореферата по физике, кандидата физико-математических наук, Шелестов, Александр Сергеевич, Петрозаводск

1. Langmuir 1., Mott-Smith H.M., Gen.Electr.Rev.,1924,v.27, p.p. 449,538,616,762,810; Phys.Rev., 1926,v.26,p.727-763.

2. Bohm D.,ed. H.Guthric and R.K.Wakering. New York.: McCraw-Hill. 1949.

3. Kennedy R.V., Allen J.E. J. Plasma Phisics. V.67. 2002. pp.243-250.

4. Allen J.E., Annaratone B.M., U. de Angeles J. Plasma Phisics. V.63. 2000. p.299.

5. Bernstain I., Rabinowitz I. Phys. Fluids, 1959. V.2. P.l 12-121.

6. Laframboise J.G. Rarefield Gas Dynamics. New York, Academic Press, 1966. V.II. P22. Univ. of Toronto, Institute of Aerospace Studies, 1966. Report 100.

7. Hall L.S., Fries R.R. Proc. 7th Int. Conf. Phen. Ion. Gases. 1966. V.3. P.15-19.

8. Kennedy R.V., Allen J.E. J. Plasma Phisics. V.69. 2003. pp.485-506

9. Швейгерт B.A., Швейгерт И.В., Богданов B.M. и др. ЖЭТФ. Т.115. 1999. с.877

10. Зобнин А.В., Нефёдов А.П., Синельников В.А., Фортов В.Е.ЖЭТФ. т. 118. 3(9) 2000 с. 554-559

11. Паль А.Ф., Старостин А.Н., Филиппов А.В. Физика плазмы. Т.27. 2001 2.с.155-164. т.28. 2002. №1. с.32-44.

12. Morfill G. Е. "Plasma Crystal" Preprint Max Planck Institute for Extraterrestrial Physics Germany, Garching, 1994

13. Thomas H, Morfill G et al. Phys. Rev. Lett. 73 5 652 (1994)

14. Chu J H, Lin I, Phys. Rev. Lett. 72 25 4009 (1994)

15. Allen J et al. Proc. Ann. Plasma Phys. Conf. (Scotland: Pearth, 1996)

16. Chu J.H., Lin I, Physica A 205 183 (1994)

17. Fortov V E, Nefedov A F Phys. Lett. A 218 89 (1996)

18. Fortov V.E., Nefedov A.P., Torchinskij V.M., Molotkov V.I.,Khrapak A.G., Petrov O.F., Volikhin K.F., J.E. T. F. Lett., 64, 2 25 (1996)

19. Khakhaev A. D., Luizova L. A., Piskunov A. A., Podryadchikov S. F., Pushcarev A. V., Shtykov A. S., XVI International Conference on Gas Discharges and their Applications, Xi'an, China, September 11-13, 2006. 2006. V. 1. P. 341— 344.

20. Melzer A Homann, A. Piel Phys. Rev. E. v. 53 3 2757 (1996)

21. Walch В, Horanyi M., Robertson S., Phys.Rev. Lett., 75, 838, (1995)

22. Sitkafoose A. A., et. al., Phys.Rev. Lett., 84, 6034, (2000)

23. Фортов B.E. и др.ДЭТФ 114, 2004, (1998)

24. Fortov V.E. et al, Phys.Rev. Lett., A 229 317 (1997)

25. Нефедов А.П. и др., письма в ЖЭТФ, 72 313 (2000)

26. Barkan A., D'Angelo N., Merlino R., Phys.Rev. Lett., 73, 3093, (1994)

27. Fortov V.E. et al, Phys. Rev. Lett., 87 205002 (2001)

28. Фортов B.E., Храпак C.A., Храпак А.Г., Молотков В.И., Петров О.Ф., УФН, т.174, 5 (2004)

29. Цытович В.Н., УФН, т.167 1 (1997)

30. Filinov V. S., Petrov О. F., Fortov V. Е., Molotkov V. I., Khahaev A. D., Podryadchikov S. F., Contrib. Plasma Phys. 2005. V. 45. № 3-4. P. 177-185.

31. Игнатов A M Физика плазмы 27 323 (1996)

32. Morfill G, Khodataev Ya K, Tsytovich V N Comm. Plasma Phys. Contr. Fusion 17 342 1996)

33. Лифшиц E.M., Питаевский Л.П., Физическая кинетика, М.Наука (1979)

34. Epstein P.S, Phys.Rev, 23 710 (1924)

35. Draine B.T., Salpeter E.E., Astrophys J., (1979)

36. Nitter Т., Plasma Sources Sci.Technol. 5, 93, (1996)

37. Daugherty J.E., Porteous R.K., Graves D.B., J.Appl.Phys., 73, 1617, (1993)

38. Hamaguchi S., Farouki R. Т., Phys. Rev.,E 49, 4430 (1994)

39. Hamaguchi S., Farouki R. Т., Phys. Plasmas, 1, 2110 (1994)

40. Ходатаев Я. К., Бингхем Р., Тараканов В. П., Цытович В. Н., Физика плазмы. 1996. №22. С. 1028

41. Winske D., Murillo М., Selwyn G., Bull. Amer. Phys. Soc., 1997.

42. Otani N., Bhattacharjee A., Phys. Rev. Lett. 1997. № 78. P. 1468.

43. Nairn C.M.C., Annaratone B.M., Allen J.E. Plasma Sources Sci. Technol. V.7. 1998. p.478.

44. Иньков Л. В., Левченко В. Д., Сигов Ю. С., Тез. докл. XXV Звенигородской конф. по физике плазмы и УТС, 178, Звенигород, 2—6 марта, 1998

45. Иньков Л. В., В. Д. Левченко, 10. С. Сигов., Прикладная Физика, 2000 3, с 138-145.

46. J. P. Baeuf, Physical Review A., v.46, 12, p.7910, 1992.

47. A.B. Зобнин, А.П. Нефедов, B.A. Синельщиков, B.E. Фортов, ЖЭТФ, 2000, т. 118, 3(9), стр. 554-559

48. Ваулина О.С., Репин АЛО., Петров О.Ф., Физика плазмы, 2006, т.32, 6, с.528 531.

49. Tonks L., Langmnuir I. Phisycal Review. 1929. v.34. p.876-922

50. Lam S.H. Phys. Fluids, 1965. V.8. P. 1002. P.73-87

51. Allen J.E., Boyd R.L., Reynolds R., Proc. Phys. Soc. 1957. V.B70. P.297-304.

52. Мартынов H. H. Matlab 7. Элементарное введение. M., 2005. 416 c.

53. Swift I.D., Schwar M. JR., Electric probes for plasma diagnostics, Iriffe Books, pp. 321, London, 1971.

54. P. Бэдсел, А. Ленгтон, «Физика плазмы и численное моделирование», Москва, Энергоатомиздат, 1989 г., 456 с

55. Р. Хокни, Дж. Иствуд, «Численное моделирование методом частиц», Москва, Энергоатомиздат, 1987 г., 640 с.