close

Вход

Забыли?

вход по аккаунту

?

Программный комплекс численного решения задачи о баллистическом диоде.

код для вставкиСкачать
Том 8, №1 (январь - февраль 2016)
Интернет-журнал «НАУКОВЕДЕНИЕ»
publishing@naukovedenie.ru
http://naukovedenie.ru
Интернет-журнал «Науковедение» ISSN 2223-5167 http://naukovedenie.ru/
Том 8, №1 (2016) http://naukovedenie.ru/index.php?p=vol8-1
URL статьи: http://naukovedenie.ru/PDF/15TVN116.pdf
DOI: 10.15862/15TVN116 (http://dx.doi.org/10.15862/15TVN116)
Статья опубликована 12.03.2016.
Ссылка для цитирования этой статьи:
Шевченко А.С. Программный комплекс численного решения задачи о баллистическом диоде // Интернетжурнал «НАУКОВЕДЕНИЕ» Том 8, №1 (2016) http://naukovedenie.ru/PDF/15TVN116.pdf (доступ свободный).
Загл. с экрана. Яз. рус., англ. DOI: 10.15862/15TVN116
УДК 519.688
Шевченко Алеся Сергеевна
ФГБОУ ВПО «Алтайский государственный университет»
Рубцовский институт (филиал), Россия, Рубцовск1
Кандидат физико-математических наук, доцент
E-mail: Ibragimova.a.s@mail.ru
РИНЦ: http://elibrary.ru/author_profile.asp?id=603476
Программный комплекс численного решения задачи
о баллистическом диоде
Аннотация. К настоящему моменту существует довольно много математических
моделей, описывающих с той или иной степенью достоверности физические явления в
полупроводниковых приборах. В процессе поиска приближенных решений задач физики
полупроводников возникает проблема построения численных алгоритмов для этих моделей.
Актуальность конструирования подобных алгоритмов не вызывает сомнений, поскольку на
сегодняшний день полупроводниковые устройства являются неотъемлемой частью многих
электронных приборов.
В данной статье рассматривается, хорошо известная из физики полупроводников,
задача о баллистическом диоде и предлагаются некоторые вычислительные алгоритмы,
построенные в отличие от обычных конечно-разностных схем, на других идеях.
Предложенные алгоритмы основаны на методе установления, применении сглаживающих
регуляризующих операторов и идей схем без насыщения.
В качестве основной математической модели взята недавно предложенная
гидродинамическая модель переноса заряда в полупроводниках в одномерном случае.
С этой целью написан программный комплекс на языке Object Pascal в среде Delphi 7,
реализующий вычислительные алгоритмы.
Приведены подробное описание программного комплекса и результаты численных
экспериментов. Полученные результаты являются физически правдоподобными и в целом
согласуются с результатами других авторов.
Ключевые
слова:
программный
комплекс;
гидродинамическая
модель;
баллистический диод; стационарные решения задачи; уравнение Пуассона; метод
установления; нестационарная регуляризация; схема предиктор-корректор; техника сплайнфункций; интегральные уравнения
1
658225, г. Рубцовск, пр-т Ленина 200-Б
1
http://naukovedenie.ru
15TVN116
Том 8, №1 (январь - февраль 2016)
Интернет-журнал «НАУКОВЕДЕНИЕ»
publishing@naukovedenie.ru
http://naukovedenie.ru
Введение.
Математическое
моделирование
физических
процессов
в
полупроводниковых устройствах имеет большое значение в условиях современного
стремительного развития микроэлектроники, сопровождающегося возрастающим спросом на
повышение функциональности, производительности и улучшение характеристик микросхем и
в последнее время превратилось в быстро развивающуюся область прикладной математики.
К настоящему моменту существует достаточно много математических моделей,
описывающих физические явления в полупроводниковых приборах. Встает вопрос о
конструировании численных алгоритмов и их обосновании для нахождения приближенных
решений таких моделей. В качестве математического описания процесса переноса заряда в
полупроводниках используется гидродинамическая модель (Anile и Romano [1, 2]), которая с
математической точки зрения представляет собой квазилинейную систему уравнений,
записанных в форме законов сохранения, и которая получается из системы моментных
соотношений для уравнения переноса Больцмана с помощью принципа максимума энтропии.
Постановка задачи. Рассмотрим модель переноса заряда в одномерном случае в
безразмерном виде (процесс обезразмеривания подробно описан в работах [3]):
Rt  J x =0,


2

J t   RE  =RQ +c11 J  c12 I ,


3
x

RE

I
=
JQ
+
cP
,
 t x

5
 10
2
I t   RE  = REQ +c21 J  c22 I , 

9
x 3
xx  R  .
(1)
(2)
 − электрический потенциал, R − электронная плотность, E − энергия электрона,
2
J  Ru , I  Rq − потоки, u − электронная скорость, q − поток энергии,   E  1 , P  R ,
3
Q   x ,   0 − диэлектрическая постоянная,    ( x) − плотность легирования (заданная
функция на отрезке [0,1]), коэффициенты c, c11, c12 , c21, c22 − гладкие функции от энергии E ,
выражения для которых приведены, например, в [3, 4].
Для системы (1), (2) сформулируем известную в физики полупроводников задачу о
баллистическом диоде. Это одномерная задача, описывающая полупроводниковый прибор,
разделенный на три части. Первая и третья области представлены высоколегированным
материалом, поэтому они называются n  областями. В средней же части (в канале) − зона
низкой концентрации легировании (так называемая n -область). В качестве характерных
величин L и N  будем рассматривать соответственно ширину n   n  n  − канала и
плотность легирования n  - зоны.
1, если x принадлежи n   зоне,

 
N


, если x принадлежи n  зоне,

N

N плотность легирования n - области. Типичный сглаженный профиль функции  ( x )
изображен на Рис. 1.
2
http://naukovedenie.ru
15TVN116
Том 8, №1 (январь - февраль 2016)
Интернет-журнал «НАУКОВЕДЕНИЕ»
publishing@naukovedenie.ru
http://naukovedenie.ru
Рис. 1. Схематическое представление n   n  n  баллистического диода
Граничные условия для системы (1) при x  0, 1, t  0 , соответствующие известной
физики полупроводников задаче о баллистическом диоде следующие:
R(t , 0)  R(t ,1)  1, 

3
E (t , 0)  E (t ,1)  .
2
(3)
Начальные данные при t  0, 0  x  1 :
R(0, x)  R 0 ( x), 
J (0, x)  J 0 ( x), 

E (0, x)  E 0 ( x), 
I (0, x)  I 0 ( x), 
(4)
R 0 ( x), E 0 ( x)  0 .
Краевые условия для уравнения Пуассона (2) при x  0,1 (t  0) следующие:
 (t , 0)  0,  (t ,1)  B,
(5)
где B − напряжение смещения.
При попытке поиска численных решений этой модели мы столкнулись с рядом
существенных сложностей, обусловленных нелинейностью уравнений модели и наличием в
ней малых параметров. С помощью традиционной теории разностных схем эту задачу решить
не удалось.
Задача (1) − (5) в стационарном случае. В стационарном случае математическая
модель (1), (2) сводится к смешанной задаче, состоящей из трех уравнений Пуассона, с
краевыми условиями на границах x  0, 1 [5-7]:
3
http://naukovedenie.ru
15TVN116
Том 8, №1 (январь - февраль 2016)
Интернет-журнал «НАУКОВЕДЕНИЕ»
publishing@naukovedenie.ru
http://naukovedenie.ru
 d 2
 2
 dx
 d 2
 2
dx
 2
d 
 dx 2




( )

( )
  X  Q    a1 2  a2X  a3Q  a4 XQ  a5Q 2  bc ,

( )
  X  Q   X      X
(    )   (e    ),
(6)
 b5Q 2 
2
 b1  b2X  b3Q  b4 XQ 
2

e     nc ,
1 


     0 при x  0,1,

  0 при x  0,
  B при x  1.

(7)
Компоненты векторов u q находятся из следующих соотношений:
u  F ( E ) Q  1    X  F0 ( E ) ,
(8)
q  G ( E ) Q  1    X  G0 ( E ).
1
  ,   ln R ,


a2  1  1   b G( E)  aF ( E) ,
Здесь
b1  mF ( E ) F0 ( E )  nG( E )G0 ( E ) ,
a5  bF ( E ) ,
b3  
1
1   
d
d
X
,
, a1  aF ( E ) F0 ( E )  bG( E )G0 ( E ) ,
dx
dx
a3  aF ( E )  bG( E )  bF ( E ) F0 ( E ) , a4  b 1    F ( E) ,

2
 mF ( E )  nG ( E )  nF ( E ) F0 ( E ) ,
b4 
b2  1   nG( E)  mF ( E) ,
1
 bF ( E ) F0 ( E ) ,
1 
b5  nF ( E ) ,
2 c21
2 c22
c a
c b
 c11 ,
b  b( E ) 
 c12 ,
m  m( E )  11
n  n( E )  12
,
,
5 1 
5 1 
1 
1 
5
5
5
5
c22  Ec12
c21  Ec11
Ec12
Ec11
3
3
3
3
, G( E )  
, F0 ( E )  1 
, G0 ( E )  1 
,
F (E)  
5
5
det
det
c22  Ec12
c21  Ec11
3
3
da 3 da

и т.д.
det  c11c22  c21c12 , a 
d 2 dE
a  a( E ) 
Также для исходной задачи (1) – (5) в стационарном случае удалось получить систему
интегральных уравнений для определения  , E , I , J при R   , которую можно
воспринимать как приближенную математическую модель для поиска стационарных решений
задачи о баллистическом диоде при малых значениях величины  .
Для поиска стационарных решений задачи (1) – (5) с помощью системы интегральных
уравнений необходимо:
Задать начальные данные для  , E . Начальные данные для  рассчитываются
2
по формуле   E  1 .
3
1.
4
http://naukovedenie.ru
15TVN116
Том 8, №1 (январь - февраль 2016)
Интернет-журнал «НАУКОВЕДЕНИЕ»
publishing@naukovedenie.ru
http://naukovedenie.ru
x
Вычислить
2.
коэффициенты:
K ( x)    c d ,
c, c11, c12 , c21, c22 ,
0
x
a11 ( x)  
c11  c12
0

x
a12 ( x )  
d ,
0
c12

d ,
c21  c22
d ,
 1   
0
x
a21 ( x)  
det  a11 1 a22 1  a12 1 a21 1 .
I0 
0
c22
d ,
 1   
 5


c22
d 2 ( x)    1    
K  d ,
  1    
0
2


 c
d1 ( x)   1     12 K  d ,
  
0
Вычислить
a22 ( x)  
x
x
3.
x
величины:
J
1 
5 

a22 1  d1 1  B   a12 1  d 2 1  B   ,
det 
2 

1 
5 

a21 1  d1 1  B   a11 1  d 2 1  B   .
det 
2 

4.
Вычислить
2
5
значения
2
5
,  :
  x   d1  x   a11  x  J  a12  x  I 0 ,
2
5
  x     x   d 2  x   a21  x  J  a22  x  I 0 ,   x     x     x  , E  x  
3   x   1
2
.
Пункты 2−4 повторяем до тех пор, пока значения неизвестных не будет удовлетворять
заданной точности.
Нестационарная регуляризация краевой задачи (6) − (7). В работах [5−7, 9] для
нахождения решения краевой задачи (6−7) использовался метод установления [8]. С этой
целью, вместо краевой задачи рассматривалась ее нетривиальная нестационарная
регуляризацию Соболева [10], т.е.






 1   2    2 


2
2
 1       

 1   2    2  


( )
(    ),
( )
( ,  ,  ,  ),
( )
( ,  ,  ,  ,  ,  ),
(9)


, 
− дифференциальные операторы, переменная t играет роль времени.
t
x
Замечание. Также рассматривалась нами и параболическая регуляризация [10]
уравнений (6):
   2 

2
    

2
    
( )
(    ),
( )
( )
( ,  ,  ,  ),
( ,  ,  ,  ,  ,  ),
но в ходе численных экспериментов было установлено, что скорость сходимости
метода установления при использовании регуляризации Соболева (9) намного выше скорости
сходимости при применении параболической регуляризации, поэтому в дальнейшем при
разработке алгоритмов используем только систему уравнений (9).
Решение системы уравнений (9) ищем в области t  0, 0  x  1 с краевыми условиями
(7) и начальными данными
5
http://naukovedenie.ru
15TVN116
Том 8, №1 (январь - февраль 2016)
Интернет-журнал «НАУКОВЕДЕНИЕ»
publishing@naukovedenie.ru
http://naukovedenie.ru
  0, x   0  x  ,

2

  0, x    0  x   E0  x   1,
3

   0, x    0  x   ln R0  x  .
(10)
Для нахождения решения смешанной задачи (7), (9) − (10) разработано несколько
алгоритмов, основанные на сведение задачи к интегральным уравнениям, использование
техники сплайн – функций и схемы предиктор – корректор [5-7, 9]. Вкратце опишем эти
алгоритмы.
Модельные задачи. Нетрудно заметить, что в системе (9) уравнения различаются
только правой частью. Поэтому для построения алгоритмов мы использовали модельную
задачу


 1   2  y   2 y  f  y  , t  0, 0  x  1;


 y  0 при x  0,1, t  0;
 y  y x при t  0, 0  x  1,
0 

f   − известная правая часть. Подставляя в (11) ( ) ,
 ,  вместо y соответственно, получим систему уравнений (9).
y
( )
(11)
,
( )
y
вместо f   и  ,
Нам потребуется также другая форма записи модельной задачи (11):
 y   y  f  y  , t  0, 0  x  1;


2

 y   y0  x  =y0  x    y0  x  при t  0, 0  x  1,
(12)
 y   y  t , x  =y  t , x    2 y  t , x  , f  y   f  y   t , x  =y  t , x   f  y   t , x  .
Для нахождения стационарных решений (1) − (5) с помощью метода установления в
y
уравнении 1   2  y   2 y  f   проведём дискретизацию по переменной t. Для этого введём


обозначения: y n  x   y  n, x   y, n  0,1,..., Y  y n1  x  ,  - шаг разностной сетки по
времени t. Аппроксимируя производную  y выражением
 2Y  BY  F, B 
Yy
получаем:

1
 2 y  y  f  y 
, F
.
1 
1 
(13)
Техника сплайн-функций. Приближённое решение модельной задачи (11) ищем в
виде интерполяционного кубического сплайна класса С 2 (см. [11]):
S ( x)  (1   )Yk   Yk 1 
 
h2
 (1   ) (2   )mk  (1   )mk 1  ,
6
(14)
x  xk
, x   xk , xk 1  , k  0, N  1, xk  k  h , N  h  1, Yk  Y ( xk ), mk   2Y ( xk ).
h
Первая производную от кубического сплайна следующая:
6
http://naukovedenie.ru
15TVN116
Том 8, №1 (январь - февраль 2016)
Интернет-журнал «НАУКОВЕДЕНИЕ»
publishing@naukovedenie.ru
http://naukovedenie.ru
 S ( x) 
Вычисляя
 S ( xk  0) 
Yk 1  Yk h
 (2  6  3 2 )mk  (1  3 2 )mk 1  , x   xk , xk 1  , k  0, N  1 .
h
6
 S ( xk  0) ,
 S ( xk  0) ,
где
 S ( xk  0) 
Yk 1  Yk h
  2mk  mk 1  ,
h
6
Yk  Yk 1 h
  mk 1  2mk  и приравнивая их, получаем
h
6
1
1
3
mk 1  2mk  mk 1  2 Yk   Yk  , k  1, N  1,
2
2
h
  1  1
 1Yk  Yk 1  1    .
    1,
−
разностные
операторы,
,
(15)
 1
операторы
сдвига:
Полагая в (13) x  xk и подставляя  2Yk в (15), получаем следующую трехслойную
разностную схему:
 h2 
 h2 
 h2 
1  B Yk 1  2 1  B Yk  1  B Yk 1 
3 
6 
 6 


h2
 Fk 1  4Fk  Fk 1 , k  1, N  1, Fk  F  xk  .
6
(16)
Краевые условия при k  0 и k  N для нахождения функции Yk , k  0, N с помощью
системы уравнений (16) следуют из краевых условий для уравнения (11):
Y0  0, YN  0 .
(17)
Систему (16) − (17) решаем методом прогонки [8].
Замечание. Значения производной  y , которые потребуются для расчётов, в узлах
разностной сетки (т.е. при x  hk  xk , k  0, N ), вычисляются с помощью сплайна (14).
Yk 1  Yk h
 (2  6  3 2 )mk  (1  3 2 )mk 1  , x   xk , xk 1  , k  0, N  1 ,
h
6
Y Y h
Y Y h
то  S ( xk 1 )  k 1 k   mk  2mk 1  , x   xk , xk 1  , k  0, N  1,  S (0)  1 0   2m0  m1 .
h
6
h
6
Так как  S ( x) 
Сведение модельной задачи (12) к системе интегральных уравнений. Задача (12) −
неоднородное дифференциальное уравнение (с параметром x). Решение данного уравнения
следующее:
t
 y  t , x   et y0  x    et    y  , x   f  y   , x  d  ,
0


(18)
где  y0  x  =y0  x    2 y0  x  − начальные данные для соотношения (12).
Из  y  t , x  =y  t , x    2 y  t , x  и y  0 при x  0, 1, t  0 получаем краевую задачу с
параметром t:
7
http://naukovedenie.ru
15TVN116
Том 8, №1 (январь - февраль 2016)
Интернет-журнал «НАУКОВЕДЕНИЕ»
publishing@naukovedenie.ru
http://naukovedenie.ru
 2 y  t , x   y  t , x    y  t , x  ,

 y  t , 0   y  t ,1  0.
(19)
Решение полученной краевой задачи (19) запишем так:
1
y  t , x     G  x, s  y  t , s  ds,
(20)
0
 sh s sh  x  1
, 0  s  x,

sh1
− функция Грина.
G  x, s   
 sh x sh  s  1 , x  s  x.

sh1
Дифференцируя (20) по x, получаем выражение для  y , которое будет использовано
нами для расчётов правых частей уравнений в системе (9):
x
 y  t , x     ch  x  s   t , s ds 
0
chx 1
 ch  s  1  t , s ds .
sh1 0
(21)
Численные расчеты для модельной задачи (11) проводятся c переменным шагом
hi  xi 1  xi , i  0, N 1 на сетке по x. В интегральных уравнениях (18), (20), (21) проводим
 i = y  n, xi  , yin  y  n, xi  ,
дискретизацию по переменной t и вводим обозначения:  y
n
−
шаг разностной сетки по переменной t. В результате переходим к соотношениям:
1
 n
 yi    G  xi , s  y   n  1 , s  ds,
0

xi

chxi 1


y
n

,
x


ch
x

s

n

1

,
s
ds

 
i
  i     
 ch  s  1   n  1 , s ds, (22)
sh
1
0
0


n

n


 y   e n    0, xi    e  y  , xi   f  y   , xi  d   ,
i

 

0

 y  0, x  =y0  xi    2 y0  xi  , n  1, 2,....
Интегралы в системе (22) вычисляем методом трапеции.
Схема предиктор−корректор. Нам понадобятся уравнение (12) и интегральные
соотношения, полученные выше, в которых проводим дискретизацию по переменной t.
Для уравнения (12) схема предиктор − корректор запишется так:
 y*   yn
2
 
y
+ y* = f  
n
,
 yn1   yn

 
*
y
+ yn1 = f   .
Выражая  y* из первого, а  yn 1 − из второго уравнения, получаем
8
http://naukovedenie.ru
15TVN116
Том 8, №1 (январь - февраль 2016)
Интернет-журнал «НАУКОВЕДЕНИЕ»
publishing@naukovedenie.ru
http://naukovedenie.ru
 y*

 yn 
 
  y
f
2

1
2
  
n
,
 yn1
y
 f 

*
1 
n
y
.
Используя уравнения (20) − (21), переходим к следующей системе выражений:
 
n
  y n



f
y
  y n
i
 y
n
*
i
2
f

y

f
n

,
x
,


,


i
i
y
i

i

1
2

1

 yi*    G  xi , s  y*  s  ds,

0

n
 y *


f


*
y
(23)
n 1
  y
i
 y
*
i
f

y

f
n

,
x
,


,


i
i
y

i
i
1 

1
 n 1
 yi    G  xi , s  y   n  1 , s  ds,
0

xi

chxi 1
 y   n  1 , xi     ch  xi  s     n  1 , s ds 
 ch  s  1   n  1 , s ds.
sh1 0

0


 
 
 
 
 
   
Интегралы в системе (23) вычисляем методом трапеции.
Модификация алгоритмов. Теперь рассмотрим модификации алгоритмов. С
вычислительной точки зрения, перспективным оказалось использование для расчёта
потенциала  и его производной следующих соотношений:




1
 n
   n 1  , s 



   s  ds,
 i
 G  xi , s  e

0

x
1
 n, x   i e    n 1 , s    s ds   1  s e    n 1 , s    s ds,
 
 
i

 
 
0
0




 s  xi  1 , 0  s  x,
G  xi , s   

 xi  s  1 , x  s  x.
Решение регуляризованной системы с помощью описанных методов.
Для осуществления расчётов нужно определить начальные значения переменных ( E ,
R ,  ),а затем выполнить следующие шаги:
1.
В начале вычисляем
( )
(  n   ) , где
 n    n, x  − значение

с
предыдущего слоя. Затем находим решение  n1     n  1 , x  первого уравнения задачи (9)
и его производную  n1 на следующем слое.
2.
Рассчитываем
величину
( )
  n, x  ,   n, x  ,   n  1 , x   B,  .
n
Далее решаем второе уравнение задачи (9). Находим  n1     n  1 , x  ,  n 1 .
9
http://naukovedenie.ru
15TVN116
Том 8, №1 (январь - февраль 2016)
Интернет-журнал «НАУКОВЕДЕНИЕ»
publishing@naukovedenie.ru
http://naukovedenie.ru
3.
Вычисляем
( )
  n  1 , x  ,   n, x  ,   n  1 , x   B,
Далее решаем третье уравнение задачу (9). Находим 
n 1
    n  1 , x  ,  n1 .
n 1

,  n,  .
Выполнив 1-3 шаги, алгоритм переходит на следующий временной слой. Данные шаги
выполняются до тех пор, пока не получим решение заданной точности.
При реализации вычислительных алгоритмов, мы использовали идею метода
установления и проводили итерации по временной переменной. В ходе вычислений возникал
скачкообразный рост при некоторых значений переменных задач, обусловленный ее
нелинейностью: норма решений для энергии становилась большой, что вызывало
переполнение буфера и остановку работы программы. Для решения этой проблемы
применили итерации по нелинейности.
Основная идея алгоритма, базирующегося на итерациях по нелинейности, состоит в
том, чтобы пересчитать значения параметров и переменных задач по формулам,
предназначенным для перехода на следующий временной слой, но оставаться при этом на
текущем слое.
Опишем, например, технику сплайн-функций, с применением техники итераций по
нелинейности.
Трехслойная схема (16) преобразуется к виду:
 h 2   m
 h 2   m  h 2   m
1  B Yk 1  2 1  B Yk  1  B Yk 1 
3 
 6 

 6 
h 2 m1

Fk 1  4Fkm1  Fkm11 , k  1, N  1, m  1, M .
6

m 1
Fk

(24)
Здесь M − количество итераций по нелинейности на каждом временном слое, значения
, k  1, N  1 , вычисленные на  m  1 - й итерации по нелинейности.
Алгоритм работы схемы следующий. Значения
Fk0 , k  1, N  1 берутся с предыдущего
m
временного слоя. На m - й итерации по нелинейности вычисляются значения Yk  , k  1, N  1 ,
затем по этим значениям по (24) вычисляются
Fkm и программа переходит на  m  1 -у
итерацию по нелинейности. Как только m  M
временной слой.
программа переходит на следующий
Описание программного комплекса. Для нахождения приближенных решений в
задачи о баллистическом диоде, написан программный комплекс на языке Object Pascal в
среде Delphi 7 [12], реализующий описанные выше вычислительные алгоритмы.
Программный комплекс состоит из следующих модулей:
Unit_main.pas – основной модуль программы.
Unit_matrix.pas − содержит процедуры и функции: сложение, вычитание и умножение
двух матриц; умножение матрицы на число; умножение матрицы на вектор; сложение двух
векторов; вычисление определителя и обратной матрицы; методы Гаусса и прогонки для
решение систем линейных алгебраических уравнений; вычисление первых производных с
помощью кубического сплайна класса С 2 , вычисление вторых производных.
Unit_c.pas − содержит значения физических параметров для баллистического диода;
константы электронно-фононного взаимодействия, энергии фононов; функции для
10
http://naukovedenie.ru
15TVN116
Интернет-журнал «НАУКОВЕДЕНИЕ»
http://naukovedenie.ru
Том 8, №1 (январь - февраль 2016)
publishing@naukovedenie.ru
вычисления обычных и модификационных функций Бесселя; процедуры для вычисления
коэффициентов c, c11, c12 , c21, c22 и их производных; процедуры для вычисления параметров,
представленных выше алгоритмов.
Unit_decision.pas – содержит процедуры реализации предложенных алгоритмов.
Главная форма программного комплекса содержит меню из 3 пунктов: «Файл»,
«Численное решение», «Справка» (см. рис. 2).
Меню «Файл» содержит подменю: «Считать данные из файла», «Сохранить в файл»,
«Выход».
Рис. 2. Главная форма программного комплекса
При выборе пункта меню «Численное решение» появляется окно «Численное решение
задачи о баллистическом диоде», см. рис. 3.
11
http://naukovedenie.ru
15TVN116
Интернет-журнал «НАУКОВЕДЕНИЕ»
http://naukovedenie.ru
Том 8, №1 (январь - февраль 2016)
publishing@naukovedenie.ru
Рис. 3. Окно расчета численного решения задачи о баллистическом диоде
12
http://naukovedenie.ru
15TVN116
Том 8, №1 (январь - февраль 2016)
Интернет-журнал «НАУКОВЕДЕНИЕ»
publishing@naukovedenie.ru
http://naukovedenie.ru
Чтобы найти стационарные решения задачи о баллистическом диоде необходимо в
программе выполнить несколько шагов.
Шаг 1. Задать параметры:
a)
физические параметры, характеризующие исходную задачу;
б)
численные параметры, определяющие вычислительный процесс;
в)
начальные данные для искомых функций  , E , R , I , J .
Начальные значения переменных задачи лучше задавать равным значениям этих
переменных в состоянии глобального термодинамического равновесия.
Необходимые для расчетов физические параметры и их описание содержаться в табл.
1.
Таблица 1
Значение физических параметров
Величина
Описание
N
Плотность легирования n  - области
N
Плотность легирования n - области
L
Длина n   n  n  - канала


