Научная статья на тему 'О решении двумерного уравнения Пуассона, когда облако объемного заряда имеет гладкие криволинейные границы. Ii. Аксиальная геометрия'

О решении двумерного уравнения Пуассона, когда облако объемного заряда имеет гладкие криволинейные границы. Ii. Аксиальная геометрия Текст научной статьи по специальности «Математика»

CC BY
135
16
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Научное приборостроение
ВАК
RSCI
Область наук

Аннотация научной статьи по математике, автор научной работы — Шевченко С. И.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

On the Solution of Two-Dimensional Poisson's Equation for the Case of a Space-Charge Cloud with Smooth Curvilinear Boundaries. II. AXIAL GEOMETRY

The paper presents a solution algorithm to the axial-geometry two-dimensional Poisson's equation by the boundary integral equation method for the case when the space charge cloud has smooth curved boundaries.

Текст научной работы на тему «О решении двумерного уравнения Пуассона, когда облако объемного заряда имеет гладкие криволинейные границы. Ii. Аксиальная геометрия»

ISSN 0868-5886

НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2007, том 17, № 1, c. 73-82

ИССЛЕДОВАНИЯ, ПРИБОРЫ, МОДЕЛИ-И МЕТОДЫ АНАЛИЗА

УДК 517.925.4: 62 - 408.64

© С. И. Шевченко

О РЕШЕНИИ ДВУМЕРНОГО УРАВНЕНИЯ ПУАССОНА, КОГДА ОБЛАКО ОБЪЕМНОГО ЗАРЯДА ИМЕЕТ ГЛАДКИЕ КРИВОЛИНЕЙНЫЕ ГРАНИЦЫ. II. АКСИАЛЬНАЯ ГЕОМЕТРИЯ

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

ВВЕДЕНИЕ

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

Как и в указанных работах, для вычисления потенциала электростатического поля (ЭСП) будем пользоваться интегральным уравнением, эквивалентным уравнению Пуассона [3],

0(ro) =\aGdL +-!- f QGdS

(L) ™0 (S)

(1)

где £0 = 8.8542 -10 12 Ф/м — электрическая постоянная; Q — функция плотности объемного заряда (ПОЗ); <г(г) — функция плотности поверхностного заряда (ППЗ); О(г, г0) — ядро интегрального уравнения (ядро для задачи Дирихле для уравнения Лапласа); точка г0, в которой ищется потенциал, — точка наблюдения (ТН); первый интеграл берется по длине контура всех электродов (границы), а второй — по площади, которую занимает объемный заряд (ОЗ); коэффициент 1/4яе0, который должен быть перед первым членом в правой части (1), включен в функцию о.

В отличие от случая плоской геометрии в случае аксиальной геометрии ядро для задачи Дирихле для уравнения Лапласа имеет вид [4, 5]

G(r,Го) =

4r

<J(r + Го)2 + (z - zo)

(2)

где К(£)— полный эллиптический интеграл первого рода,

k2 =■

4rr0

(r + ro)2 + (z - Zo)2

В [6] после некоторых преобразований ядра (2) было получено выражение, в котором удалось выделить логарифмическую сингулярность:

О (г, Гз) = Л (г, г) + /2( г, г )1п[(г - Г0)2 + (7 - ^)2], (3)

где

fl( z, r) = -

4r

гХ

7(r + Го)2 + (z - Zo)2

x{g (z,r) + g2(z, r)ln[(r + ro)2 + (z - zo)2]}; 4r

f2( zr)=-1 2 2;

yl(r + ro)2 + (z - zo)2

4

&i( zr) = E aM ;

i=0 4

g2( ^r)=X bm;

m1( z, r) = 1 - m = 1 - k2 =

(r - ro)2 + (z - zo)2 (r + ro)2 + (z - zo)2

значения коэффициентов а{ и bi приведены в справочнике [7]. Входящие в (3) функции /1, /2 являются гладкими везде в области г Ф 0 .

Подробности решения уравнения (1) с ядром (3) можно найти в [2]. Отличие второго интеграла в правой части (1), представляющего собой вклад в потенциал от всего облака ОЗ, в аксиальном случае от аналогичного интеграла в плоском случае определяется видом ядра (3). Этот интеграл, в котором Q и О считаются известными, можно находить численным методом (например, повторным интегрированием методом Гаусса), методом аналитического интегрирования или аналитически-численным методом.

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

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

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

1. ОСНОВЫ АЛГОРИТМА

Для вычисления вклада в потенциал от всего облака пространственного заряда необходимо вычислить второй интеграл в правой части (1). Этот интеграл с помощью введения прямоугольной сетки разбивается на сумму интегралов по площади тех ячеек, которые содержат ОЗ. К таким ячейкам относятся "полные" прямоугольные ячейки, вклад от которых рассмотрен в [2], и "граничные" ячейки, т. е. ячейки, которые пересекаются границей облака ОЗ и которые поэтому находятся частично вне, а частично внутри облака ОЗ, вклад в потенциал от которых будет рассмотрен в данной работе.

Для случая плоской геометрии вклад в потенциал от некоторой граничной ячейки (т. е. от ОЗ, содержащегося в этой граничной ячейке) рассмотрен в [8], где были выделены следующие проблемы:

1) анализ,

2) позиционирование,

3) интерполяция,

4) интегрирование по площади граничной ячейки.

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

1.1. Анализ

Вопросы разбиения всех граничных ячеек на разные типы, отнесения каждой граничной ячейки к одному из установленных типов и относительного положения ТН и рассматриваемой граничной ячейки подробно рассмотрены в [9]. Показано, что всего возможны 12 вариантов пересечения ячеек сетки границей расчетной области. Выписан алгоритм анализа полученных вариантов пересечений ячеек сетки границей расчетной области, позволяющий для каждой ячейки установить реализуемый вариант пересечения. В данной работе в качестве границы расчетной области используется граница облака ОЗ.

1.2. Позиционирование

При рассмотрении второго интеграла в правой части (1) важно знать, где расположена ТН (как она позиционируется в пределах граничной ячейки). Возможны следующие варианты [9]: ТН лежит вне, внутри, точно на границе ячейки сетки.

1.3. Интерполяция

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

В отличие от рассмотренного в [1, 8] случая плоской геометрии для аксиальной геометрии ядро (3) является более сложным. Дополнительно появляются функции f1 и f2, для которых требуется осуществить некоторую интерполяцию. В [2], где для аксиального случая рассматривались только вклады в потенциал от ячеек, полностью заполненных объемным зарядом, и использовалась интерполяция в пределах каждой ячейки произведений fQ и f2Q билинейным полиномом, были получены не вполне удовлетворительные результаты.

В данной работе применялось раздельное интерполирование функций f1 и f2 , ядра (3) и функции Q(х, у) билинейным полиномом

fl,2 ■ Q = (4,2 + Bl,2х + CuУ + D12ху)X

x (A + Bx + Cy + Dxy) (4)

для граничной ячейки четырехугольной формы. Для треугольной формы граничной ячейки или той (треугольной) части пятиугольной граничной ячейки, которая получается после разбиения пятиугольной граничной ячейки на четырехугольную ячейку и треугольную, — D = 0, D12 = 0 .

После простейших преобразований получаем для произведений fl2 ■ Q :

f1,2 ■ Q = 2 xmy". (5)

m,n=0

Коэффициенты связаны с коэффициентами

A, B, C, D, A12, B12, C12, D12 соотношениями, приведенными в Приложении I.

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

Q(x, у) = со,о + ci,o x + С01 y + с1Д xy + c2 0x2 +

+ C21X2y + C0 2 y2 + Ci,2 xy2 + C2 2x2y2.

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

1.4. Интегрирование по площади граничной ячейки

Так же, как в работе [8], посвященной плоскому случаю, для нахождения вклада в потенциал от граничных ячеек будем использовать три метода вычисления второго интеграла в правой части (i), определяющего вклад в потенциал от облака объемного заряда:

