Компьютерное моделирование извлечения ионов из плазмы Текст научной статьи по специальности «Физика»

ISSNÜS6S-5SS6

НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2ÜÜ1, том 11, № 2, c. 4Ü-43

ОРИГИНАЛЬНЫЕ СТАТЬИ

УДК 519.688: 537.534.7 © А. В. Ворошилов, А. П. Щербаков

КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ИЗВЛЕЧЕНИЯ ИОНОВ ИЗ ПЛАЗМЫ

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

ВВЕДЕНИЕ ВЫЧИСЛИТЕЛЬНАЯ ПРОЦЕДУРА

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

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

д// дt + V д// дг + Г д// др = 0 (1)

и уравнения Пуассона для электростатического потенциала и

Ди = -4лв(пг - пе). (2)

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

пг = п0г | /а3Р И

пе = П0е • ехр[е(и - ир )/кт;] .

Здесь ир и п0е — соответственно среднее

значение потенциала и средняя концентрация электронов в области плазмы (ионообразования), Те — электронная температура, п0г — средняя концентрация ионов.

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

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

Ф = е{и - ир )/кТе ,

x = ГKe , q=n>oe,

где XDe = (kTe/Ane2n0e ) — электронный дебаев-ский радиус, получаем уравнение Пуассона (2) в виде

А Ф = -Ч + ехр(ф).

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

^0 = С^ + С2^2 + С3Ф3 + С4Ф4 - С5 [ехр(<?0 Ь ч],

где (р0 — значение потенциала в рассматриваемом узле, ф1,.. .,ф4 — значения потенциала в соседних узлах, с1,.,с5 — коэффициенты, зависящие от геометрии и структуры сетки.

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

,л т+1

Ф0 =

= (ф0 — і)'5 ЄХр(0 )+ С1ф1 + С2ф2 + С3ф3 + С4ф4 + С5Ч

1 + С5 ехр(ф0т)

Здесь Ф0т+1 — значение потенциала в рассматриваемом узле на (т + і) -й ньютоновской итерации, ф0т — на предыдущей итерации.

Как показали численные эксперименты, итерационная схема оказалась устойчивой и обеспечила сходимость при относительной точности 10-4 за

а

б

в

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

без учета объемного заряда (а), пг ~ пе = 7.5 х 1010см-3(б), пг ~ пе

= 1.5 х 1011 см-3 (в),

пе = 3.0 х 1011 см-3 (г). Приведена верхняя половина системы

г

п

40-50 ньютоновских итераций. При относительной точности 10-3 требовалось всего несколько итераций.

Разработанная вычислительная процедура,

алгоритмы и программа в настоящей работе используются для моделирования источников с экстракцией ионов из плазмы. Такие системы используются для формирования интенсивных ионных пучков [2]. Схема источника изображена на рис. 1. Он состоит из трех осесимметричных электродов, к которым приложены потенциалы и1 = 40 000 В, и2 = 0 и и3 = 4000 В.

ИНТЕРПРЕТАЦИЯ РЕЗУЛЬТАТОВ КОМПЬЮТЕРНЫХ ЭКСПЕРИМЕНТОВ

Область слева от первого электрода — область плазмы. Эта область является наиболее "тяжелой" в вычислительном отношении. Концентрации ионов и электронов здесь приблизительно одинаковы, и в правой части уравнения Пуассона (2) стоит разность двух больших и приблизительно равных величин. Поэтому в большинстве имеющихся алгоритмов моделирование движения ионов начинается с границы плазмы, которая определяется вне вычислительной процедуры исходя из некоторых модельных представлений. Используемая в настоящей работе вычислительная методика позволила моделировать движение ионов как в области плазмы, так и вне ее. Граница плазмы в этой методике получается автоматически как переходная область от области с практически постоянным потенциалом и высокой концентрацией электронов, равной кон-

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

Результаты моделирования представлены на рис. 1. Здесь изображены электроды, эквипотен-циали поля и трубки ионного тока для различных значений плотности плазмы и соответственно для различных значений плотности ионного тока. Рис. 1, а дает картину траекторий пучка без учета объемного заряда — это первая итерация решения уравнения Власова. Пучок имеет широкую угловую расходимость. Последующие рисунки представляют результаты, когда итерационный процесс сошелся с заданной точностью по распределению поля и объемного заряда. Учет объемного заряда приводит к кардинальным изменениям в структуре поля в области первого электрода и в форме пучка. Сгущение эквипотенциалей поля в области первого электрода (здесь они проведены с шагом 5 В) следует трактовать как границу плазмы. Ее размеры составляют несколько радиусов Дебая, значения которого для рассматриваемых плотностей плазмы составляют 6х10-3см (рис. 1, б), 4х10-3см (рис. 1, в) и 3х10-3см (рис. 1, г). Распределение потенциала поля на оси для случая пг ~ пе = 3.0х 1011 см-3 представлено на рис. 2 (кривая 1): в области границы плазмы начинается падение потенциала поля. Здесь же для сравнения представлено распределение лапласова (без учета объемного заряда) поля.

С увеличением плотности плазмы происходит

г —-

■*ч ч

’

• ч "" Ч

24 \1

\

\

\ '

\

\

\ ч

\ '' ^

4

\

'

\

,

\ ' ч |

\\ ■

Чх '

1 ' ^

44 - _ 1 -— == * * *= ---

Рис. 2. Распределение поля на оси системы для случаев: 1 — пг ~пе = 3.0х 1011 см-3, 2 — без учета объемного заряда

сдвиг границы плазмы вниз по пучку (рис. 1, б, в, г), что соответствует основным экспериментально наблюдаемым закономерностям извлечения ионов из плазмы [2].

ЗАКЛЮЧЕНИЕ

Разработанная вычислительная процедура обеспечивает сходимость процесса за 20-70 итераций по ионному объемному заряду (внешние итерации). При этом требуется несколько десятков итераций метода SOR. Все это обеспечивает приемлемое время счета (доли минуты на один цикл итераций) на персональных компьютерах среднего уровня.

СПИСОК ЛИТЕРАТУРЫ

1. Щербаков А.П. Компьютерная модель формирования пучка ионов в источниках с электронным ударом при повышенных давлениях с учетом биполярного объемного заряда // Научное приборостроение. 2000. Т. 10, № 2. С.35-42.

2. Физика и технология источников ионов / Ред. Я. Браун. М.: Мир, 1998. 496 с.

Институт аналитического приборостроения РАН,

Санкт-Петербург

Материал поступил в редакцию 13.04.2001.

SIMULATION OF ION EXTRACTION FROM PLASMA

A. V. Voroshilov, A. P. Shcherbakov

Institute for Analytical Instrumentation RAS, Saint-Petersburg

Computer simulation of ion extraction from plasma is performed on the basis of the earlier developed [1] techniques for modeling ionic and electronic processes with account of bipolar space charge and ion interaction with a neutral gas component.