ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ
«БЕЛГОРОДСКИЙ ГОСУДАРСТВЕННЫЙ НАЦИОНАЛЬНЫЙ
ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ»
(НИУ «БелГУ»)
ИНСТИТУТ ИНЖЕНЕРНЫХ ТЕХНОЛОГИЙ И ЕСТЕСТВЕННЫХ НАУК
КАФЕДРА ОБЩЕЙ МАТЕМАТИКИ
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ЭЛЕКТРОННОЙ ПУШКИ
С УПРАВЛЕНИЕМ КОНФИГУРАЦИЕЙ ПУЧКА
Выпускная квалификационная работа
обучающегося по направлению подготовки
01.03.02 Прикладная математика и информатика
очной формы обучения,
группы 07001305
Гончарова Дмитрия Викторовича
Научный руководитель
к.ф.-м.н.
Никуличева Т.Б.
БЕЛГОРОД 2017
Оглавление
Введение ................................................................................................................... 3
Глава 1. Теоретические основы предметной области ......................................... 8
1.1 Основное назначение электронной пушки .................................................. 8
1.2 Проектирование пушек, формирующих цилиндрические и плоские
пучки.................................................................................................................... 12
1.3 Процедура математической параметризации электронной пушки ........ 12
Глава 2. Математические методы анализа электронных пушек ...................... 20
2.1 Расчеты электронных пушек методом синтеза. ........................................ 20
2.2 Расчеты электронных пушек методом анализа. ....................................... 22
2.3 Способы расчета реальных магнитных полей в электронных пушках .. 26
2.4 Расчет электростатических электронных пушек ...................................... 30
Глава 3. Проведение модельных экспериментов по исследованию параметров
пушки ...................................................................................................................... 39
Основные результаты и выводы. ......................................................................... 49
Заключение ............................................................................................................ 50
Список использованных источников .................................................................. 51
2
Введение
Электронная пушка (ЭП) – это устройство, предназначенное для
получения пучка электронов определенной конфигурации, обладающих
определенным интервалом разброса значений кинетической энергии. ЭП
служат для генерации и формирования электронного пучка, иногда
называемого пучком. Чаще всего данное приспособление используется в
кинескопах и других электронно-лучевых трубках, СВЧ-приборах (например,
в лампах бегущей волны), а также в приборах, в основе которых лежит
использование вместо светового потока пучка электронов, например,
электронные микроскопы и ускорители заряженных частиц.
Типичными задачами моделирования ЭП являются задачи определения
характера движения электронов в потоке, формируемом электродами
заданной конфигурации, обычно без учета пространственного заряда. Путем
последующего изменения формы и расположения электродов добиваются
требуемых параметров электронного пучка. Часто желательно бывает
определить геометрические формы, расположение электродов и потенциалы
на них, считая известными физические параметры пучка.
Основными недостатками существующих устройств ЭП являются
ограничение функциональных возможностей, надежности и ресурса работы,
обусловленное наличием термоупругих напряжений в основных узлах,
обеспечивающих рабочий режим электронно-лучевой пушки, что затрудняет
их юстировку.
Актуальной задачей в настоящее время является создание компактных,
технологичных
и
удобных
в
эксплуатации
электронных
пушек,
формирующих протяженный электронный поток заданной ширины с
равномерным распределением плотности тока по его сечению и обладающих
высокой надежностью в течение длительного срока службы. Обеспечение
этих параметров электронной пушки во многом зависит от конструкции ее
катодного узла. Оптимальная геометрия катодного узла позволит повысить
3
равномерность распределения протяженного электронного потока по
сечению и повысить надежность и долговечность катодного узла.
В том случае, если при конструировании электронной пушки
используются массивные и громоздкие элементы устройства для крепления
катодной нити, происходит рассеивание тепловых потоков, уменьшается
температура на концах катодной нити, что приводит к уменьшению тока
эмиссии и удельной мощности линейного электронного потока и создает
неравномерное распределение плотности тока по длине катода. Кроме того,
при высокой рабочей температуре снижается предел прочности материала
катода, что ограничивает ресурс его работы. Также в случае деформации
крепежных элементов в результате нагрева может происходить искажение
геометрии ЭОС электронной пушки прибора и может нарушаться
токопрохождение в ней. Все это приводит к снижению надежности и
долговечности работы катода, катодного узла и электронной пушки в целом.
Катодный узел содержит ряд элементарных ленточных катодов,
расположенных последовательно вдоль общей оси, параллельной продольной
оси окна вывода электронов. Элементарные ленточные катоды закреплены с
помощью держателей элементарных ленточных катодов на общем основании
катодного узла, расположенном параллельно плоскости торцевого окна
вывода электронов. Обращенные в сторону окна вывода электронов
эмитирующие
поверхности
катодов
могут
иметь
плоскую
или
цилиндрическую (выпуклую, вогнутую) форму. Элементарные ленточные
катоды отделены друг от друга технологическими зазорами. При подаче тока
накала на катоды вследствие их теплового расширения происходит
увеличение линейных размеров катодов, в результате чего элементарные
катоды могут сближаться друг с другом.
При выборе оптимальной величины зазоров между элементарными
ленточными катодами можно получить при нагреве общую эмитирующую
поверхность,
близкую
к
непрерывной
поверхности,
и
обеспечить
4
формирование протяженного электронного потока с более равномерным
распределением плотности тока.
В реальных конструкциях катодных узлов длина общей эмитирующей
поверхности может составить несколько сотен миллиметров. Поэтому при
использовании в них проволочных держателей для крепления элементарных
ленточных катодов возникают технологические проблемы, связанные с
тепловыми деформациями элементов катодного узла. Под воздействием
механических усилий и высокой температуры проволочные держатели
деформируются и изгибаются во всевозможных направлениях, изменяя
положение закрепленных на них элементарных ленточных катодов друг
относительно друга, а также относительно других электродов электронной
пушки. В результате происходит искажение формы элементарных ленточных
катодов, катодного узла и геометрии ЭОС электронной пушки. При этом
искажается
форма
элементарных
отдельных
ленточных
электронных
катодов,
и
потоков,
общего
получаемых
электронного
с
потока
электронной пушки, происходит неравномерное смещение электронных
потоков относительно друг друга и электродов электронной пушки и
оседание электронов на этих электродах. В результате ухудшаются
равномерность
распределения
электронного
потока
по
сечению
и
токопрохождение, увеличиваются потери мощности, снижается КПД и срок
службы электронной пушки. Усилению этих недостатков способствуют
также имеющие место технологические разбросы параметров и геометрии
отдельных катодов и держателей.
Для практической реализации эффективных электронных пушек в
качестве источника электронов исследователи рекомендуют использование
полевого электронного катода (ПЭК), обладающего рядом неоспоримых
преимуществ по сравнению с другими видами источников свободных
электронов. К этим преимуществам можно отнести:
- отсутствие напряжения накала;
- отсутствие температурного градиента;
5
- устойчивость ключевых узлов к внешней радиации;
- безынерционность;
-
крутизна
(вольт-амперных
характеристик)
ВАХ
по
экспоненциальному закону.
Совокупность
этих
свойств
обусловливает
перспективность
использования ПЭК в различных электронных приборах таких, как
электронно-лучевые трубки, плоские дисплейные экраны и т.д. Однако, на
настоящий момент перечисленные преимущества ПЭК и их использование в
приборах вакуумной электроники сложны в практической реализации.
Экспериментальные методы исследования требуют определенных
экономических
исследования
численного
и
временных
являются
затрат,
методы
эксперимента
с
поэтому
основными
математического
применением
методами
моделирования
и
средств
современной
расчета
характеристик
компьютерной техники.
Теоретическое
исследование
в
области
электронных полевых пушек представлено в работах Добрецова Л.Н.,
Егорова Н.В., Шешина Е.П. Следует отметить, что в последние годы идет
интенсивный процесс создания новых устройств на основе ЭП.
Поэтому, тема представленной ВКР «Математическое моделирование
электронной пушки с управлением конфигурацией пучка», несомненно,
является актуальной.
Целью ВКР является исследование электронно-лучевой пушки с
изменяемой конфигурацией пучка заряда.
Исходя из поставленной цели, необходимо решить следующие задачи:
1. Изучить библиографические источники по теме ВКР
2. Рассмотреть аналитическое представление задачи моделирования
формы пучка электронов
3. Применить метод анализа для решения задачи
4. Применить метод синтеза для решения задачи
5. Построить модель электронно-лучевой пушки
6
6. Провести тестирование разработанной модели
Научная новизна работы
1. Математическое моделирование электрического поля в системе с
катодным узлом произвольной конфигурации.
2. Алгоритм
вычисления
электрического
поля
по
заданному
распределению потенциала с автоматическим выбором шага.
3. Прогнозирование структуры и свойств электронного пучка на основе
вычислительного эксперимента.
Практическая значимость работы
1. Разработан программный комплекс для моделирования электронной
пушки.
2. Возможность
использования
результатов
моделирования
для
практической реализации.
Структурно ВКР состоит из введения, трех глав, заключения, списка
используемых библиографических источников и приложений.
Во
введении
сформулирована
цель,
рассмотрена
актуальность
поставлены
задачи,
выбранной
приведено
темы,
содержание
пояснительной записки ВКР.
В
первой
главе
рассмотрена
предметная
область,
приведены
разновидности схем электронных пушек, проанализирована типичная
процедура проектирования пушки, проведен обзор математических пакетов и
программных средств для моделирования.
Во второй главе рассмотрены аналитический и синтетический методы
решения поставленной задачи, выбраны исходные типы данных, построены
зависимости.
В третьей главе приведена практическая реализация моделирования
электронной пушки и тестирование построенной модели.
В заключении приведены основные результаты и выводы по работе.
7
1. Теоретические основы предметной области
1.1 Основное назначение электронной пушки
Современные
электронно-лучевые
и
радиотехнические
приборы
высокой точности в основе своей работы используют различные по
интенсивности и пространственной конфигурации электронные пучки
(рис.1.1).
Рис. 1.1. Схема электронной пушки: 1 — катод; 2 — модулятор; 3 —
первый анод; 4 — второй анод; е — траектории электронов.
Одной из наиболее важных и сложных в техническом исполнении
задач, которые необходимо решить при изготовлении и практическом
использовании ЭП – это формирование и фокусировка интенсивных
электронных пучков. Основная задача электронной пушки заключается в
формировании электронного пучка определенной интенсивности, заданной
конфигурации,
с
управляемым
лучом,
что
достигается
заданием
управляющего напряжения, как следствие, регулировкой скорости движения
электронов (в идеальном случае – ламинарным) (рис.1.2).
Методы управления конфигурацией и фокусировкой электронных
пучков,
связаны
с
заданием
определенных
значений
управляющих
напряжений, силы тока, а также скорости электронов, имеющих, по
возможности, движение. Техническая реализация таких приборов оказалась
связанной с большим количеством проблем и вопросов, и изначально не дала
положительных результатов. Поэтому разработка электронных пушек тесно
связана с поисками новых формирующих и фокусирующих систем.
8
Рис. 1.2. Структурная схема осесимметричной электронной пушки,
используемой в клистронах (показана в разрезе).
Существенный прорыв в области проектирования электронных пушек
был сделан Дж. Р. Пирсом, который предложил совершенно новый метод
формирования прямолинейных ламинарных электронных пучков простой
конфигурации.
Данный
метод
был
положен
в
основу
разработки
электронных пушек, отличающихся высокой эффективностью, которые в
настоящее время широко известны под названием пушек Пирса (рис. 1.3).
Конструкция типичной пушки Пирса, аналог которой представлен на рис.1.3,
состоит из вогнутого сферического эквипотенциального катода 2 с
подогревателем 1, широко раскрытого прикатодного электрода 4 и анода 3 с
центральным отверстием.
Для того чтобы получить узкий электронный пучок, падающий через
отверстие анода определенного диаметра на регистрирующий экран,
требуется создать в камере такую конфигурацию электрического поля, при
которой электроны со всей поверхности катода равномерно сходятся в узкий
электронный пучок. Для этого необходимо рассчитать параметры формы
электродов. Наиболее эффективным, с точки зрения экономических и
временных затрат, методом такого расчета является аналитический метод
или метод математического моделирования.
9
Пушки
Пирса,
рассмотренной
ниже
конфигурации,
позволяют
получать узкие сходящиеся электронные пучки, обладающие малой
интенсивностью.
Рис. 1.3. Экспериментальная электронная пушка Пирса, Ua=50 кВ,
P=0.5 мкА/В3/2, Сj=100.
Чтобы получить сфокусированные электронные пучки, обладающие
более высокой интенсивностью существует очень широкое многообразие
модифицированных пушек Пирса (рис. 1.4). Их основными особенностями
являются относительно небольшие расстояния между катодом и анодом, а
также применение прикатодных электродов закрытой конструкции. Это
позволяет получить в электронных пушках электрические поля, которые
могут
скомпенсировать
поперечные
расфокусирующие
силы
пространственного заряда более эффективно. Но такая конфигурация
обладает рядом недостатков, например, при малых расстояниях катод-анод
между катодом и фокусирующим электродом отсутствует необходимый
тепловой зазор.
10
Рис.
1.4.
Экспериментальная
модифицированная
пушка
Пирса,
Ua=20кВ, Р=3 мкА/B3/2, Сj=6: 1—-подогреватель; 2—-катод; 3—прикатодный
электрод; 4— -анод; 5—коллектор.
Один из примеров математического расчета пучка с управляемой
конфигурацией профиля, показан на рис. 1.5. Здесь: 1 – катод, 2 –
экранирующий электрод, 3 – управляющий электрод, 4 – анод, 5 – выходной
электрод (заземленный). Недостатком данной модели является сильная
замагниченность пучка. Регулируемым параметром, позволяющим задавать
необходимую конфигурацию профиля пучка в данной модели является
управляющее напряжение на электроде 3, которое меняется от –2 кВ до +2
кВ в зависимости от напряжения на катоде.
Рис. 1.5. Геометрия электронной пушки с управляемым профилем
пучка.
11
1.2 Проектирование пушек, формирующих цилиндрические и плоские
пучки.
Для
математического
анализа
электронных
пушек
различной
конфигурации, необходимо сначала изучить физические принципы, которые
ложатся в основу разработки пушек Пирса. Например, очень важным
является тот факт, что пушки Пирса позволяют формировать не только
сплошные сходящиеся электронные пучки, называемые коническими, а
также цилиндрические, плоские (ленточные и клиновидные) и трубчатые
(полые) пучки. Установление формы конфигурации пучка является важной,
т.к. от этого зависит выбор расчетной системы координат.
Поэтому при моделировании электронной пушки, которая формирует
параллельный ленточный пучок, используется декартова система координат,
а модель пучка задается как часть плоскопараллельного потока.
При
моделировании
электронной
пушки,
которая
формирует
клиновидный, т.е. сходящийся пучок, используется цилиндрическая система
координат, а модель пучка задается как часть радиального цилиндрического
потока.
Чтобы сформировать модель параллельного цилиндрического пучка
можно поступить следующим образом: из бесконечного параллельного
потока «вырезается» область в виде цилиндра.
1.3 Процедура математической параметризации электронной пушки
Разработка математического описания электронных пушек включает в
себя два основных этапа. На первом этапе важной задачей является
определение геометрии электронной пушки в первом приближении, что
позволит определить основные параметры электронного луча, такие как
мощность, компрессия, первеанс, интенсивность электронного пучка. Как
12
правило, на этом этапе используется, так называемый, метод синтеза,
который заключается в упрощенной реализации математической модели.
Второй этап – метод анализа, на котором с помощью численного
эксперимента проводится решение самосогласованных задач по расчету
конфигурации электронного пучка и оптимизации геометрии электронной
пушки.
Одним из основных методов синтеза электронной пушки является
метод синтеза по Овчарову. Данный метод расчета заключается в получении
специального решения внутренней задачи, которое удовлетворяет некоторой
совокупности заранее требуемых характеристик формируемого пучка.
Построение решения осуществляется в параксиальном приближении с
помощью аналитической реализации параксиального уравнения, которое
записывается в специальной криволинейной системе координат. В качестве
координатной системы выбирается семейство линий, связанных между собой
критерием подобия (семейство а) и семейство линий ортогональных к ним
(семейство b). В системе координат q1q2 (рис. 1.6) линии, принадлежащие
этим семействам, описываются уравнениями:
r z
qR z ;
r
z
Rz
,
Rz
r
z
где R z –одна произвольная линия из базового семейства; q –параметр
подобия.
Такой выбор координатной системы позволяет численно реализовать
аналитическую модель потока на траекториях по критерию подобия.
Параметрическое моделирование ортогональной сетки выполняется по
следующему алгоритму. Для задания линий семейства а используется
параметр подобия
q2
q
r z R z ,который определяет относительное
расстояние данной линии семейства от оси симметрии. В качестве параметра
подобия, который лежит в основе построения линий семейства b,
принимается расстояние, на котором данная линия пересекает ось
13
симметрии, q1
z .Построенная на таких допущениях ортогональная сетка,
лежащая в основе координатной системы, предполагает, что параметры q1
и q 2 имеют значение криволинейных координат, причем, первая из них имеет
размерность длины, а вторая является безразмерной.
Рис. 1.6 Ортогональная координатная сетка, образованная семейством
подобных линий (линии q2=const) и семейством линий, ортогональных к ним
(линии q1=const)
В случае, если в качестве линии базового семейства R z
берется
граница электронного пучка, то параксиальное волновое уравнение, является
уравнением движения, и записывается как
R2
d 2U 0
dq12
2R
dR dU 0
dq1 dq1
4U 0 R
d 2R
dq12
1
0
e
2 U0
m
,
где U 0 –значение потенциала на оси электронного пучка; I – ток пучка.
Введем
нормированные
характеристики
R
R Rн , u U 0 U н и
q1 lн , получаем
q1
d 2u
R
dq12
dR du
2R
dq1 dq1
2
d 2R
4uR
dq12
1
,
u
(1.1)
где Rн , U н и l н – нормирующие множители, a
i
0
e
2 U н3 / 2
m
1
2
;
Rн
.
lн
(1.2)
14
Аналитическая модель электронной пушки строится по следующему
алгоритму. На первом этапе задается распределение потенциала на оси,
удовлетворяющее условию (1.1), которое учитывает ограничение тока с
катода пространственным зарядом, а также движение электронного пучка в
канале равного потенциала на выходе из пушки.
Картина,
соответствующая
описанному
выше
распределению
потенциала показана на рис. 1.7, б. В качестве начального приближения
распределение потенциала может задать аналитическим выражением вида
u
k1
f , где f
4/3
2
5
1
an q1 n .
n 1
Рис. 1.7. Форма электродов (а) и осевое распределение потенциала (б) в
пушке, рассчитанной по методу Овчарова.
На
d 2R
dq12
результирующее
решение
накладывается
требование
0 ,чтобы описать, что катод имеет сферическую форму. Это
требование выполняется, если коэффициенты
k , a1
и a 2 задаются по
следующему закону:
k
9
i
4
2/3
;
a1
8 dR
15 dq1 q
;
1 0
a2
361 dR
900 dq1
2
.
(1.3)
q1 0
Чтобы выполнить условия (1.3), необходимо наложить дополнительные
связи на коэффициенты a 3 , a 4 и a 5 :
15
119 1
9 k
a3
187 1
9 k
a4
3a2 ;
10 6a1
15 8a1 3a2 ;
a5
77 1
9 k
На
следующем
3a1
6.
a2
этапе
выполняется
численное
интегрирование
уравнения (1.1), при заданном значении параметра i и начальных условиях
R
q1 0
1и
dR
dq1
q1 0
dR
dq1
. В результате вышеописанных манипуляций
0
определяется конфигурация пучка электронов, которая соответствует
заданному распределению потенциала.
В криволинейно системе координат с ортогональной сеткой проводится
также и решение внешней математической задачи. Сначала требуется задать
начальное распределение потенциала снаружи пучка, это можно сделать,
например, по приближенной формуле:
U
где
2
u
U
d 2R
q 2 uR
dq12
U
Uн
–
определенным
2
i
1 q22
4 u
ln q22 ,
нормированное
шагом
U
const ,
значение
можно
потенциала.
рассчитать
Положив
с
распределение
эквипотенциальных линии в криволинейной системе координат q1 , q2 . Для
расчета цилиндрического пучка необходимо перейти к цилиндрическим
координатам с помощью следующей нормировки:
r
lн
1
q2 R 1
2
2 2
q2
dR
dq1
2
;
z
lн
q1
1
2
2 2
q2
dR
.
dq1
Математическая параметризация и оптимизация электронной пушки
заключается
в
последовательном
изменении
расчетной
геометрии
управляющей системы электродов катод-анод электронной пушки, а также
16
порождаемого магнитного поля до тех пор, пока конфигурация получаемого
электронного пучка не будет близкой к предполагаемым на практике. Этот
процесс включает состоит из следующих основных шагов:
- задание начальной геометрической формы пушки;
- начальное распределение потенциала и конфигурации магнитного
поля;
-расчет траектории формируемого пушкой электронного пучка;
- внесение поправок в исходную геометрию и траектории электронов
на новом шаге и т.д.
Построение математической модели электронной пушки заключается в
решении самосогласованной задачи электронной оптики, описываемой
системой уравнений:
- уравнения потенциала 2U
,
(1.4)
0
dV
- уравнение движения частиц m
dt
eE
- уравнение неразрывности потока div V
V
0.
B ,
(1.5)
(1.6)
Здесь U –величина потенциала, m , e – масса и заряд электрона
соответственно,
– плотность заряда, E – вектор напряженности
электрического поля, B – вектор индукции магнитного поля, V – скорость
электрона.
Совместное решение системы уравнений (1.4) – (1.6) выполняется
методом простой итерации. На первой итерации рассчитывается поле
потенциала электронной пушки, полагая величину пространственного
заряда
= 0. На следующих итерациях внешнее электрическое поле и
траектории движения электрона рассчитываются с учетом пространственного
заряда. Увеличение количества итераций продолжается до того момента,
пока значения, полученные на n-ой итерации не будут достаточно близки к
значениям последующей (n+1)-ой итерации. В качестве критерия сходимости
17
процесса рассматриваются относительные изменения координат
и
r
скоростей Vr на выходе из анодного узла электронной пушки:
rn
rn
1
тр ,
rn
где
тр
Vn Vn
Vn
1
тр
.
(1.7)
– заданная погрешность расчета траекторий.
Чтобы учесть распределение пространственного заряда в электронном
пучке используется, так называемая, дискретная модель потока из
деформируемых элементов. Для этого пучок электронов разбивается на слои
в поперечном сечении. При формировании каждого деформируемого
элемента учитывается условие равного тока на каждом слое, которое
получается из условия одинакового заряда на каждом элементе. При этом с
учетом неравномерного распределения плотности тока на катоде площади
слоев будут различными.
Заряд каждого деформируемого элемента рассчитывается по формуле:
q
I 0 t N сл ,
где N сл – число слоев; I 0 – ток луча; t – шаг интегрирования по времени.
Распределение по узлам сетки заряда на каждом деформируемом
элементе выполняется при помощи алгоритма «частица в ячейке», который
заключается в разбиении элемента на отдельные частицы.
Суммарный ток потока электронов складывается из общего вклада
токов с элементарных участков катода:
M
I0
Im ,
(1.8)
m 1
где М – число элементарных площадок на эмитирующей поверхности катода.
Из уравнения (1.3) по распределению потенциала вблизи поверхности
катода вычисляются токи с элементарных участков с использованием закона
степени трех вторых для плоского диода:
Im
2.33 10
3/ 2
6 U m Sm
2
d
,
(1.9)
18
где U m
– потенциал точки, лежащей по нормали напротив центра
элементарного участка на расстоянии d от катода,
Sm
– площадь
элементарного участка катода.
Чтобы увеличить скорость сходимости для итеративного процесса,
используется корректировка тока на катоде с использованием метода нижней
релаксации:
In
In
1
i
где I n , I n
1
In
In
1
,
(1.10)
– величина тока на n-ой и (n-1)-ой итерациях;
i<1
– весовой
коэффициент, называемый параметром нижней релаксации.
Электрическое поле, задаваемое уравнением (1.4), рассчитывается
методом конечных разностей с применением формулы верхней релаксации
на каждом приближении. В граничных узлах сетки, находящихся на
электродах пушки, задается условие Дирихле, а на открытых участках –
условие Неймана.
Уравнение движения (1.5) для расчета траектории электронов
решаются методом Рунге-Кутта 4-го порядка шаг интегрирования задается
автоматически. Начальные условия – координаты электронов на катодном
узле пушки и их начальные скорости (в первом приближении – тепловое
распределение) задаются при разбиении на слои.
Учитываемые
при
решении
уравнения
(1.5)
компоненты
напряженности электрического поля находятся из уравнения Лапласа (или
Пуассона) по значениям потенциала в ближайших узлах сетки методом
конечно-разностного дифференцирования. Компоненты вектора индукции
магнитного поля рассчитываются в параксиальном приближении по
формулам:
Bz r , z
где B z
Bz
2
r 2 B|| z , Br r , z
r |
B z ,
2
– распределение компоненты магнитной индукции в осевом
направлении электронной пушки.
19
2. Математические методы анализа электронных пушек
2.1 Расчеты электронных пушек методом синтеза.
Методы решения задач фокусировки электронных пучков можно
разделить на два класса:
1. В начальном приближении считаются известными геометрия и
управляющее напряжение на электродах, а также магнитное поле (если
требуется) системы формирования. Необходимо вычислить траектории
движения
электронов
с
учетом
или
без
учета
собственного
пространственного заряда потока электронов.
2. В начальном приближении считаются известными требуемые
траектории движения электронов. Необходимо вычислить геометрию и
управляющее напряжение внешних электродов, а также конфигурацию
магнитного поля (если требуется), обеспечивающие создание траекторий
движения.
Первый класс методов называется методом анализа, второй – методом
синтеза.
Одним из классических примеров использования метода синтеза
считается метод расчета геометрии электродов пушек Пирса, с учетом того,
что электроны движутся прямолинейно. Этот пример явно демонстрирует,
что метод синтеза подразделяется на две задачи – внешнюю и внутреннюю
задачи теории формирования. Это следует из того, что изначально задаются
траектории
движения
электронов,
затем
находится
распределение
потенциала внутри потока электронов (по методу Пирса – внутри катодной
системы), а затем рассчитывается (или последовательно подбирается) форма
фокусирующего электрода и анода снаружи пучка электронов, которые
обеспечивают требуемое распределение потенциала.
Недостатком решения задачи по методу Пирса является тот факт, что
отсутствует отверстие в аноде. Это вносит существенные погрешности в
20
расчет: приближение, на котором в рассмотрение задачи вводится учет
отверстия в аноде, картина резко меняется: вблизи анода существенно
меняются
траектории
электронов
(отклоняются
от
прямолинейного
движения) и распределение потенциала. Возникает так называемая анодная
линза, которая приводит к изменению распределения потенциала в поле
электронной пушки.
Поэтому большой практический и теоретический интерес представляет
задача
разработки
для
систем
формирования
электронных
пучков
последовательных методов синтеза, которые могут позволить рассчитывать
устройства, предсказывающие движение электронных пучков по заданным
траекториям электронов.
Основоположником метода синтеза считается Г. А. Гринберг. Метод
Гринберга состоит в следующем: уравнения движения заряженной частицы
записываются, в так называемой, натуральной системе координат, т.е. в
такой ортогональной системе, оси которой совпадают с направлениями
касательной, главной нормали и бинормали к траектории в каждой ее точке.
Такой метод задания координатной сетки позволяет решать и прямую и
обратную задачи формирования электронных пучков, т.е. или найти
траектории пучка по заданным электрическому и магнитному полям
внешних фокусирующих электродов, или найти внешние фокусирующие
поля по заданным траекториям.
Уравнения движения по методу Гринберга являются частным случаем
уравнений параксиальной электронной оптики. Они дают возможность
выполнить достаточно подробный теоретический анализ для систем
фокусировки с криволинейной осью и используется в ряде расчета
практических задач. По методу Гринберга в системе электронной пушки
анализируются только прямолинейные пучки заряженных части без учета его
собственного объемного заряда.
Крупный вклад в развитие метода синтеза сделал В.Т. Овчаров. Он
предложил для решения внешней и внутренней задач системы формирования
21
электронных пучков использовать криволинейную ортогональную систему
координат. Система выбирается таким образом, чтобы либо электронные
траектории лежали на одной из координатных поверхностей, либо чтобы
одна из ее координатных линий совпадала с заданными траекториями.
Использование
такой
криволинейной
системы
координат
позволяет
преобразовать задачу о нахождении управляющих напряжений внутри
электронного потока
к решению обыкновенного дифференциального
уравнения второго порядка. Такой подход позволяет рассматривать всю
систему формирования в виде единого целого.
На теорию синтеза для
систем формирования накладываются
определенные ограничения, которые связаны с введением необходимых
упрощающих
предположений.
Способы
решения
также
связаны
с
трудностями, как при расчетах, так и при решении внешней задачи, то есть с
геометрией электродов и конфигурацией магнитных полей.
Основным недостатком моделей электронных пушек, рассчитанных
методом синтеза является сложная геометрия вычисляемых фокусирующих
электродов, что в свою очередь связано с проблемами их технической
реализации.
Упростить
сложную
форму
фокусирующих
электродов
позволяет расчет электронных пушек методом анализа, который представлен
ниже.
2.2 Расчеты электронных пушек методом анализа.
При
расчете
электронных
пушек
методом анализа
в
первом
приближении считаются известными форма электродов, образующих
систему электродов катод-анод, распределение напряжения в них, а также
величина плотности объемного заряда в пространстве, которое ограничено
контурами
электродов.
Для
решения
задачи
по
распределению
напряженности электрического поля в электронной пушке, используются
22
различные методы, самым распространенным из которых является метод
конечных разностей.
Метод состоит в следующем: дифференциальное уравнение заменяется
конечно-разностным аналогом, который получается если провести замену
производных
конечно-разностными
Предположим,
что
рассчитываемое
приближенными
поле
потенциала
выражениями.
удовлетворяет
уравнению Пуассона в декартовой системе координат:
2
2
U
y2
U
z2
.
(2.1)
0
Пусть требуется найти в некоторой точке О рассматриваемой области
вторые производные потенциала. Из математического анализа известно, что
можно записать через значения первых производных в соседних точке О
точках а, b, с, d:
2
U
z2
1
h
U
a
z
U
с,
z
2
1
h
U
b
y
U
d.
z
U
y2
(2.2)
Первые производные в уравнении (2.2) можно выразить через
конечные разности:
U
a
z
1
U1 U 0 ,
h
U
c
z
1
U0 U3 ,
h
U
b
y
1
U2 U0 ,
h
U
d
y
1
U0 U4 ,
h
(2.3)
Здесь U1 , U 2 , U 3 , U 4 – величина потенциала в соседних точках 1, 2, 3,
4, окружающих точку О.
Подставив (2.3) в (2.2), получим:
2
U
z2
2
U
y2
1
U1 U 0
h2
1
U2 U0
h2
U0 U3 ,
(2.4)
U0 U4 ,
23
затем найдем сумму:
2
2
U
z2
U
y2
1
U1 U 2
h2
U3
U4
(2.5)
4U 0 .
Таким образом, получим для уравнения Пуассона его конечноразностный аналог:
U1 U 2
U3
U4
h2
4U 0
.
(2.6)
0
Аналогично,
можно
получить
конечно-разностный
аналог
для
уравнения Лапласа в декартовой системе координат
U1 U 2
U3 U4
4U 0
0.
(2.7)
Соответственно для цилиндрической системы координат конечноразностный аналог для уравнения Пуассона можно записать, как:
2
U
r2
1
r
U1 U 2
U
r
2
U
z2
U3 U4
;
(2.8)
0
h
U2 U4
2r0
4U 0
h2
(2.9)
0.
где r0 – расстояние от рассматриваемой точки до оси симметрии.
Если точка, лежит на оси симметрии, то вместо (2.9) для нее
записывается уравнение:
U1 U 3
4U 2
Дискретный
6U 0
аналог
h2
(2.10)
0.
связывает
значения
потенциала
в
каждой
отдельной дискретной точке. Таким образом, чтобы вычислить поле в
расчетной области, она разбивается на квадратную сетку с шагом h.
Дискретный аналог записывается в каждом узле, который лежит внутри
области рассмотрения. Конечно-разностное уравнение связывает потенциал в
каждом узле с четырьмя прилежащими данному узлу соседними узлами
сетки. Для узлов, лежащих на границе расчетной области, задаются
начальные значения потенциала, например, по контуру электродов.
24
Конечно-разностные уравнения (2.6) – (2.10) записываются в узлах
сетки, и образуют систему линейных алгебраических уравнений, число
которых равно числу неизвестных. Т.е., в такой постановке можно свести, к
решению системы алгебраических уравнений решение краевой задачи.
Система дополняется граничными условиями, которые задают значения
потенциала в опорных точках и граничных узлах.
Чтобы минимизировать погрешность, связанную с заменой конечноразностным уравнением дифференциального, следует уменьшить шаг сетки,
т.е., соответственно увеличить количество узлов и, как следствие, увеличить
порядок системы уравнений. Количество узлов в расчетах иногда может
превосходить
десять
тысяч,
что
приводит
к
невозможности
непосредственного решения системы уравнений методом исключений.
Поэтому для решения используется метод простых итераций, или иначе
последовательных приближений. На данный момент этот метод, имеющий
ряд разновидностей, получил широкое применение при расчетах полей на
ЭВМ.
Для расчета траекторий движения электронов в электронных пушках,
широко используется метод последовательных приближений, который
заключается
в
следующем.
В
первом
приближении
используется
конфигурация электрического поля без учета собственных полей пучка
электронов. На втором приближении рассчитываются поля и траектории с
приближенным учетом собственных полей потока электронов. Подбор
пошагового процесса последовательных приближений применяется до того
момента, пока значения функций текущего п-го приближения не будут
достаточно
близкими
к
значениям
функции
предыдущего
(n – l)-гo
приближения. Чтобы оценить сходимость процесса используется критерий,
на котором можно анализировать сходимость, например, координат и углов
расходимости траекторий движения электронов в некоторой выбранной
плоскости анализируемой системы. Если процесс сходится, то для получения
25
конечного результата с заданной точностью обычно необходимо 5 – 10
приближений.
Для решения самосогласованных систем методом последовательных
приближений применяется дискретная модель пучка электронов в виде
траекторий – токовых трубок. На входе в анализируемую систему пучка
электронов разделяется в поперечном направлении на N элементарных слоев
– токовых трубок. Тогда в каждой трубке
Ik парциальный ток
рассчитывается по вычисленной площади поперечного сечения трубки и
распределению
плотности
тока
в
поперечном
сечении
пучка
(в
предположении, что последнее считается известным). Этот ток объясняется
одной «центральной» траекторией трубки, ход которой рассчитывается в
будущем. В этом случае решение системы самосогласованных уравнений
приводится
к совместному решению уравнений движения, поля и
непрерывности тока. Последний, примененный к данной модели пучка,
имеет вид
Ik = const. Известное распределение заряда вычисляет поле
следующего приближения и т.д.
2.3 Способы расчета реальных магнитных полей в электронных пушках
Наиболее популярными в применении на данный момент являются, так
называемые датчики э.д.с. Холла, т.е. полупроводниковые измерители
магнитных полей. Датчики э.д.с. Холла позволяют измерять как переменные,
так и переменные магнитные поля.
Если через боковые грани образца проходит постоянный ток, то в
присуствии внешнего магнитного поля возникнет поперечная разность
потенциалов, которая называется эффектом Холла. Если образец из
полупроводникового материала сделан в форме параллелепипеда, то эта
разность потенциалов задается формулой
26
Uy
R
ix H z
10 8 в
d
(2.11)
где Нz – напряженность магнитного поля, ix – сила тока в образце, R –
константа Холла, d – толщина образца.
Поэтому можно по формуле (2.11) оценить, что если пропускается
постоянный ток через образец, то в образце будет возникать напряжение
электрического поля, которое будет пропорциональным напряженности
магнитного поля с точностью до десятых долей процентов.
В качестве примера рассмотрим задачу, которая позволяет рассчитать
форму
и
геометрию
электронного
пучка
от
катода
до
анода,
проанализировать возможность сходимости пучка, а затем на основе
рассчитанных параметров, оптимизировать электронную пушку.
Исходные данные:
1. Напряжение на аноде – 52 кВ;
2. Число пучков электронов – 40;
3. Конфигурация пучков электронов:
а) 21 луч с диаметром 84 мкм,
б) 19 луч с диаметром 64 мкм;
4. Канал системы катод-анод 6,5 – 8 мм;
5. Суммарный первеанс
20
10-6 А/В3/2;
6. Диаметр катода – 8,6 мм.
Для расчета электронной пушки на основе метода синтеза будем
использовать теорию Овчарова. Внутренние траектории рассчитываются из
крайней по формуле
r
q2
0
где
Z
,
l
(2.12)
– функция, задающая крайнее значение траектории потока электронов в
цилиндрической системе координат; Ф0 – радиальные размеры пучка; l –
единица измерения продольных размеров пучка; q2 – криволинейная
27
ортогональная координата; Z – продольная координата цилиндрической
системы координат; r – радиальная координата цилиндрической системы
координаты.
В криволинейной системе координат крайняя траектория потока имеет
q2 = 1, осевая q2 = 0, а остальные 0< q2 <1.
Чтобы
получить
аксиального
рассчитанную
электронного
луча,
конфигурацию
необходимо
симметричного
решить
следующее
дифференциальное уравнение:
4
2
u
2
Здесь
u
4u
2
h2
4
k
hk2
2
i
.
u
(2.13)
x – функция, задающая крайнее значение траектории потока
электронов и совпадающая c выражением (2.12); и(x) – задает осевое
распределение потенциала; h(x) – задает осевое распределение магнитного
поля; h = h(0) – значение функции h(x) на катоде;
функции
=
0 – значение
x на катоде.
Т.к. в осевом направлении пучка криволинейная система координат
является цилиндрической, следовательно, функции и(х) и h(x) совпадают с
функциями, которые описывают, соответственно, осевое распределение
потенциала и магнитного поля в цилиндрической системе координат.
Дифференцирование в уравнении (2.13) проводится по продольной
координате х. Постоянная в (2.13) задается выражением
i 0,297
H0 l
.
V0
(2.14)
здесь Н0 – напряженность магнитного поля; l – продольные размеры пучка;
V0 – потенциал.
Постоянная
i
в
(2.12)
определяет
ток
потока
электронов.
Микропервеанс пучка определяется по i и по значению потенциала V0
формулой:
28
0,0605 P
i
где
2
(2.15)
,
= (Ф0 / l); P – микропервеанс пучка.
Параксиальная
теория
позволяет
решать
внешнюю
задачу
в
криволинейной системе координат используя трансцендентное уравнение
4
2
V
2
u
q22
(u
h2
4
k
hk2
2
4
)
(2.16)
2
i
1 q22 ln q22 ,
4 u
где V = U /U0 – потенциал по заданному значению эквипотенциальной линии.
Уравнение (2.16) разрешается относительно функции q2(x) для каждого
значения x.
По результатам анализа находится значение функции q2*(x), по
которому определяется форма искомой эквипотенциальной линии в
криволинейной ортогональной системе координат.
Затем требуется вернуться от криволинейной системы координат к
цилиндрической с помощью параметризации
2
d
dq 2
1
q2
2
q2 ,
(2.17)
которая выполняется с начальными условиями:
q2
0;
Интегрирование выполняется до q 2
x.
q 2* , где q2* - решение уравнения
(2.17) для текущего значения x.
Соответствующее q2* значение переменной
есть
*, которая
используется дня вычисления цилиндрических координат r и z:
r
l
Z
l
q 2*
;
(2.18)
*
29
На практике уравнения (2.17) и (2.18), описывающие внешнюю задачу,
можно решить только численно при помощи электронных вычислительных
машин.
Начальное приближение распределения функции потенциала внутри
пучка в первом приближении параксиальной теории формировании в
криволинейной системе координат задается формулой
2
где
V1
u
V1
–
2
q 22
(u
потенциал
4
h2
2
4
заданной
4
k
hk2
),
(2.19)
эквипотенциальной
линии.
Функция
распределения плотности тока внутри пучка электронов в криволинейной
системе координат задается однородной.
2.4 Расчет электростатических электронных пушек
В качестве примера рассмотрим задачу, со следующими единицами
измерения:
- Ф0 – радиус в начале пучка, который определяет радиальные размеры
системы,
- l – расстояние от катода до оси пучка, который определяет продольные
размеры системы,
- U0 – значение потенциала.
Для решения внутренней задачи по расчету электростатических
электронных пушек используются два варианта: или известно задание
траекторий
электронного
рассчитывается
осевое
пучка,
а
по
распределение
нему
из
потенциала,
уравнения
или,
(2.13)
известна
конфигурация осевого распределения потенциала, а из уравнения (2.13)
рассчитываются траектории движения электронов.
Для распределения потенциала, которое задается функцией и(х), а
также для траектории электронов, которая задается функцией (x) в расчете
30
электронной пушки применяются определенные условия, называемые
граничными условиями для функции и(х):
0, u x
0, и x
x 1, u 1, и
u
x
Условия
(2.20)
0;
(2.21)
0.
(2.20)
определяют
работу
катода
при
наличии
пространственного заряда, а условия (2.21) – определяют нулевой потенциал
на оси электронной пушки.
Для функции (x) при х=0 используется условие:
x
(2.22)
0.
Условие (2.22), определяет сферичность эмитирующей поверхности
катода.
Будем рассчитывать электронную пушку задавая функцию и(х), и
вычисляя функцию (x). Тогда функция и(х) задается таким образом, чтобы
выполнялись условия (2.20), (2.21), но также необходимо еще задавать
функцию и(х) в области малых значений х, при которой функция
(x),
определяемая из выражения (2.12), соответствовала бы выполнению условия
(2.22).
Если найден такой способ задания функции и(х), то, после проведения
расчета нескольких вариантов внутренней задачи, можно построить схему
рекомендаций как проводить расчет электронных пушек, и формировать на
выходе пучок электронов с определенными параметрами.
Чтобы решить уравнение (2.22) требуется задание начальных условий.
Расчетная модель внутренней задачи электронной пушки выстраиваться в
направлении от катода, с заданием граничных условий и(х) и и'(х) при х = 0.
Но в этом случае, когда на катоде и(х) = 0, то правая часть уравнения (2.22)
стремится к бесконечности. Эту проблему можно разрешить следующим
образом.
Задается
функция
и(х),
затем
находится
приближенное
аналитическое решение уравнения (2.12), которое справедливо в области
малых х. Это решение принимается за начальное приближение, и уже с этого
31
аналитического решения выполняется первая итерация с катода в точку, где
значение функции и(х) уже не равно нулю. Далее проводится решение
уравнения (2.12) методом простых итераций до достижения сходимости.
Выберем в качестве начального аналитического решения вид функции и(х)
следующей формулой:
u
kx4 / 3 f 2 ,
f
1
(2.23)
где
5
an x n ,
(2.24)
n 1
k, an (n = 1,2, …..5) – некоторые постоянные коэффициенты.
Немаловажен тот факт, что функция и(х), которая задается выражением
(2.23), остается всегда положительной (при положительных значениях х) и
удовлетворяет граничному условию (2.20). А при малых значениях х
функция и(х) совпадает с функцией kx4/3, описывающей осевое распределение
потенциала.
Коэффициенты а1, а2, многочлена (2.23) выбираются так, чтобы
выполнялось условие (2.22), а коэффициенты а3, а4, a5 должны удовлетворять
условию (2.21).
Чтобы найти соответствующие коэффициенты а1, а2, определим для
функции и(х), которая задается выражением (2.13), приближенное решение
уравнения (2.12), являющееся справедливым при малых х.
В такой постановке задачи решение для функции
(х) ищется в виде
полинома
5
x
1
bn x n ,
(2.25)
n 1
Из этого многочлена вытекает, что значение
"(x) при х = 0 будет
определяться значением b2. Таким образом, чтобы выполнялось условие
(2.22) нужно получить такие значения для коэффициентов an, при которых b2
будет равно нулю. Для этого подставим формулы (2.23), (2.25) в уравнение
(2.12) и, затем приравняем к нулю коэффициенты при одинаковых степенях
32
х, и выразим bn через an. После полученных манипуляций получим, что
коэффициенты а1, а2, учитывая также, что должно выполняться условие b2=
0, должны соответствовать следующим выражениям:
a1
8
b1 ;
15
(2.26)
a2
361 2
b1 .
900
(2.27)
Из уравнения (2.25) следует, что коэффициент b1 однозначно
определяет значение первой производной от функции (x) в точке x = 0, т.е.
на катоде. Поэтому введем обозначение b1
k , с учетом которого
формулы (2.26) и (2.27) запишутся:
a1
8
15
k;
a2
361
900
k
2
(2.28)
(2.29)
.
Данный расчет позволяет также установить, что при малых х
коэффициенты к, b3, b4 связаны с постоянными коэффициентами i, а3, а4
выражениями:
k
9 2/3
i ;
4
33 74377
37 222750
b3
b3
(2.30)
0,228771
3
k
4
k
a3 ;
1,154518
(2.31)
k
b3
(2.32)
0,783582 a4
Эти соотношения позволяют приближенно вычислить решение
уравнения (2.12), которое справедливо при малых х, если значения
коэффициентов а3, а4 известны.
Имея
начальное
приближение,
необходимо
найти
значения
коэффициентов а3, а4, а5, которые удовлетворяют условиям (2.21). Вычислим
первую и вторую производные от функции и(х) и в точке х = 1 наложим
33
следующие условия: u 1
0. В результате будет получена
0, и 1
1, u 1
система трех уравнений, решив которую относительно а3, а4, а5, получим:
a3
119
9
9
i
4
a4
187
9
a5
77
9
9
i
4
9
i
4
1/ 3
48
15
10
1/ 3
10
1/ 3
24
15
64
15
k
361
300
k
k
361
900
361
900
k
2
k
;
(2. 33)
2
k
2
;
(2.34)
6.
(2.35)
Уравнения (2.23), (2.28), (2.29), (2.33) – (2.35) лежат в основе
определения способа задания функции и(х), для которого выполняются как
граничные условия (2.20), (2.21), накладываемое на функцию и(х), так и
граничное условие (2.22), накладываемое на функцию (х).
Как только найдена функция и(х), находится решение внутренней
задачи для электростатической электронной пушки, т.е. решение уравнения
(2.12).
Найдем решение уравнения (2.12) с помощью конечно-разностного
аналога при следующих начальных условиях: x
x0 ,
0,
0.
Величину параметра х0 выбирается достаточно малым (0,0001 + 0,01), а
значения
и
0
’0 для точки х = х0 вычисляется в соответствии с (2.25) по
следующим выражениям:
0
0
1
k
k
x0
x0
b3 x03
3 b3 x02
b4 x04 ;
(2.36)
4 b4 x03 .
(2.37)
Значения коэффициентов b3, b4 при малых значениях х, вычисляются
по выражениям (2.31), (2.32), а входящие в них значения а3, а4, а5,
определяются формулами (2.32) – (2.35).
Решение
уравнения
(2.12)
конечно-разностным
методом
проводиться до критической точки xкр, в которой производная
будет
’(х) равна
нулю, т.е. до узловой точки пучка.
34
Для решения внутренней задачи моделирования электронной пушки
требуется задать значения параметров i,
. Параметр i, как вытекает из
уравнения (2.14), является характеристикой первеанса электронной пушки.
Параметр
задает радиус кривизны катода пушки (Rкp), вычисляемый из
выражения:
Rкр
1
l
;
k
Можно кончено-разностным методом решить также внешнюю задачу.
Для этого из уравнения (2.15) сначала находим решение внешней задачи в
криволинейной системе координат, и используя его как начальное
приближение, находим решение уравнения (2.16), осуществляя таким
образом переход к цилиндрической системе координат. Решая внешнюю
задачу требуется задание параметра V = U/U0, где U – значение потенциала
того электрода, геометрическую форму которого требуется вычислить. При
анализе геометрической формы прикатодного фокусирующего электрода
параметр V считается равным нулю, а если требуется вычислить форму анода
пушки, то значение параметра V рассчитывается по формуле
2
V
1
4
i
(1 ln b 2 ),
(2.38)
где b = rn/r – коэффициент заполнения канала пучком; rn, r – соответственно
радиусы пучка и пролетного канала.
Из уравнения (2.29) вытекает характеристика провисания потенциала
по отношению к трубе дрейфа прибора, которая заполняется потоком
электронов с микропервеансом Р и коэффициентом заполнения b. Это
вытекает из уравнений (2.12), (2.15) при условии (2.21).
После того, как найдены решения для внутренней и внешней задач по
описанному выше алгоритму, требуется при помощи (2.28) найти радиус
кривизны катода Rкp. Этот параметр характеризует площадь эмитирующей
поверхности и находится из точки пересечения дуги радиуса Rк с графиком
функции (х).
35
Обобщив полученные значения внутренней задачи для электронной
пушки,
можно
составить
методику
расчета
пушки
с
ожидаемыми
характеристиками.
Выполнение граничного условия (2.12) сводит функцию
(х) для
малых значений х к прямой линии, которая образует с осью х угол
. Это
позволит вычислить радиус катода r по результатам решения задачи по
пересечению данной прямой с дугой окружности, радиус которой
вычисляется из выражения (2.28). Выполнив все расчеты, будем иметь:
rк
l
1
2
2
1
k
(2.39)
;
Пусть радиус кривизны пучка обозначается как rкр. Тогда rкр – это
значение функции
(х), получаемое из функции кривизны
p,
который
вычисляется из формулы:
rкр
(2.40)
кр ,
l
Определим линейную сходимость потока электронов, как отношение
радиуса катода rк к радиусу кривизны пучка rкр. Из выражений (2.31), (2.40)
можно определить линейную сходимость электронного потока следующим
выражением:
1
S
кр
1
2
(2.41)
;
k
Характеристика зависимости
p
от
, получается как результат
реализации внутренней задачи в зависимости от набора значений параметра
, которые лежат в интервале 1,2
2,4. В данном расчете величина
значения параметра i остается постоянным и равно 0,4. Полученную
зависимость аппроксимируем формулой
кр
1,05 0,709
k
0,125
2
k
.
(2.42)
Подставляя (2.42) в (2.41), получаем:
36
b
S
1,05 0,709
k
0,125
2
k
1
2
(2.43)
;
k
Выражение (2.43) может использоваться для расчета интервала
значений параметра
, который позволяет построить модель пушки со
сфокусированным пучком при заданных значениях сходимости S.
При создании методики расчета электростатической электронной
пушки полагаем, что заданы первеанс пучка электронов Р , линейная
сходимость электронного пучка S и коэффициент заполнения канала
пучком b.
Чтобы
решить
внутреннюю
приближение параметров i,
задачу
нужно
задать
начальное
, а чтобы решить внешнюю задачу –
дополнительно указать значения коэффициента
и потенциалов V1, V2. По
значению потенциала V1 определяется форма прикатодного фокусирующего
электрода пушки, а значению потенциала V2 форма анода пушки. Поэтому
величина V1 полагается равной нулю, а величина V2 вычисляется по
начальным значениям Р и b из выражения (2.31). Величина параметра i
выбирается равной 0,4. Величина параметра
находится из начального
приближения S и рассчитанному значению
из выражения (2.34). Это
уравнение трансцендентное и его можно решить только численно.
Когда результирующие величины параметров i, V1,
рассчитаны,
проводится
окончательной
расчет
,
параметров
численно
пушки,
результатом которого является сформированный потому электронов с
определенными параметрами.
Описанный алгоритм расчета реализован в программе «Расчет
параметров электронной пушки». Эта программа вычисляет геометрию
электронной пушки, а затем конфигурацию пучка электронов. Для
вычисления необходимо задать шесть параметров:
H – ширина области моделирования,
R – радиус управляющего электрода,
37
– угол раствора управляющего электрода,
r – глубина погружения катода,
L – расстояние между катодом и анодом,
D – диаметр выходного отверстия анода.
В результате расчета определяется теоретическая и технологическая
геометрия
электронной
пушки
и
по
результатам
рассчитывается
коэффициент пропускания.
38
3. Проведение модельных экспериментов по исследованию параметров
пушки
Для
проведения
модельных
экспериментов
был
разработан
и
использован программный модуль «Расчет параметров электронной пушки».
Исследования проводились для определения оптимальных параметров
настройки пушки.
Анализ литературных источников показал, что при определении
оптимальных
параметров
оптимальную
величину
настройки
зазоров
пушки
между
необходимо
элементарными
выбрать
ленточными
катодами. Если выбранная величина зазоров будет слишком велика, то при
высоких температурах вследствие линейного теплового расширения кернов с
размещенным на них эмиссионным слоем, образующим эмитирующую
поверхность каждого катода, произойдет увеличение длины катодов на
величину, меньшую, чем величина разделяющих их зазоров. При этом между
элементарными ленточными катодами остаются не покрытые эмиссионным
слоем промежутки, что препятствует получению общей эмитирующей
поверхности, близкой к непрерывной поверхности, и снижает равномерность
распределения электронного потока по сечению. Если выбранная величина
зазоров будет слишком мала, то при высоких температурах элементарные
ленточные катоды могут увеличить свою длину на величину, превышающую
величину разделяющих их зазоров, вследствие чего в местах смыкания
элементарных
катодов
могут
возникнуть
механические
напряжения,
приводящие к деформации и даже к разрушению катодов. При этом может
произойти изменение положения эмитирующей поверхности элементарных
катодов относительно других электродов электронной пушки, то есть
искажение геометрии ЭОС пушки. Это также может привести к снижению
равномерности распределения электронного потока по сечению, а также
уменьшению надежности и долговечности катодного узла и электронной
пушки.
39
Таким образом, для обеспечения надежной работы пушки и увеличения
продолжительности ее работоспособности нужно:
1.
Создать однородный тепловой режима катодного узла
2.
Получить более однородную общую эмитирующую поверхность
3.
Сохранить ее заданное геометрическое расположение в ЭОС
электронной пушки
Анализ библиографических источников показал, что величина зазора
между элементарными ленточными катодами зависит от коэффициента
теплового расширения материала катодного цилиндра при максимальной
температуре элементарного ленточного катода, длины катодного цилиндра
при комнатной температуре и максимальной температуры элементарного
ленточного катода
Минимальная величина зазора определяется величиной максимально
возможного удлинения катодного цилиндра, выполненного из материала с
заданным коэффициентом теплового расширения, при максимальной
температуре,
до
которой
катод
может
быть
нагрет
в
процессе
технологического цикла изготовления, испытания и работы электронной
пушки. Если данное условие выполняется, то удается при достижении
катодом максимальной температуры избежать деформации элементарных
ленточных катодов и выхода их из строя и сохранить геометрию ЭОС
электронной пушки, что обеспечит равномерность распределения по сечению
протяженного электронного потока и приведет к повышению надежности и
долговечности катодного узла и электронной пушки.
Выбор
элементарными
максимально
ленточными
возможной
катодами
величины
зазоров
определяется
между
назначением
электронной пушки, в которой будет использован данный катодный узел.
При работе электронной пушки в результате термического расширения
происходит увеличение длины каждого катодного цилиндра элементарного
40
ленточного катода, в результате эти сближаются друг с другом, образуя
общую эмитирующую поверхность, близкую к непрерывной поверхности.
Имеет место зависимость размеров эмитирующей поверхности катода
от выходной мощности в импульсе, подаваемого на анод пушки
ускоряющего напряжения и тока на выходе.
Так, например, в работе [] эксперименты на макетах электронных
пушек показали, что на расстоянии 3-4 см от выходного окна электронной
пушки, выходящий из нее протяженный электронный поток заданной
мощности будет однородным по своей длине на уровне 90% при условии, что
зазоры между элементарными ленточными катодами не превысят 2,3 мм.
Для формирования формы пучка на практике достаточно выдержать
форму электродов вблизи границы пучка. Если заданы напряжение и ток
пучка, а также его поперечный размер на выходе из пушки, то тогда расчет
пушки сводится к определению расстояния анод - катод.
Площадь катода легко определить по заданным размерам пучка, что
позволяет оценить плотность тока на катоде и найти расстояние анод - катод.
Наличие
отверстия
в
аноде
пушки
приводит
к
образованию
рассеивающей линзы-диафрагмы, фокусное расстояние которой находится
как
отношение
потенциала
анода
с
коэффициентом
к
разности
напряженностей слева и справа от анода. И рассматриваемые пушки будут
давать на выходе расходящиеся пучки, у которых углы расхождения
пропорциональны отношению радиуса анода к расстоянию анод – катод. Из
этого следует, что расфокусирующее действие анодного отверстия возрастает
по мере увеличения отношения толщины потока к междуэлектродному
расстоянию. При практическом моделировании пушки была использована
упрощенная
форма
электродов,
которая
обеспечивает
необходимое
распределение потенциала лишь вблизи границы пучка.
На основании всего вышесказанного было решено исследовать
зависимость
моделирования,
пропускающей
радиуса
способности
управляющего
от:
ширины
электрода,
угла
области
раствора
41
управляющего электрода, глубины погружения катода, расстояния между
катодом и анодом, диаметра выходного отверстия анода на предмет
выявления тех значений указанных параметров, при которых пропускающая
способность достигает максимального значения.
При последовательном подборе параметров исследуется влияние
каждого отдельно взятого фактора на пропускающую способность при
некоторых фиксированных значениях других факторов. По мере получения
оптимальных значений ими заменяются соответствующие фиксированные.
Сначала был определен диапазон угла раствора управляющего
электрода.
Первоначально угол был взят в диапазоне от 0 до 45 градусов.
Проведенные эксперименты (рисунки ) показали, что диапазон можно сузить,
так как именно при углах от 21 до 24 градусов коэффициент пропускания
принимает максимальные значения.
Рисунок 3.1. Первоначальный подбор угла раствора
42
Рисунок 3.2. Изменение зависимости коэффициента пропускной способности
от угла при изменении расстояния между катодом и анодом
Рис. 3.3 Изменение зависимости коэффициента пропускной способности от
угла при изменении напряженности
Из графиков зависимостей коэффициента пропускной способности от
угла при изменении расстояния между катодом и анодом при различных
значениях последнего (рис. 3.2) можно увидеть, что наилучшее расстояние
должно находиться в пределах [60; 80].
Зафиксировав значение расстояния между катодом и анодом на уровне
70мм, исследовали влияние напряженности (рис. 3.3). В результате получили
слабую степень влияния изменения коэффициента пропускной способности
43
от угла при значениях напряженности от 5 до 20кэВ. Наилучшим значением
напряженности было выбрано равное 15кэВ.
Далее исследования проводились при уменьшении диапазона углов
сначала от 18 до 27 градусов, а затем от 21 до 24 (рис. 3.4 – 3.6) для
определения наилучшего значения диаметра выходного отверстия анода D.
Рис. 3.4 Значения коэффициента пропускной способности от угла
(18<=α<=27) при изменении диаметра выходного отверстия анода D
44
Рис. 3.5 Значения коэффициента пропускной способности от угла
(21<=α<=24) при изменении диаметра выходного отверстия анода D
Рис. 3.6 Изменение зависимости коэффициента пропускной способности от
угла при различных значениях диаметра выходного отверстия анода D
45
Следующим шагом явилось определение оптимального значения
глубины погружения катода. Результаты проведенных экспериментов
представлены на рисунке 3.7.
Рис. 3.7 Изменение зависимости коэффициента пропускной способности от
угла при различных значениях глубины погружения катода
При изучении влияния радиуса управляющего электрода на изменение
коэффициента пропускной способности было обнаружено, что наилучшим
диапазоном значений данного показателя является [70;80] (рисунок 3.8).
46
Рис. 3.8 Изменение зависимости коэффициента пропускной способности от
угла при различных значениях радиуса управляющего электрода
Предварительное предположение о том, что при любых значениях
ширины области моделирования, большей порогового значения, дальнейшее
ее увеличение не оказывает влияния на результаты моделирования,
экспериментально подтвердились (рис. 3.9).
Рис. 3.8 Изменение зависимости коэффициента пропускной способности от
угла при различных значениях ширины области моделирования.
47
Вся проведенная работа позволила выявить наилучшие диапазоны
изменения влияющих на конфигурацию пушки параметров.
Следующим этапом стал подбор конкретных точечных значений
параметров из рассмотренных диапазонов их оптимальных значений.
В результате были получены следующие наилучшие значения
влияющих параметров (рис.3.9):
радиус управляющего электрода R = 75 мм,
угол раствора управляющего электрода
= 22.5 градуса,
глубина погружения катода r = –2 мм,
расстояние между катодом и анодом L = 80 мм,
диаметр выходного отверстия анода D = 1.5 мм.
а)
б)
Рис. 3.9 Изменение зависимости коэффициента пропускной способности от:
а) Напряжения на катоде при фиксированных значениях;
б) диаметра отверстия на аноде при различных значениях угла
По результатам проведенного моделирования был создан практический
образец.
Экспериментальные
результаты
полностью
соответствуют
рассчитанным параметрам.
48
Основные результаты и выводы.
1. Проведено математическое моделирование методом последовательных
приближений с нижней релаксацией;
2. Модель является устойчивой для большого количества управляющих
внешних параметров и демонстрирует хорошую сходимость;
3. Разработан
комплекс
программ
и
проведен
вычислительный
эксперимент, который позволил установить оптимальные параметры
электронной пушки;
4. По результатам проведенного моделирования был создан практический
образец. Экспериментальные результаты полностью соответствуют
рассчитанным параметрам.
49
Заключение
В работе было проведено математическое моделирование электронной
пушки с управлением конфигурацией пучка. На основе изучения и анализа
библиографических
аналитическое
источников
представление
по
задачи
теме
ВКР
было
моделирования
рассмотрено
формы
пучка
электронов, определены недостатки существующих ЭП, рассмотрены
вопросы применения методов анализа и синтеза для решения задачи. Были
исследованы оптимальная геометрия катодного узла; влияние внешних
параметров модели на конфигурацию электронного пучка; влияние на
фокусировку пучка; движение и нагрев дисперсной частицы в потоке
плазмы. На основе проведенных исследований было выявлено, что для
определения конфигурации электронной пушки необходимо провести
эксперименты по выявлению зависимостей пропускающей способности от:
ширины области моделирования,
радиуса управляющего электрода,
угла раствора управляющего электрода,
глубины погружения катода, расстояния между катодом и
анодом,
диаметра выходного отверстия анода
на предмет выявления тех значений указанных параметров, при
которых пропускающая способность достигает максимального значения.
Была разработана компьютерная программа и проведены расчеты
оптимальных параметров пушки. Была построена модель ЭЛП и проведено
тестирование разработанной пушки.
В дальнейшем планируется изучение траектории электронов в
различных модификациях электронной пушки с учетом особенностей
внешней среды.
50
Список использованных источников
1 ―Математические
электронных
модели
приборов‖.
и
автоматизированное
Методические
указания
проектирование
к
лабораторным
работам. Под редакцией В.К. Федяева. №2143. Рязань 1993. 60 с.
2 С.И. Молоковский А.Д. Сушков. ―Интенсивные электронные и ионные
пучки‖. М.: Энергоатомиздат, 1991. – 304 с.
3 А.А. Жигарев Г.Г. Шамаева. ―Электронно-лучевые и фотоэлектронные
приборы‖. М.: Высшая школа, 1982. – 463с.
4 Климаков
А.А.,
Моделирование
электронной
пушки
с
системой
фокусирующих линз с учетом распределения пространственного заряда //
Вестник СПБГУ. серия 10: Прикладная математика. Информатика.
Процессы
управления,
Издательство:
Санкт-Петербургский
го-
сударственный университет (Санкт-Петербург), 2013, С. 225—229.
5 Рошаль А.С. Моделирование заряженных пучков. М.: Атомиздат, 1979.
6 К.Г. Симонов. Электронные отпаянные пушки. - М.: Радио и связь, 1985. Гл.3.
7 А.Б. Киселев. Металлооксидные катоды электронных приборов. -М.: изд.
МФТИ, Физматкнига. 2002. - Гл.2.
8 В. Эспе. Технология электровакуумных материалов. - Том I. - M.:
Госэнергоиздат, 1962. - С.146 и 159.
9 Диденко А.Н., Лигачев А.Е., Куракин И.Б. Воздействие пучков
заряженных частиц на поверхность материалов и сплавов. – М.:
Атомиздат, 1987. – 187 с.
10 Диденко А.Н., Григорьев В.П., Усов Ю.П. Мощные электронные пучки и
их применение. – М.: Атомиздат, 1977. – 210 с.
11 Кулешов А.А., Земсков С В., Позняк Ю.В. Электронное учебное пособие
по высшей математике на базе системы Mathematica
12 Шмидский Я.К. Mathematica 5. Самоучитель. – М.: Диалектика, 2004.
13 Матросов А.В. Maple 6. Решение задач математики и механики. – СПб.:
БХВ-Петербург, 2001.
51
14 Рыжиков Ю.И. Решение научно-технических задач на персональном
компьютере. – СПб.: Корона принт, 2000.
15 Эдвардс Ч.Г., Пенни Д.Э. Дифференциальные уравнения и краевые
задачи: моделирование и вычисление с помощью Mathematica, Maple и
MATLAB. – 3-е издание. – Киев: Диалектика-Вильямс, 2007.
16 Дьяконов В.П. Mathematica 4.1/4.2/5.0 в математических и научнотехнических расчетах. – М.: СОЛОН-Пресс, 2004.
17 Макаров Е.Г. Самоучитель Mathcad 14, 2012. [Электронный ресурс]
Режим доступа: http://mirknig.com/2012/02/08/samouchitel-Mathсad-14.html
18 Очков В.Ф. Mathcad 14 для студентов и инженеров: русская версия. СПб.:
BHV, 2012 г. - 368 с.
19 Пискунов В. В Работа в Mathсad [Электронный ресурс] Режим доступа:
elib.ispu.ru/library/lessons/pekunov/index.html
20 Яньков В. Ю., Якушина Н.А. Решение прикладных задач в пакете
"Маткад". – М.: Спутник, 2011г. – 155с.
52
Отзывы:
Авторизуйтесь, чтобы оставить отзыв