next up previous
Next: 5 Behandlung von Grenzflächenmodellen Up: 4 Behandlung von Volumsmodellen Previous: 4.3 Quasi-FermienergieFermienergie

4.4 Transformation der Differentialgleichungen

 

Um die Kopplungen der Halbleitergleichungen zu verringern, wird in [1] eine Transformation für das stationäre Drift-Diffusionsmodell vorgeschlagen, welche eine Diagonalisierung nach führenden Differentialtermen zweiter Ordnung bedingt. Diese Transformation wird hier für das hydrodynamische Modell erweitert. Die Transformation für das Drift-Diffusionsmodell ergibt sich dann als Spezialfall.

Die fünf nichtlinearen Differentialgleichungen (4.14, 4.17, 4.18, 4.21, 4.22) werden zum nichtlinearen Differentialgleichungssystem

equation1527

zusammengefaßt. Das entspricht dem nichtlinearen Gleichungssystem (3.1) im Kapitel 3 bevor die Diskretisierung durchgeführt wird. Für das Verfahren nach Newton muß das Gleichungssystem

  eqnarray1532

gelöst werden.

Die Komponenten der Matrix

equation1542

sind die partiellen Ableitungen nach den Variablen, die mit den Ortsableitungen vertauscht werden dürfen. Die Komponenten von tex2html_wrap_inline8760 sind daher Differentialoperatoren bezüglich des Ortes, die auf die Komponenten des Vektors in (4.43), der die linearen Abweichungen enthält, angewendet werden. tex2html_wrap_inline8762 kann als Summe einer Teilmatrix tex2html_wrap_inline8764 , welche die Differentialoperatoren zweiter Ordnung enthält, und tex2html_wrap_inline8766 welche Terme mit Differentialoperatoren erster Ordnung enthält,

equation1603

geschrieben werden. Ziel ist es nun, die Matrix tex2html_wrap_inline8774 mit Hilfe einer Transformationsmatrix tex2html_wrap_inline8552 in Diagonalform zu bringen,

equation1610

Nachdem für die Einträge der Matrix tex2html_wrap_inline8774 nur die Ableitungen zweiter Ordnung benötigt werden, können alle Größen in den Differentialgleichungen mit Ausnahme von tex2html_wrap_inline8652 , n, p, tex2html_wrap_inline8790 und tex2html_wrap_inline8792 zufolge der Regel tex2html_wrap_inline8794 als konstant angenommen werden. Damit ergibt sich folgendes Gleichungssystem tex2html_wrap_inline8796 :

eqnarray1617

mit den Flußdichten

eqnarray1636

Damit ergibt sich unter Beachtung, daß eine zeilenweise Skalierung von tex2html_wrap_inline8774 keinen Einfluß auf die zu bestimmende Transformationsmatrix tex2html_wrap_inline8552 hat,

  equation1690

Daraus folgt die Diagonalisierung

equation1713

und weiter, daß die Matrix tex2html_wrap_inline8558 der Matrix in (4.58) entspricht, in der die Operatoren tex2html_wrap_inline8810 durch 1 ersetzt sind,

equation1742

mit

equation1774

Daraus folgt tex2html_wrap_inline8552 für das hydrodynamische Modell:

Für das Drift-Diffusionsmodell reduziert sich die Matrix tex2html_wrap_inline8558 auf den linken oberen tex2html_wrap_inline8822 -Block:

equation1825

Daraus folgt tex2html_wrap_inline8552 für das Drift-Diffusionsmodell:


next up previous
Next: 5 Behandlung von Grenzflächenmodellen Up: 4 Behandlung von Volumsmodellen Previous: 4.3 Quasi-FermienergieFermienergie

IUE WWW server
Mon Jul 1 14:21:36 MET DST 1996