V
Значение
1
5 1017 3
сm
1
2 1015 3
сm
7
6 10 m или 3 107 m
10749.32 или 26873.33
Значение функции   x  в n -зоне
баллистического диода
Напряжение смещения
119170 105
m
s
0.5  2Volt
Функция  в расчетах является кусочно-линейной версией функции плотности
легирования   x  (см. рис. 1) и задается в программе так:
1, если 0  x  d  1 ,

   x  , если d  1  x  d  1 ,
  x   1
 , если 1  d  1  x  1  d  1 ,
1, если 1  d    x  1,

1
15
1
2
1 
xd
,
,x 
1     x5  x3  x  
16
3
2
1
5

x  1  d 
15
2
1
 1 
2  x   1     x5  x3  x  
,x 
.
16
3
2
1
5

1  x   
Безразмерная величина B вычисляется автоматически, как только задали значение
напряжение смещения V .
В табл. 2 приведены необходимые для расчетов численные параметры и диапазоны их
значений.
Алгоритмы прекращают работу, если достигается необходимая точность 1 :
13
http://naukovedenie.ru
15TVN116
Интернет-журнал «НАУКОВЕДЕНИЕ»
http://naukovedenie.ru
Том 8, №1 (январь - февраль 2016)
publishing@naukovedenie.ru
  Rin1  Rin   in1   in  in1  in   1.
N
i 0
Таблица 2
Значение численных параметров
Величина
N
t
M
Описание
Количество узлов сетки по x
Шаг разностной сетки по t
Количество итераций по нелинейности
Значение
100  500
0.02  105
1 10
1
1
или
6
3
d
1
Точность вычислений
1
Напряжение смещения
1 1

8 12
После того как указали все необходимые параметры необходимо нажать кнопку
«Утвердить начальные данные».
Шаг 2. Выбор алгоритма решения задачи (интегральные уравнения, модификация
интегральных уравнений, схемы предиктор-корректор, техника сплайн-функций,
приближенная математическая модель), т.е. нужно выбрать соответствующую вкладку с
алгоритмом решения задачи.
Управление расчетом осуществляется с помощью кнопок «Начать расчет» и «Прервать
расчет» на выбранной вкладке.
Кнопка «Начать расчет» становится доступной только после того как будет нажата
кнопка «Утвердить начальные данные».
При нажатии кнопки «Начать расчет», ее название меняется на «Идет расчет».
Кнопка «Прервать расчет» становится доступной только после того как будет нажата
кнопка «Начать расчет».
Все полученные результаты отражаются на вкладках: «Данные расчета», «График
плотности», «График зависимости энергии E от t», «График приближенного решения для
потенциала», «График приближенного решения для энергии E», «График приближенного
решения для электронной плотности R», «График приближенного решения для I», «График
приближенного решения для J», «График приближенного решения для u» и «График
скорректированного u». Только для алгоритма «Приближенная математическая модель» есть
вкладка «Вольт-амперная характеристика баллистического диода».
Вкладка «График зависимости энергии E от t» (наибольшее значение энергии на n-ом
временном слое) позволяет нам следить за поведением энергии (см. Рис. 4).
Кнопки «Сохранить результаты расчета в файл» (пункт меню «Файле» → «Сохранить в
файл» → «Результаты расчета»), «Считать данные из файла» (пункт меню «Файле» →
«Считать данные из файла») позволяют осуществлять сохранение и загрузку результатов
расчета (приближенные решения для электрического потенциала, плотности электронов,
энергии электронов, электрического тока и потока энергии). Данные сохраняются в файлах c
расширением .txt, что позволяет использовать полученные результаты в качестве начальных
условий для дальнейших расчетов.
14
http://naukovedenie.ru
15TVN116
Интернет-журнал «НАУКОВЕДЕНИЕ»
http://naukovedenie.ru
Том 8, №1 (январь - февраль 2016)
publishing@naukovedenie.ru
Графики полученных решений можно сохранить в любом графическом формате с
помощью панели «Сохранение графиков» или пункта меню «Файле» → «Сохранить в файл»
→ «Графики».
Рис. 4. Вкладка «График зависимости энергии E от t»
Результаты тестовых расчетов. Приведем теперь графики численных решений,
полученные в результате расчетов.
На рис. 5. приведены графики энергии, потенциала и вольт-амперная характеристика
баллистического диода (зависимости силы тока от напряжения смещения B).
Рис. 5. Результаты расчетов для случая   0 с параметрами V  1.5Volt ,   0.004 ,
L  6 107 m , N  100 , 1  106
Расчеты проводились с помощью различных алгоритмов.
15
http://naukovedenie.ru
15TVN116
Интернет-журнал «НАУКОВЕДЕНИЕ»
http://naukovedenie.ru
Том 8, №1 (январь - февраль 2016)
publishing@naukovedenie.ru
Вычислительный алгоритм для системы интегральных уравнений и его модификация
дают хорошие результаты при малом   0.004 .
Приближенные решения, полученные с помощью модификации интегральных
уравнений, представлены на Рис. 6.
Рис. 6. V  1.5 Volt ,   0.004 , L  3107 m , N  100 , 1 
1
, 1  106
12
Другие подходы (с использованием техники сплайн-функций и схемы предикторкорректор) дают хорошие результаты при не больших  .
В результате численных экспериментов было выявлено, что скорость сходимости
алгоритма, использующего технику сплайн-функций выше, чем у алгоритма, использующий
схему предиктор-корректор.
C целью верификации предложенных алгоритмов для задачи (9) проводились расчеты с
использованием метода ортогональной прогонки, который позволяет в принципе проводить
расчеты с гарантированной точностью. На рис. 7. приведены результаты сразу трех
алгоритмов (схема предиктор-корректор, техника сплайн функций и метод ортогональной
прогонки) и видно, что результаты практически совпадают.
16
http://naukovedenie.ru
15TVN116
Интернет-журнал «НАУКОВЕДЕНИЕ»
http://naukovedenie.ru
Том 8, №1 (январь - февраль 2016)
publishing@naukovedenie.ru
Рис. 7. V  1 Volt ,   0.2 , L  6 107 m , N  200 , 1 
1
, 1  106
12
Заключение. Проведено множество численных расчетов. Полученные результаты
являются физически правдоподобными и в целом согласуются с результатами других авторов.
Данные результаты носят как теоретический, так и вычислительный характер. Данная работа
представляет интерес для специалистов в области вычислительной математики, может быть
использована для моделирования физических процессов в полупроводниковых устройствах и
решения различных прикладных задач.
17
http://naukovedenie.ru
15TVN116
Том 8, №1 (январь - февраль 2016)
Интернет-журнал «НАУКОВЕДЕНИЕ»
publishing@naukovedenie.ru
http://naukovedenie.ru
ЛИТЕРАТУРА
1.
Anile A.M., Romano V. Non parabolic band transport in semiconductors: closure of
the moment equations // Cont. Mech. Thermodyn, 1999. ‒ V.11. ‒ P. 307–325.
2.
Romano V. 2D simulation of a silicon MESFET with a non-parabolic hydrodynamical
model based on the maximum entropy principle // J. Comp. Phys, 2002. ‒ V.176. ‒ P.
70–92.
3.
Blokhin A.M., Bushmanov R.S., Romano V. Nonlinear asymptotic stability of the
equilibrium state for the MEP model of charge transport in semiconductors, //
Nonlinear Analysis, 2006. ‒ V. 65. ‒ P. 2169‒2191.
4.
Blokhin A.M., Ibragimova A.S. Numerical method for 2D Simulation of a Silicon
MESFET with a Hydrodynamical Model Based on the Maximum Entropy Principle //
SIAM Journal on Scientific Computing, – 2009. – V.31. – Issue 3. – P. 2015–2046.
5.
Blokhin A.M., Ibragimova A.S. 1D Numerical Simulation of the MEP Mathematical
Model in ballistic diode problem // Journal of Kinetic and Related Models, − 2009. −
V. 2. − N.1. − P. 81−107.
6.
Блохин А.М., Ибрагимова А.С., Семисалов Б.В. Конструирование
вычислительных алгоритмов для задачи о баллистическом диоде //
Вычислительная математика и математическая физика, – 2010. – Т.50. – №1. – c.
1−21.
7.
Blokhin A.M., Ibragimova A.S., Semisalov B.V. Design of Numerical Algorithms for
the Ballistic Diode Problem // Computational Mathematics and Mathematical Physics,
– 2010. – V. 50. – N. 1. – P. 180 – 200.
8.
Бабенко К.И. Основы численного анализа. ‒ Москва‒Ижевск: НИЦ "Регулярная
и хаотическая динамика", 2002. ‒ 848 c.
9.
Блохин А.М., Ибрагимова А.С., Семисалов Б.В. Численный анализ задач
переноса заряда в полупроводниковых устройствах. Saarbrucken, Germany:
Palmarium Academic Publishing, 2012. – 209 c.
10.
Шевченко А.С. Алгоритм поиска приближенных решений уравнения Пуассона //
Молодой ученый, 2015. − №3. − С. 18-23.
11.
Завьялов Ю.С., Квасов Б.И., Мирошниченко В.Л. Методы сплайн‒функций.
−М.: Наука, 1980. – 355 c.
12.
Белов В.В., Чистяков В.И. Программирование в Delphi: процедурное, объектноориентированное, визуальное. Учебное пособие для вузов. − М.: Горячая линияТелеком, 2014. ‒ 240 с.
18
http://naukovedenie.ru
15TVN116
Интернет-журнал «НАУКОВЕДЕНИЕ»
http://naukovedenie.ru
Том 8, №1 (январь - февраль 2016)
publishing@naukovedenie.ru
Shevchenko Alesya Sergeevna
Altay State University
Rubtsovsk branch, Russia, Rubtsovsk
E-mail: Ibragimova.a.s@mail.ru
Software complex of numerical solution of the ballistic
diode problem
Abstract. At present there are quite a few of mathematical models describing physical
phenomena in semiconductor devices with one or another degree of reliability. The problem of
constructing numerical algorithms for these models arises in process of finding approximate
solutions of semiconductor physics. Undoubtedly, the designing of such algorithms is actual, because
today semiconductor devices are essential parts of many electronic devices.
In this paper, we consider the ballistic diode problem, which is well known in the physics of
semiconductors, and propose computational algorithms based on ideas other than those underlying
usual finite-difference schemes. The proposed algorithm is based on the stabilization method, the
application of regularized smoothing operators and ideas of schemes without saturation.
As the basic mathematical model, we use a recently proposed hydrodynamic model of charge
transport in semiconductors in one-dimensional case.
The computational algorithms are implemented as a software package using Object Pascal in
Delphi 7 environment.
The detailed description of a software complex and results of numerical experiments are
provided. The obtained results are physically plausible and, in general, agree with the results of other
authors.
Keywords: software complex; hydrodynamic model; ballistic diode; stationary solutions of
problem; the Poisson equation; stabilization method; nonstationary regularization; predictorcorrector scheme; spline functions technique; integral equations
19
http://naukovedenie.ru
15TVN116
Том 8, №1 (январь - февраль 2016)
Интернет-журнал «НАУКОВЕДЕНИЕ»
publishing@naukovedenie.ru
http://naukovedenie.ru
REFERENCES
1.
Anile A.M., Romano V. Non parabolic band transport in semiconductors: closure of
the moment equations // Cont. Mech. Thermodyn, 1999. ‒ V.11. ‒ P. 307–325.
2.
Romano V. 2D simulation of a silicon MESFET with a non-parabolic hydrodynamical
model based on the maximum entropy principle // J. Comp. Phys, 2002. ‒ V.176. ‒ P.
70–92.
3.
Blokhin A.M., Bushmanov R.S., Romano V. Nonlinear asymptotic stability of the
equilibrium state for the MEP model of charge transport in semiconductors, //
Nonlinear Analysis, 2006. ‒ V. 65. ‒ P. 2169‒2191.
4.
Blokhin A.M., Ibragimova A.S. Numerical method for 2D Simulation of a Silicon
MESFET with a Hydrodynamical Model Based on the Maximum Entropy Principle //
SIAM Journal on Scientific Computing, – 2009. – V.31. – Issue 3. – P. 2015–2046.
5.
Blokhin A.M., Ibragimova A.S. 1D Numerical Simulation of the MEP Mathematical
Model in ballistic diode problem // Journal of Kinetic and Related Models, − 2009. −
V. 2. − N.1. − P. 81−107.
6.
Blokhin A.M., Ibragimova A.S., Semisalov B.V. Konstruirovanie vychislitel'nykh
algoritmov dlya zadachi o ballisticheskom diode // Vychislitel'naya matematika i
matematicheskaya fizika, – 2010. – T.50. – №1. – c. 1−21.
7.
Blokhin A.M., Ibragimova A.S., Semisalov B.V. Design of Numerical Algorithms for
the Ballistic Diode Problem // Computational Mathematics and Mathematical Physics,
– 2010. – V. 50. – N. 1. – P. 180 – 200.
8.
Babenko K.I. Osnovy chislennogo analiza. ‒ Moskva‒Izhevsk: NITs "Regulyarnaya i
khaoticheskaya dinamika", 2002. ‒ 848 c.
9.
Blokhin A.M., Ibragimova A.S., Semisalov B.V. Chislennyy analiz zadach perenosa
zaryada v poluprovodnikovykh ustroystvakh. Saarbrucken, Germany: Palmarium
Academic Publishing, 2012. – 209 c.
10.
Shevchenko A.S. Algoritm poiska priblizhennykh resheniy uravneniya Puassona //
Molodoy uchenyy, 2015. − №3. − S. 18-23.
11.
Zav'yalov Yu.S., Kvasov B.I., Miroshnichenko V.L. Metody splayn‒funktsiy. −M.:
Nauka, 1980. – 355 c.
12.
Belov V.V., Chistyakov V.I. Programmirovanie v Delphi: protsedurnoe, ob"ektnoorientirovannoe, vizual'noe. Uchebnoe posobie dlya vuzov. − M.: Goryachaya liniyaTelekom, 2014. ‒ 240 s.
20
http://naukovedenie.ru
15TVN116
Документ
Категория
Без категории
Просмотров
4
Размер файла
1 365 Кб
Теги
комплекс, решение, баллистического, программное, диод, задачи, численного
1/--страниц
Пожаловаться на содержимое документа