PDE-Löser für Drift-Diffusion und verwandte Modelle

Ich versuche, grundlegende Halbleitermodelle für pädagogische Zwecke zu simulieren - ausgehend vom Drift-Diffusionsmodell. Obwohl ich keinen handelsüblichen Halbleitersimulator verwenden möchte - ich werde andere (gängige, aktuelle oder undurchsichtige) Modelle erlernen, möchte ich einen handelsüblichen PDE-Löser verwenden.

Aber auch für den einfachen 1D-Fall besteht das Driftdiffusionsmodell aus einer Reihe gekoppelter nichtlinearer PDEs:

Stromdichtegleichungen

J_{n} = q n (x) μ_{n} E (x) + q D_{n} \nabla n

$J_n = q n(x) \mu_n E(x) + qD_n \nabla n$

J_{p} = q p (x) μ_{p} E (x) + q D_{p} \nabla p

$J_p = q p(x) \mu_p E(x) + qD_p \nabla p$

Kontinuitätsgleichung

\frac{\partial n}{\partial t} = \frac{1}{q} \nabla \cdot J_{n} + U_{n}

$\frac{\partial{n}}{\partial{t}} = \frac{1}{q} \nabla \cdot J_n + U_n$

\frac{\partial p}{\partial t} = \frac{1}{q} \nabla \cdot J_{p} + U_{p}

$\frac{\partial{p}}{\partial{t}} = \frac{1}{q} \nabla \cdot J_p + U_p$

Poissongleichung weiterempfehlen

\nabla \cdot (ϵ \nabla V) = - (p - n + N_{D}^{+} - N_{A}^{-})

$\nabla \cdot (\epsilon \nabla V) = -(p - n + N_D^+ - N_A^-)$

und eine Reihe von Randbedingungen.

Ich habe einige Python-FEM-Löser, FEniCS / Dolfin und SfePy ausprobiert , , aber ohne , weil ich sie nicht in der schwachen Variationsform mit Testfunktionen formulieren konnte.

Es gibt natürlich die Möglichkeit, die numerische Lösung von Grund auf neu zu implementieren, aber ich habe mich noch nicht eingehend mit FEM / Numerical befasst. Ich hoffe, dass dies nicht meine einzige Option ist, da ich nicht mit numerischen Problemen überfordert sein möchte.

Gibt es also ein Paket (vorzugsweise Open Source), das diese Gleichungen in dieser Form aufnimmt und löst? Oder ist die von den Werkzeugen geforderte Variationsform vielleicht nicht so schwer? Welche Möglichkeiten habe ich auf jeden Fall?

Vielen Dank

Edit: Versuch, die schwache Variationsform für FEniCS / Dolfin oder SfePy zu formulieren

Unter Verwendung von drei PDEs (Poisson + zwei Kontinuitätsgleichungen mit J substituiert) suchen wir nach V, n und p. Die Poisson-Gleichung (unter Verwendung einer Testfunktion ) ist einfach. Ich habe jedoch Schwierigkeiten mit den Kontinuitätsgleichungen. $u_V$

Die zweite PDE (starke Form)

\frac{\partial n}{\partial t} = \nabla \cdot (C_{1} n \nabla V + C_{2} \nabla n) + U

$\frac{\partial{n}}{\partial{t}} = \nabla \cdot (C_1 n \nabla V +C_2\nabla n) + U$

C_{1}, C_{2}

$C_1, C_2$

U, n, p, V

$U, n, p, V$

$f_n$

\int_{Ω} f_{n} \frac{n - n_{1}}{Δ t} d Ω - C_{1} \int_{Ω} f_{n} \nabla \cdot (n \nabla V) d Ω - C_{2} \int_{Ω} f_{n} \nabla^{2} n d Ω - \int_{Ω} f_{n} U d Ω

$\int_{\Omega}f_n\frac{n - n_1}{\Delta t}d\Omega - C_1 \int_{\Omega}f_n \nabla \cdot ( n \nabla V) d\Omega - C_2 \int_{\Omega}f_n \nabla^2 n d\Omega - \int_{\Omega}f_n U d\Omega$

C_{1} \int_{Ω} f_{n} \nabla \cdot (n \nabla V) d Ω

$C_1 \int_{\Omega} f_n \nabla \cdot ( n \nabla V) d\Omega$

