Реферат: Методы численного моделирования МДП-структур
В данной работе рассматриваются модели МДП-структур, относящиеся по введённой классификации к первому классу. В этом методе производные
неизвестных функций, входящие в исходные дифференциальные уравнения и краевые условия, заменяются конечно-разностными отношениями (построение разностной схемы), в результате чего получается система алгебраических уравнений, которая затем решается прямыми или итерационными методами.
3.1. Алгебраизация ФСУ,
На первом этапе решения системы дифференциальных уравнений необходимо осуществить алгебраизацию задачи путём аппроксимации на сетке множества точек, которыми моделируется область изменения неизвестных.
Каждое из трёх основных уравнений математической модели в интегральной форме выражает закон, который выполняется как в элементарной ячейке, так и во всей области определения,что является следствием фундаментальных физических свойств непрерывности электрического смещения и тока.конечно-разностная схема предполагает сохранение этих свойств и для алгебраических уравнений.
Рассмотрим некоторые вопросы касающиеся построения сеток дискретизации[2]. Соображения удобства реализации алгоритма решения основной системы на ЭВМ, а так же требование его экономичности обуславливают применение регулярных сеток, расположение узлов в которых подчиняется определённым закономерностям. В практике численного моделирования микроэлектронных структур примеяются как непрерывные прямоугольные (неравномерные), так и треугольные сетки (рис.2.). Треугольная сетка позволяет с меньшим количеством дополнительных узлов сгущать сетку в областях локальных неоднородностей (рис.2.б).
При автоматическом построении сетки нужно знать где необходимо сгущение узлов и как интерполировать различные величины для вновь введённых точек сетки. Пространственная сетка должна быть такой, чтобы ошибка дискретизации была распределена по ней равномерно, т.е чтобы частные производные по пространству аппроксимировались с заданной точностью. Метод конечных разностей наиболее удобно реализуется на непрерывных
        |
|
|
   |
 |
    |
а)
 |
           |
|
|
                  |
    |
     |
б)
Рис.2.Виды сеток и
их локальные уточнения.
прямоугольных сетках.Он является точным, если значения величин в каждой точке сетки могут быть описаны полиномом второго порядка.
3.1.1 Дискретизация уравнения Пуассона
В настоящее время для конструирования разностных схем, аппроксимирующих исходные дифференциальные уравнения, применяют различные методы. Например, с помощью интегро-интерполяционного метода, предложенного А.Н.Тихоновым и А.А.Самарским естественным образом можно получить следующую схему [1].
В области V={(x,y),0<x<Lx ,0<y<Ly } (рис.1.) вводится неравномерная сетка
x0 =xi +hi+1 , yj+1 =yj +rj+1
hi ,rj -шаги сетки.
Проинтегрировав уравнение Пуассона [2], представленного в интегральной форме (1.53), по области Vij ={(x,y),xi-1/2 <=x<=xi+1 , yj-1/2 <=y<=yj+1/2 } (рис.3.)
получим разностный аналог уравнения Пуассона в точках, принадлежащих полупроводниковой среде,в виде :
 |
 |
|
|
rj * +hi * =rk hi * rj * (3.11)
Здесь приняты обозначения :
rk =pk -nk +Nk , rj * =(rj-1 +rj )/2, hi * =(hi-1 +hi ),
В узлах сетки,лежащих в диэлектрике, разностный аналог полевого уравнения совподает с (3.11) если положить rk =0.
Если поделить обе части уравнения на -hi * rj * то можнопереписать левую часть уравнения при помощи разностного оператора Dhr [1], тогда разностный аналог уравнение Пуассона примет вид :
 |
Рис.3. Ячейка алгебраизации уравнения Пуассона.
 |
Рис.4. Ячейка алгебраизации уравнения непрерывности.
(Dhr j)k =-rk (*)
Схемы такого типа выражают на сетке законы сохранения и называются консервативно разностными схемами.
3.1.2.Дискретизация уравнения непрерывности.
Аппроксимируя интегральную форму электронного уравнения непрерывности (1.51) на элементарной ячейке Vi (рис.4.), отвечающей лежащему в полупроводниковой среде k-му узлу сетки пространственной дискретизации
|