1) аналитическое вычисление интегралов;

2) численное вычисление;

3) аналитически-численное вычисление интегралов.

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

В случае аксиальной геометрии вклад в потенциал от всего облака пространственного заряда имеет вид (второй интеграл в правой части (i))

Vp =

1

4П£,

J dx dy Q(x,y) • G(x,y;xo,Уо) . (6)

0 (S)

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

т. е. сумму вкладов в потенциал фр(¿х,\ ) как от

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

Вклад в потенциал от одной граничной ячейки имеет вид

Фр (ix , i У ) =

1

4пе0

-X

хк Ук(х)

х| 4х | 4у • х,у) • О(х,у;х0,у0), (7)

XN уИ (х)

где пределы "внутреннего" интегрирования: yN(х)— нижний предел и ук (х)— верхний предел и пределы "внешнего" интегрирования: xN — нижний предел и хк — верхний предел определяются конфигурацией граничной ячейки, по площади которой проводится интегрирование. Один из пределов "внутреннего" интегрирования (по переменной у) является зависящим от "внешней" переменной (х).

При любом из методов нахождения вклада в потенциал можно указать два пути повышения точности решения уравнения Пуассона:

• повышение точности интерполяции функции плотности ОЗ;

• повышение точности интерполяции границы облака ОЗ.

1.4.1. Аналитическое интегрирование

Если в выражение (7) подставить ядро (3) и применить интерполяцию (5), то для вклада в потенциал от одной граничной ячейки получим:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

ф (i , i ) = У ( • / О) + С(2) . / (2) ) ; (8)

i pv x' y' л / j у m ,n m,n m ,n m,n J ' v '

4n^o m n = 0

(9)

o m,n=0

Ук(x)

С = J xmdx J yndy;

xN yN (x)

=

Лк ук(x)

J xmdx J yndy • ln[(x - x0)2 + (y - У0)2]. (10)

yN (x)

Более сложный вид ядра определяет то, что вклад в потенциал от некоторой граничной ячейки в случае аксиальной геометрии содержит (в сравнении с вкладом в потенциал от этой же граничной ячейки в случае плоской геометрии) дополнительные члены, содержащие интегралы от гладких функций, и дополнительные члены, содержащие интегралы 1,022), 11(2), 1:222), 12220), 12? с логарифмической сингулярностью и другими степенями переменных интегрирования в подынтегральной функции.

Выражения для интегралов 1{^п являются простыми, а выражения для тех интегралов , которые не использовались в [8], приведены в Приложении II для случая, когда переменным является верхний предел "внутреннего" интегрирования

(у к = ук(х) = «• х +в ).

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

При повышении точности интерполяции границы облака ОЗ (т. е. переходе от кусочно-линейного представления границы к более точному представлению (кусочно-квадратичному или сплайн-представлению)) "внутреннее" интегрирование в интегралах (8) проводится без изменения (по сравнению с кусочно-линейным представлением границы). Во "внешнем" интегрировании возникают интегралы, являющиеся практически не-интегрируемыми. Все попытки преодолеть эту трудность приводили к такому увеличению времени счета, которое делало бессмысленным применение аналитического интегрирования.

1.4.2. Численное интегрирование методом Гаусса

Подробности вычисления вклада в потенциал от плоской граничной ячейки, заполненной объемным зарядом, методом повторного численного интегрирование по формуле Гаусса можно найти в [8].

При вычислениях этим методом, вклад в потенциал от ОЗ, содержащегося в треугольной или четырехугольной граничной ячейке, определяется выражением, реализующим в (7) повторное численное интегрирование методом Гаусса [10]:

Ф = —

Р 4п£о

Л х х

хУ А

/ XI

1=1

Лу (X )У Аут X, Ут ) °(X, Ут I X0, Уо)

, (11)

где Лх = Рх [,х + 1] - F_SX [/у ] для треугольной граничной ячейки или Лх = Рх[/х +1] - Рх[,х] для четырехугольной граничной ячейки; Рх [/х ] — линия сетки, параллельная оси х и имеющая номер ,х; F_SX [,у ] — х-координата точки пересечения линии сетки р [,у ], параллельной оси у, и линии границы облака пространственного заряда; Л у (х1) = = ук (х,) - Ру^у ]; , И^у — число узлов интегрирования вдоль направлений х и у ; Ах1, Аут — коэф-

фициенты квадратурной формулы Гаусса; хг = Р [,х ] + Лх • х^ (,) — координаты узлов интегрирования вдоль направления х; ут = Ру [,у ] + +Л () • хN (т) — координаты узлов интегрирования вдоль направления у; х^ (,) и х^ (т) — корни каждого узла интегрирования; Q(х1, у.) — значения плотности объемного заряда в узлах интегрирования (х1, yj).

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

Для аксиального случая формула двукратного численного интегрирования (11) имеет практически такой же вид, как такая же формула для плоского случая. Отличие заключается лишь в том, что в аксиальном случае в эту формулу следует подставлять значение ядра, вычисляемое по формуле (3).

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

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

1.4.3. Смешанное (аналитически-численное) интегрирование

Если взглянуть на формулу (10), то видно, что "внутреннее" интегрирование (по переменной у) не содержит никаких особенностей и сводится к вычислению табличных интегралов [11]

ук(х)

Л(х) = / у^у • 1п [х2 + у2 ] ;

уN (х)

формула записана в системе координат, привязанной к ТН.

Вторая фаза интегрирования ("внешнее" интегрирование) осуществляется численным интегрированием по методу Гаусса.

Положительным качеством смешанного интегрирования (или преимуществом по сравнению с аналитическим интегрированием) является то,

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

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

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

2. ТЕСТОВЫЙ ПРИМЕР И ОБСУЖДЕНИЕ

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

Выберем треугольную граничную ячейку, которая относится к типу 9 по классификации работы [9]. Для равномерной сетки 10^10 соотношение времен расчетов трех обсуждаемых алгоритмов нахождения вклада в потенциал (аналитическое интегрирование, аналитически-численное интегрирование и двукратное численное интегрирование) составляет примерно 1 : 2.1 : 2.5. Соотношение времен выполнения расчетов вклада в потенциал от одной граничной ячейки не должно зависеть от разбиения сетки. Только общее время вычисления вклада в потенциал от всех граничных

1) В предыдущей работе [8] данного цикла была допущена опечатка: вместо написанного в ней "сферический конденсатор" следует читать "цилиндрический конденсатор"

ячеек должно зависеть от разбиения, т. к. от этого разбиения зависит общее число граничных ячеек.

В [8] для плоского случая приведено сравнение точности вычисления интегралов. В аксиальном случае часть интегралов являются такими же,

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

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

Результаты расчета распределения потенциала и вклада в потенциал от всего облака пространственного заряда вдоль радиуса приведены в табл. 1 и табл. 2. При расчете потенциала использован метод аналитического интегрирования, а расчет вклада в потенциал от всего облака пространственного заряда выполнен тремя методами: А — аналитический, Б — численный, АБ — аналитически-численный. Разбиение сетки — 10*10 для данных табл. 1 и 100*100 в табл. 2. Вычисления проделаны для равномерного распределения объемного заряда, расположенного между обкладками сферического конденсатора.

Сравнение результатов, представленных в табл. 1, 2, показывает, что при переходе к более частой сетке ошибки в вычислении потенциала и вклада в потенциал от всего облака пространственного заряда значительно уменьшаются, что можно объяснить повышением точности представления криволинейной границы. Уже при разбиении 10*10 потенциал получается с вполне удовлетворительной точностью.

ЗАКЛЮЧЕНИЕ

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

Табл. 1. Распределения потенциала и вклада в потенциал от всего облака пространственного заряда (колонки А, Б, АБ) вдоль радиуса. Разбиение сетки 10*10.

№ — номер точки вдоль радиуса; г — радиальное расстояние точки, рчисл. — значения потенциала, полученные в результате работы алгоритма; ртеор. — теоретические значения потенциала; дР = (рчисл. - ртеор.) / ртеор.тах — величина относительной ошибки; А, Б, АБ — методы вычислений: соответственно аналитический, численный, аналитически-численный

№ г Р 1 числ. Р 1 теор. дР А Б АБ

0 0.070 0.00012 0.00000 0.00012 0.0895844 0.0898250 0.0898250

1 0.073 0.13919 0.13963 -0.00044 0.0900915 0.0903257 0.0903211

2 0.076 0.26714 0.26773 -0.00059 0.0901326 0.0903484 0.0903517

3 0.079 0.38490 0.38559 -0.00069 0.0897511 0.0899660 0.0899603

4 0.082 0.49358 0.49431 -0.00074 0.0889794 0.0891765 0.0891792

5 0.085 0.59412 0.59485 -0.00073 0.0878415 0.0880325 0.0880323

6 0.088 0.68734 0.68803 -0.00068 0.0863549 0.0865337 0.0865371

7 0.091 0.77395 0.77455 -0.00060 0.0845341 0.0847147 0.0847082

8 0.094 0.85456 0.85503 -0.00047 0.0823913 0.0825530 0.0825576

9 0.097 0.92970 0.93002 -0.00032 0.0799350 0.0800950 0.0800939

10 0.100 0.99973 1.00000 -0.00027 0.0771680 0.0773198 0.0773198

Табл. 2. Распределения потенциала и вклада в потенциал от всего облака пространственного заряда вдоль радиуса. Разбиение сетки 100*100. Обозначения те же, что в табл. 1

№ г Р числ. Р теор. дР А Б АБ

0 0.070 0.00012 0.00000 0.00012 0.0960336 0.0960336 0.0960336

1 0.073 0.13960 0.13963 -0.00002 0.0966684 0.0966684 0.0966684

2 0.076 0.26768 0.26773 -0.00004 0.0968104 0.0968104 0.0968104

3 0.079 0.38553 0.38559 -0.00005 0.0964926 0.0964926 0.0964926

4 0.082 0.49425 0.49431 -0.00006 0.0957435 0.0957434 0.0957434

5 0.085 0.59479 0.59485 -0.00006 0.0945876 0.0945876 0.0945876

6 0.088 0.68797 0.68803 -0.00006 0.0930464 0.0930464 0.0930464

7 0.091 0.77450 0.77455 -0.00005 0.0911384 0.0911384 0.0911384

8 0.094 0.85499 0.85503 -0.00004 0.0888799 0.0888798 0.0888799

9 0.097 0.93000 0.93002 -0.00002 0.0862851 0.0862851 0.0862850

10 0.100 0.99989 1.00000 -0.00011 0.0833666 0.0833665 0.0833665

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

Проведено сравнение аналитического, численного и смешанного (аналитически-численного) методов интегрирования при нахождении вклада в потенциал.

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

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

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

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

Численный метод (двукратного численного интегрирования по формуле Гаусса) является самым медленным из трех рассмотренных. И по точности он несколько уступает другим рассмотренным методам. Однако в рамках этого метода возможно не только использование более точной интерполяции для функции плотности ОЗ и для границы облака ОЗ (чем в аналитическом методе), но и применение для них аналитических функций.

Приложение I. Соотношения между коэффициентами стп и А, В, С, Б, А 2, #12, С 2, Л 2

С(1'2) = А и0,0 1,2

с°'2) = С

0,1 1,2 с(1'2) = В

4,1 1,2 с(1.2) = 2,1 с(1,2) =

= В,,

= В

А ; с&2) = #1,2 • А + Аи • В ; .А + А 2 • ;

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

С + С12 • В + А12 • Б + Д,2 • А ; Б + А 2 • В; с(122) = С12 • Б + А 2 • С ;

В

(1,2)

(1,2)

= А, • Б .

Приложение II. Вид интегралов 1. 2.; I, j = 0 ^ 2

1022) = (а+1)

