Главная Учебники - Разные Лекции (разные) - часть 55
Министерство общего и профессионального образования РФ Воронежский государственный университет факультет ПММ
кафедра Дифференциальных уравнении
“Моделирование распределения потенциала в МДП-структуре” Исполнитель : студент 4 курса 5 группы Никулин Л.А. Руководитель : старший преподаватель Рыжков А.В. Воронеж 1998г. ОГЛАВЛЕНИЕ
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ РАСПРЕДЕЛЕНИЕ ПОТЕНЦИАЛА В МДП-СТРУКТУРЕ
Математическая модель - - - - - - - - - - - - - - - - - - - 3
ПРИМЕНЕНИЕ МАТЕМАТИЧЕСКИХ МЕТОДОВ К
РЕШЕНИЮ ЗАДАЧИ
Использование разностных схем для решения
уравнения Пуассона и для граничных условий
раздела сред
Уравнение Пуассона - - - - - - - - - - - - - - - - - - - - - - 5
Граничные условия раздела сред - - - - - - - - - - - - - - - - - - - - - - 8
Общий алгоритм численого решения задачи
Метод установления - - - - - - - - - - - - - - - - - - - - - - 10
Метод переменных направлений - - - - - - - - - - - - - - - - - - - - - - 13
Построение разностных схем - - - - - - - - - - - - - - - - - - - - - - 16
ПРИЛОЖЕНИЕ - - - - - - - - - - - - - - - - - - -
- - - - - - - - - - - - - - - - - -
Математическая модель распределения потенциала в МДП-структуре Математическая модель
Пустьj
(
x,y
)
- функция, описывающая распределение потенциала в полупроводниковой структуре. В области оксла (С
DEF
) она удовлетворяет уравнению Лапласа: d2
j
+
d2
j
= 0
dx2
dy2
а в области полупроводника (прямоугольник ABGH
) - уравнению Пуассона: d2
j
+
d2
j
=
0
dx2
dy2
где q
- элементарный заряд e
; e
nn
-диэлектрическая проницаемость кремния; N
d
(x,y)
-распределение концентрации донорской примеси в подложке ; N
a
(
x,y
)
-распределение концентрации акцепторной примеси в подложке; e
0
-диэлектрическая постоянная 0
D
E
y
B G
C F
A H
x
На контактах прибора задано условие Дирихле: j
|
BC
= Uu
j
|
DE
= U
з
j
|
FG
= Uc
j
|
AH
= Un
На боковых сторонах полупроводниковой структуры требуется выполнение однородного условия Неймана вытекающее из симметричности структуры относительно линий лежащих на отрезках AB
иGH
: d
j
=
0
d
j
= 0
dy
AB
dy
GH
На боковых сторонах окисла так же задается однородное условие Неймана означающее что в направлении оси OY отсутствует течение электрического тока: d
j
=
0
d
j
= 0
dy
DC
dy
EF
На границе раздела структуры окисел- полупроводник ставится условие сопряжения : j
| -0
=
j
| +0
e
ok
E
x
|-0
-
e
nn
E
x
|+0
= - Qss
где Qss
-плотность поверхностного заряда; e
ok
-диэлектрическая проницаемость окисла кремния; e
nn
-диэлектрическая проницаемость полупроводника . Под символом “+0
” и”-
0
” понимают что значение функции беретсябесконечно близко к границе CF
со стороны либо полупроводникалибо окисла кремния . Здесь первое условие означает непрерывностьпотенциала при переходе границы раздела сред а второе - указывает соотношение связывающее величину разрыва вектора напряженностипри переходе из одной среды в другую с величиной поверхностногозаряда на границе раздела. ПРИМЕНЕНИЕ МАТЕМАТИЧЕСКИХ МЕТОДОВ К
РЕШЕНИЮ ЗАДАЧИ
Использование разностных схем для решения уравнения Пуассона и для граничных условий раздела сред
Уравнение Пуассона
В области{(x,y) : 0 < x < Lx
, 0 < y < Ly
}
вводится сетка W={(x,y) : 0 < i < M1
, 0 < j < M2
}
x0
=0 , y0
=0,
x
M
1
= Lx
, y
M
2
= Ly
xi+1
= xi
+ hi+1
, yj+1
= yj
+ rj+1
i = 0,...,M1
-1
j = 0,...,M2
-1
Потоковые точки: xi+ ½
= xi
+
hi+1
, i = 0,1,...,M1
-1
2
yj+ ½
= yj
+
rj+1
, j = 0,1,...,M2
-1
2
Обозначим : U
(xi
,yj
) = Uij
I(xi+½
,yj
) = Ii+½,j
I(xi
,yj+½
) = Ii,j+½
Проинтегрируем уравнение Пуассона: Dj
= - q
(N
d
+ N
a
)
Q(x,y)
по области: V
ij
= { (x,y) : xi- ½
<
x
<
xi+ ½
, yj- ½
<
y
<
yj+ ½
}
xi+ ½
yj+ ½
xi+ ½
yj+ ½
ò ò
Dj
dxdy =
ò ò
Q(x,y)dxdy
xi- ½
yj- ½
xi- ½
yj- ½
Отсюда: yj+½
xi+½
ò
(Ex(xi+½,y) - Ex(xi-½,y) )dx +
ò
(
Ey(x,yj+½) - Ey(x,yj-½)
)
dy=
yj-½
xi-½
xi+ ½
yj+ ½
=
ò ò
Q(x,y)dxdy
xi-
½
yj- ½
Здесь: Ex
(x,y) = - d
j
(x,y)
dx (*)
Ey
(x,y) = - d
j
(x,y)
dy
x
у
-компоненты вектора напряженности электрического поля Е
. Предположим при yj-½
<
y
<
yj- ½
Ex
(xi +
½
,yj
) = Ei+ ½ ,j
= const
yj-½
<
y
<
yj- ½
Ex(xi - ½
,yj
) = Ei- ½ ,j
= const (**)
xi-½
<
x
<
xi+ ½
Ey(xi
,
yj + ½
) = Ei,j+ ½
= const
xi- ½
<
x
<
xi+ ½
yj- ½
<
y
<
yj+ ½
-
Q(x,y) = Qij
= const
Тогда (Ex
)i+ ½ ,j
- (Ex
)i -½ ,j
r*
j
+ (Ey
)ij+ ½
- (Ey
)ij- ½
h*
i
= Qij
h*
i
r*
j
где h*
i
= hi
- hi+1
, r*
j
= rj
- rj+1
2 2
Теперь Еi+
½
,j
выражаем через значение j
(x,y)
в узлах сетки: xi+1
òE
x(x,yj
)dx =
-
j
i+1,j
-
j
ij
xi
из (*
*)
при y=yj
: (Ex
)i+ ½ ,j
= -
j
i+1j
-
j
ij
hi+1
Анологично : (Ey
)i,j+ ½
= -
j
ij+1
-
j
ij
rj+1
Отсюда: (
Dj
)ij
= 1
j
i+1,j
-
j
ij
-
j
i j
-
j
i-1,j
+ 1
j
i j+1
-
j
ij
-
j
ij
-
j
ij-1
=
h*
i
hi+1
hi
r*
j
rj+1
rj
= N
dij
+ N
aij
Граничные условия раздела сред
SiO2
e1
Si y en
x Для области V
0j
yj+ ½
x ½
e
n
e
0
ò
(Ex
(x ½
,y) - E+
x
(0,y))dy +
e
n
e
0
ò
(Ey
(x,yj+ ½
) - Ey
(x,j- ½
))dx =
yj- ½
0
x ½
yj+½
= q
ò
ò
(
N
d + N
a
)
dxdy
0 yj-½
Для области V`
0j
yj+ ½
x ½
e
n
e
0
ò
(E-
x
(0,y) - Ex
(x -½
,y))dy +
e
n
e
0
ò
(Ey
(x,yj+½
) - Ey
(x,j-½
))dx = 0
yj- ½
0
где E+
x
(0,y)
и E-
x
(0,y)
-предельные значения х
компоненты вектора Е
со стороны кремния и окисла.Складывая равенства и учитывая условия: e
n
e
0
d
j
+
-
e
1
e
0
d
j
-
= -Qss
dx dx
имеем yj+½
x½
ò
(
e
n
e
0
Ex
(x½
,y) -
e
1
e
0
Ex
(x-½
,y) - Qss
(y)
)
dy
+
e
n
e
0
ò
(
Ey
(x,yj+½
) +
E
y
(x,
yj-½
)
)
dx
+
yj-½
0
0 x½
yj+½
+
e
1
e
0
ò
(Ey
(x,yj+½
) - Ey
(x,yj-½
))dx = q
ò
ò
(
N
d
+ N
a
)
dxdy
x-½
0 yj-½
Сделав относительно Ex
и Ey
предположения анологичные (**)
положив Qss
(
y
)
= Qss
= const
при yj-½
< y < yj+½
и учитывая условия : j+
=
j-
dj
+
= dj
-
dy dy
“+”- со стороны кремния “-“ - со стороны окисла Получим : e
n
e
0
(Ex
)½,j
-
e
1
e
0
(Ex
)-½,j
- Qss
r*
j
+
e
n
e
0
h1
+
e
1
e
0
h-1
.
(Ey
)0,j+½
- (Ey
)0,j-½
=
2 2 = q (N
d0j
- N
a0j
) h1
r*
j
2
что можно записать : 1
e
n
e
0
j
ij
-
j
0j
-
e
1
e
0
j
0j
-
j
ij
+
e
n
e
0
h1
+
e
1
e
0
h-1
j
0,j+1
-
j
0j
-
j
0j
-
j
0,j-1
=
h*
h1
h-1
2h*
r*
j
rj+1
rj
= - q
( Nd0j
- Na0j
) . h1
- Qss
2 h*
h*
где h*
= h1
+ h-1
2
Общий алгоритм численого решения задачи
Метод установления
Для вычисленя решений многих решений многих многих стационарных задач математической физики, описывающих равновесные состояния, рассматриватривают последнии как результат установленияразвивающегося во времени процесса, расчёт которых оказывается проще, чем прямой расчёт равновесного состояния. Рассмотрим применение метода установления на примере алгоритма для вычисления решения задачи Дирихле: L
xx
Umn
+
L
yy
Umn
=
j
(xm
,yn
) (1)
Umn
|
г =
Y
(smn
) m,n = 1,2,...,M-1
аппроксимирующий дифференциальную задачу Дирихле: d2
U
+ d2
U
=
j
(x,y) 0<= x <=1
dx2
dy2
(2)
U
|
г =
Y
(s) 0<= y <=1
Вслучае задачи (1)
удаётся провести теоретический анализ различных алгоритмов установления с помощью конечных рядов Фурье. Способыточного решения задачи (1)
выдерживающие обобщения на случай переменных коэффициенто и областей скриволинейной границей, например, метод исключения Гаусса , при сколько-нибудь больших и становится неудобным и не применяются. Решение U
(x,y)
Задачи (2)
можно понимать как не зависящую от времени температуру в точке (
x,y
)
пластинки, находящейся в теплолвом равновесии. Функция j
(x,y)
и Y
(s)
означаютв таком случае соответственно распределения источников тела и температуру на границе. Рассмотрим вспомогательную нестационарную задачу о распределении тепла: dV
= d2
V
+ d2
V
-
j
(x,y)
dt dx2
dy2
V
|
г
=
Y
(s) (3)
V
(x,y,0) =
Y
0
(x,y)
где j
иY
те же что и в задаче (2)
, а Y
0
(x,y)
- произвольная. Поскольку источники теплп j
(x,y)
и температура на границе Y
(s)
не зависит от времени, то естественно, что и решение V
(x,y,t)
с течением времени будет менятся всё медленнее, распределение температур V
(x,y,t)
в пределе при t -OO превращается в равновесное распределение тмператур U
(x,y)
, описываемое задачей (2)
. Поэтому вместо стационарной задачи (2)
можно решать нестационарную задачу (3)
до того времени t
, пока её решение перестаёт менятся в пределах интересующей нас точности. В этом состоит идеал решения стационарных задач методом установления. В соответствии с этим вместо задачи (2)
решается задача (3)
, а вместо разностной схемы (1)
для задачи (2)
рассмотрим и составим три различные разностные схемы для задачи (3)
. Именно, рассмотрим простейшую явную разностною схему: Up+1
mn
- Up
mn
=
L
xx
Up
mn
+
L
yy
Up
mn
-
j
(xm
,yn
)
t
Up+1
mn
|
г
=
Y
(smn
) (4)
U0
mn
=
Y0
xm
,yn
)
Рассмотрим так же простейшую неявную разностную схему: Up+1
mn
- Up
mn
=
L
xx
Up+1
mn
+
L
yy
Up+1
mn
-
j
(xm
,yn
)
t
Up+1
mn
|
г
=
Y
(smn
) (5)
U0
mn
=
Y
0
(xm
,yn
)
и исследуем схему применения направлений U’mn
- Up
mn
= 1
[
L
xx
U’mn
+
L
yy
Up
mn
-
j
(xm
,yn
)]
t
2
Up+1
mn
- U’mn
= 1
[
L
xx
U’mn
+
L
yy
Up+1
mn
-
j
(xm
,yn
)]
t
2 (6)
Up+1
mn
|
г
= U’mn
|
г
=
Y
(smn
)
U0
mn
=
Y
0
(xm
,yn
)
Будем считать, что Y
0
(xm
,yn
)
по уже известному Up
={Up
mn
}
для схемы (4)
оссуществляется по уже явным формулам. Вычисление Up+1
= {Up+1
mn
}
по схеме (5)
требует решения задачи : L
xx
Up+1
mn
+
L
yy
Up+1
mn
- Up+1
mn
=
j
(xm
,yn
) - Up
mn
t
t
(7)
Up+1
mn
|
г
=
Y
(smn
)
Вычисление Up+1
= {Up+1
mn
}
по уже известным Up
= {Up
mn
}
по схеме (6)
осуществляется прогонками в направлении оси OX
для вычисления решений {U’mn
}
одномерных задач при каждом фиксированом n
, а затем прогонками в направлнии осиOY
для вычисления решений {Up+1
mn
}
одномерных задач при каждом фиксированом m
. Для каждой из двух разностных схем (4)
и (6)
рассмотрим разность для счёта погрешностеи вычислений: e
p
mn
= Up
mn
- Umn
между сеточной функцией Up
= {Up
mn
}
и точным решением U = {Umn
}
задачи (1)
. Решение {Umn
}
задачи (1)
удовлетворяет уравнениям: Up
mn
- Umn
=
L
xx
Umn
-
j
(xm
,yn
)
t
Umn
|
г
=
Y
(smn
)
U0
mn
= Umn
Вычитая эти равенства из (4)
почленно, получим для погрешности e
p
mn
следующую разностную задачу: e
p+1
mn
|