Рис. 3.2. Форма электродов (а) и осевое распределение потенциала (б) в пушке, рассчитанной по методу Овчарова.
При таком задании потенциала условия (3.2) выполняются автоматически. На искомое решение налагается требование
|
Реализация условий (3.3) обеспечивается за счет связей, наложенных на коэффициенты a3, a4и a5:
Далее производится численное интегрирование уравнения (3.1), для чего необходимо выбрать значение параметра i и задать начальные значения
Решение внешней задачи также производится в криволинейной системе координат. Расчет потенциала вне пучка может быть произведен по приближенной формуле:
где
Типичная форма электродов рассчитанной таким способом пушки показана на лис 3.2а.
Анализ пушки состоит в последовательном изменении геометрии электродов пушки и формы магнитного поля до тех пор, пока параметры формируемой пушкой пучка не будут близки к заданным. Этот процесс включает в себя следующие основные этапы: выбор исходного варианта геометрии пушки и конфигурации магнитного поля, траекторный анализ, по результатам которого определяются параметры формируемого пушкой пучка, внесение изменений в исходную геометрию и последующий траекторный анализ нового варианта т.д.
Анализ электронной пушки основывается на решении самосогласованной задачи электронной оптики, математическая модель которой включает:
- уравнения поля
- уравнение движения частиц
- уравнение неразрывности потока
Здесь U – потенциал, r - плотность заряда, m, e – масса и заряд электрона, Е – напряженность электрического поля, В – индукция магнитного поля, V – скорость частицы.
Совместное решение уравнений (1) – (3) выполняется методом последовательных приближений. В первом приближении производится расчет поля электронной пушки без учета пространственного заряда. На втором и последующих приближениях внешнее поле и траектории рассчитываются с учетом пространственного заряда. Процесс последовательных приближений продолжается до тех пор, пока результаты последующего n-го приближения не будут достаточно близки к результатам последующего ( n+1 )-го приближения. В качестве критерия сходимости процесса могут служить относительные изменения радиальных координат r и скоростей Vr контрольных электронов в конце расчетной области пушки:
где
Распределение пространственного заряда в пучке учитывается с помощью дискретной модели потока из деформируемых элементов. Электронный поток разбивается в поперечном сечении на слои. Формирование элементов проводится из условия получения одинакового заряда каждого элемента, и следовательно, одинакового тока каждого слоя. Площади слоев с учетом неравномерного распределения плотности тока на катоде будут различными.
Заряд деформируемого элемента находится по формуле:
где I0 – ток луча; Nсл – число слоев; Dt – шаг интегрирования.
Разнесение по узлам сетки заряда деформируемого элемента производится с помощью разбиения его на отдельные частицы и применения к ним алгоритма “частица в ячейке”.
Ток пучка определяется суммированием токов с элементарных участков катода:
где М – число дроблений эмитирующей поверхности катода.
Токи с элементарных участков вычисляются по найденному из уравнения ( 1 ) распределению потенциала вблизи поверхности катода с использованием закона “степени 3/2” для плоского диода:
где Um – потенциал точки, лежащей по нормали напротив центра элементарного участка на расстоянии d от катода,
Sm – площадь элементарного участка катода.
Для ускорения сходимости последовательных приближений применяется корректировка катодного тока по способу нижней релаксации:
где In, In-1 – токи на n-ом и (n-1)-ом приближениях; wi<1 – параметр нижней релаксации.
Расчет электрических полей, описываемых уравнением (1), проводится методом конечных разностей с использованием итерационной формулы последовательной верхней релаксации. В граничных узлах сетки, находящихся на электродах пушки, задается условие Дирихле, а на открытых участках – условие Неймана.
Уравнение движения (2) для контрольных электронов решаются методом Рунге-Кутта 4-го порядка с автоматическим выбором шага интегрирования. Начальные условия – стартовые координаты и скорости контрольных электронов на катоде – определяются при формировании слоев.
Необходимые для решения уравнения (2) составляющие напряженности электрического поля находятся по потенциалам ближайших узлов сетки с применением численного дифференцирования. Составляющие индукции магнитного поля рассчитываются в параксиальном приближении по формулам:
где B(z) – известное распределение продольной составляющей магнитной индукции на оси пушки.
4. Предварительный расчет
ВАРИАНТ №18 | ||
Pвых,кВт | lср, см | Df/fср,% |
20,5 | 11,2 | 1,4 |
a. Определим fср :
b. Зададим максимальный КПД :
если 5см < lср < 15 см, то 50% < hmax < 70%. Выберем hmax = 60%.
c. Определим рабочий КПД :
d. Определим мощность источника питания :
e. Определим значение микропервианса, обеспечивающего требуемую полосу пропускания :
f. Рассчитаем ускоряющее напряжение :
g. Определим ток луча :
h. По выбранному значению gaопределим радиус пролетной трубы а, см.
Выберем gaиз интервала 0,4…1,0. Меньшие значения относятся к длинноволновым клистронам, большие – к клистронам сантиметрового диапазона : ga= 0,9.
Для определения радиуса пролетной трубы, рассчитаем величину g по формуле:
Найдем a :
i. Определим радиус луча :
Для этого выберем коэффициент заполнения Кз из интервала 0,5…0,8.
Кз = 0,8.
Рассчитаем радиус луча :
j. Рассчитаем плотность тока луча :
k. Выберем предварительное значение коэффициента сходимости Ак, который является тангенсом угла наклона огибающей пучка на катоде, из интервала –1,8… –2,4 :