2 1п(а2 +1) х5к -х^

+ (хк XN )

2

2а _ 2 9 9

а, 4 4 ч

+ "7"(хк — XN ) 6

^ УN-—а2 в + + 1п(а +1) в2

а 2 3 1п(а +1) „

---а + —--ав

6 9 • 4 4

2

( хК х2 )

ав2

+ ёх

2 у3 _ 2 в3 _ 2^ + - у5 9 ^ 9 в 3 + 3 ^

_ 1п( х2к + у2)

3 3

+ У^хъ

Л- с»- ~ Л- с»-

3 к 3 к

+ 1п( х2 + у2)

3 3

Укх + У^хъ

3 % 3 х

N

( хк хы ^

аг^^—^ _ arctg

V

УN

у N

1 5 _ 2 4

15 У N 3 У N

1 (2) = 2.0

а(1п(а4+1) _ 2) (хк _ х'„ )+(хК _ хх; ))((а2 + 1) _ 2) + 2у„) + 6Щ

+ ^ [ х'к 1п(х'к + У2) _ 41п(4 + у1) ] +

+ ёх

2у N У3

в

2,0 3ах1 + в 2 / „. 3 , о 2>, ---;-х2 + 2-

2 2(а +1) 2 3

-х22 _ (_ах3 + вх12)

(4 _ ^ +(4 _ 4)

_6ах1 + 2 в в

2(а +1)

ах

(4 _ 4 )-

+ 1п( 2к + х22)

_ 1п( Z2N + х22)

а, 4 4. 3ах1 + в 3 3ах1 2в / 2 ^^ , 3 о 2\

-(zK_X22) + ■

^к +■

-(ZкK + х22) + (_ах13 + вх12) zк

а

-(zN_X22) +

_3ах1 + в з -Зах1_2в

" ZN +-

(z— + х2) + (_ ахС + вх12) ZN

Г хк хN ^

aгctg—^ _aгctg "

У N

УN

"2 у 4+^ уn

^ х! ах, + в х4 ах,г + в ^ — аг^ —-—- —- aгctg —^^—-х,к '^2 ^х

N

(

arctg— - arctg-

ßx:

-3axi i ß

2(а i i)

2x231 (-ax3 i ßxi2)2x2

^ arctg ^ - ^ arctg ^

2 Xk 2 XN

Л

J(2) = (x5 - x5 )

2,i УЛК aN>

а2 i ln(a i i), 2 14

--1 — i —-- (а i i)

iG 25 iG

i (XK xn)

-aß+aßln(a21 i)

44

i ( XK XN )

б б i5 9

ln(a21 i)

i dx

4 o 5

У4 „4 2У5

"Уы"

- ln(X2 i yN)

G.i • xl i ^xl 3 K

i ln(xN i y2N)

G ix5 i v-lxN T 3 XN

75

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

xK xN arctg—1- - arctg N

V

yN

yN

2 X6 — X6

J2,2 _ -9a 4 ^ (XK - )

2 20 ln(a21 ib 2

—УЫ -— a ßi

2 i5

i5

iG

(a i i)

i(XK xn) i ( - XN )

i dx

a aß2 i ln(a21 i) б б 4 ß

9 27 Уы ( Xk Xn ) 27 27

2ß3 i A i ß2 ln(a21 i)

2yN a21 i 4 a21 i 3 2) 2 - 2 2

^X 2 i ( ^X 2 i Xi ) ^X 2 2 Xi ^X 2

- ln(x2k i y2) i & xN ln(xN i y2) i &2yNf

X к X N

arctg —- - arctg

Л

УЫ

a21 i 25

( Z5 Zn) Xi

a21i, 4 4.

-(zl - z4) i

a i i .2/ 3 3

i5

"X2( Z3 ZN )

a21 i

Уы i52

( x22 x23 1 X2 Л

yJ

13

а2 ii

Xi X2 ( ZK ZN ) i

i ln(zK i x22)

a21 i

iG

'ZK i Xi

a21 i . 4 4. a21 i

——(zK -x¿)

/ 3 2\ 3 а i i 2/2 2\ 2 ( X 2 i Xi ) Zк Xi X 2 ( Z к i X 2 ) i Xi X 2 Z к

- ln( zN i x^)

(

a21 i

iG

'ZN i Xi

а 11^4 44 a i К 3 2\3 а i i 2/2 2 \ 2

" ( zn - X2 ) i " ( X2 i Xi ) ZN Z Xi X2 ( ZN i X2 ) i Xi X2 ZN

Z Z

arctg —K - arctg—

ai i, 3 9ч 3 W TI í

( X 2 i Xi ) X 2 i X 2

a21 i

i '^2 Xi X 2

Y z z Л

arctg— - arctg N

V

a

a

У

11,2 , (хк хы) + 2. (хк хы) + в (хк хы) + „ (хк хм) + „ (хк х2) +

(

1п(а2 +1) 2

15

„3 хк

9 3

«-+ ^[а2в(хк _ х4) + ■3-ав2(х3 _ х3) + в(х2)

-У- [(хк + у2 )1п( хк + у2) _ (х2 + у2)1п( х2 + у2) ]+

9

+ а4^4 _zAN) + 3а2(в_аx1)(z3 _z3) + 3а(в_ах1)[ах1 +1](zKK _z—) + + (в _ ах1 )2 (в _ 4ах1) • ёх_х1(в_ ах1 )3,

где ёх = хк _ х,

х1 х2

х2 / (а2 +1)'

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

1. Иванов В.Я., Шевченко С.И. О расчете плоских электростатических полей в приборах, имеющих области, заполненные объемным зарядом // Научное приборостроение. 1999. Т. 9, № 4. С. 88-94.

2. Шевченко С.И. О расчете аксиально-симметричных электростатических полей в областях, заполненных объемным зарядом // Научное приборостроение. 2002. Т. 12, № 2. С. 2329.

3. Соболев С.Л. Уравнения математической физики. М., 1966. 400 ^

4. Ильин В.П. Численные методы решения задач электрофизики. М.: Наука, 1985. 334 с.

5. Иванов В. Я. Методы автоматизированного проектирования приборов электроники. Новосибирск: Изд-во СО АН СССР, 1986. 193 с.

6. Шевченко С.И. Алгоритм получения высокой точности в расчетах аксиально-симметричных электростатических полей // Научное приборостроение. 2002. Т. 12, № 1. С. 40-45.

7. Абрамовиц М., Стиган И. Справочник по специальным функциям. М.: Наука, 1979. 832 с.

8. Шевченко С.И. О решении двумерного уравнения Пуассона, когда облако объемного заряда имеет гладкие криволинейные границы. I. Плоская геометрия // Научное приборостроение. 2006. Т. 16, № 1. С. 85-93.

9. Шевченко С.И. О расчете траекторий заряженных частиц вблизи криволинейных поверхностей // Научное приборостроение. 2005. Т. 15, № 1. С. 70-79.

10. Крылов В.И., Шульгина Л.Т. Справочная книга по численному интегрированию. М.: Наука, 1976. 370 с.

11. Двайт Г.Б. Таблица интегралов. Наука, 1973. 228 с.

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

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

ON THE SOLUTION OF TWO-DIMENSIONAL POISSON'S EQUATION FOR THE CASE OF A SPACE-CHARGE CLOUD WITH SMOOTH CURVILINEAR BOUNDARIES.

II. AXIAL GEOMETRY

S. I. Shevchenko

Institute for Analytical Instrumentation RAS, Saint-Petersburg

The paper presents a solution algorithm to the axial-geometry two-dimensional Poisson's equation by the boundary integral equation method for the case when the space charge cloud has smooth curved boundaries.

i Надоели баннеры? Вы всегда можете отключить рекламу.