$\mathbf{\nabla V}$ $V, u_n, n$ $\nabla \cdot \phi \mathbf{A} = \mathbf{A} \cdot \mathbf{\nabla \phi} + \phi \nabla \cdot \mathbf{A}$

C_{1} \int_{Ω} f_{n} \nabla \cdot (n \nabla V) d Ω = C_{1} \int_{Ω} f_{n} (\nabla V \cdot \nabla n) + C_{1} \int_{Ω} f_{n} n \nabla \cdot \nabla V

$C_1 \int_{\Omega} f_n \nabla \cdot ( n \nabla V) d\Omega = C_1 \int_{\Omega} f_n (\nabla V \cdot \nabla n) + C_1 \int_{\Omega} f_n n \nabla \cdot \nabla V$

Da V durch die Poisson-Gleichung gelöst wird, können wir den kürzlich berechneten Wert verwenden, der in der Software Dolfin / FEniCS zulässig ist, und vereinfachen, wie wir V in dieser zweiten gekoppelten Gleichung behandeln? Diese Art von Techniken funktioniert beim Diskretisieren (z. B. Gummel, ...), was ich mit diesen fertigen Lösungsmitteln nicht mache!

$J_n$ $n$ $J_n, J_p, n, p, V$ $J_n$

pde Weaam
quelle

Warum kannst du ihre schwachen Formen nicht aufschreiben?

Bill Barth

@ BillBarth Ich habe meine Frage bearbeitet, bitte schauen Sie. Vielen Dank.

Weaam

Ihre Teilintegration ist falsch. Überprüfen Sie die Formel, es fehlen Zeichen, Sie haben rechts mehr Ableitungen als links und Sie haben das Randintegral vergessen.

Wolfgang Bangerth

u_{n}

$u_n$

Ja, ich hätte vorsichtiger sein sollen. Bitte überprüfe meine Bearbeitung, insbesondere meine Frage, wie wir V behandeln, da dies bereits durch die vorherige PDE hätte gelöst werden sollen. Hat dies Auswirkungen auf die Variationsform? Vielen Dank.

Weaam

Die Scharfetter-Gummel (SG) -Formulierung wird üblicherweise zur Lösung der Stromdichtegleichungen verwendet. Dies ist eine spezielle Formulierung, die Schwierigkeiten beim Lösen der nichtlinearen Abhängigkeit zwischen Potential und Stromdichte überwindet.

Ein Standardtext, der beschreibt, wie diese Gleichungen unter Verwendung von Box-Integrationsmethoden angewendet werden, ist in diesem Buch zu finden: Selberherr, S., Analyse und Simulation von Halbleiterbauelementen. Springer-Verlag 1984

Diese Art der Simulation wird als Technology Computer Aided Design (TCAD) bezeichnet. Im Gegensatz zur Finite-Elemente-Methode (FEM) wird zur Berechnung der Ströme die Finite-Volumen-Methode (FVM) verwendet. Dies liegt daran, dass es in die SG-Formulierung passt, von der gezeigt wurde, dass sie (von Praktikern dieser Methode) funktioniert, wenn die Stromdichtegleichungen gelöst werden.

Wenn Sie dies mit generalisierten PDEs lösen möchten, verfügt COMSOL über ein Halbleitermodul, das dieses Problem mit einer hybriden FEM / FVM-Methode löst.

Außerdem sind hier kommerzielle und Open Source-TCAD-Simulatoren aufgeführt: http://www.tcadcentral.com

Meines Wissens nach sind generalisierte PDE-TCAD-Löser DEVSIM, FLOOPS, PROPHET. Bei den kommerziellen Tools sind die meisten physikalischen Gleichungen in einer kompilierten Sprache wie C ++ hartcodiert.

Juan
quelle

Ich entschuldige mich für die extrem späte Antwort. Ich erkannte, dass eine solche direkte Anwendung von DD (auch mit SG) ziemlich instabil war (zumindest meine Implementierung in Fenics), und gab sie daher auf. In einem späteren VLSI-Kurs habe ich tatsächlich Comsol- und TCAD-Tools verwendet. Vielen Dank für Ihre umfassende Antwort.

Weaam,

PDE-Löser für Drift-Diffusion und verwandte Modelle

Antworten: