Mixed finite-element methods for elliptic convection-dominated problems arising in semiconductor physics

In this work we develop and analyze an adaptive numerical scheme for simulating a class of macroscopic semiconductor models. At first the numerical modelling of semiconductors is reviewed in order to classify the Energy-Transport models for semiconductors that are later simulated in 2D. In this class of models the flow of charged particles, that are negatively charged electrons and so-called holes, which are quasi-particles of positive charge, as well as their energy distributions are described by a coupled system of nonlinear partial differential equations. A considerable difficulty in simulating these convection-dominated equations is posed by the nonlinear coupling as well as due to the fact that the local phenomena such as "hot electron effects" are only partially assessable through the given data. The primary variables that are used in the simulations are the particle density and the particle energy density. The user of these simulations is mostly interested in the current flow through parts of the domain boundary - the contacts. The numerical method considered here utilizes mixed finite-elements as trial functions for the discrete solution. The continuous discretization of the normal fluxes is the most important property of this discretization from the users perspective. It will be proven that under certain assumptions on the triangulation the particle density remains positive in the iterative solution algorithm. Connected to this result an a priori error estimate for the discrete solution of linear convection-diffusion equations is derived. The local charge transport phenomena will be resolved by an adaptive algorithm, which is based on a posteriori error estimators. At that stage a comparison of different estimations is performed. Additionally a method to effectively estimate the error in local quantities derived from the solution, so-called "functional outputs", is developed by transferring the dual weighted residual method to mixed finite elements. For a model problem we present how this method can deliver promising results even when standard error estimator fail completely to reduce the error in an iterative mesh refinement process.

Hypersurfaces with many singularities

1. Teil: Bekannte Konstruktionen. Die vorliegende Arbeit gibt zunächst einen ausführlichen Überblick über die bisherigen Entwicklungen auf dem klassischen Gebiet der Hyperflächen mit vielen Singularitäten. Die maximale Anzahl mu^n(d) von Singularitäten auf einer Hyperfläche vom Grad d im P^n(C) ist nur in sehr wenigen Fällen bekannt, im P^3(C) beispielsweise nur für d<=6. Abgesehen von solchen Ausnahmen existieren nur obere und untere Schranken. 2. Teil: Neue Konstruktionen. Für kleine Grade d ist es oft möglich, bessere Resultate zu erhalten als jene, die durch allgemeine Schranken gegeben sind. In dieser Arbeit beschreiben wir einige algorithmische Ansätze hierfür, von denen einer Computer Algebra in Charakteristik 0 benutzt. Unsere anderen algorithmischen Methoden basieren auf einer Suche über endlichen Körpern. Das Liften der so experimentell gefundenen Hyperflächen durch Ausnutzung ihrer Geometrie oder Arithmetik liefert beispielsweise eine Fläche vom Grad 7 mit $99$ reellen gewöhnlichen Doppelpunkten und eine Fläche vom Grad 9 mit 226 gewöhnlichen Doppelpunkten. Diese Konstruktionen liefern die ersten unteren Schranken für mu^3(d) für ungeraden Grad d>5, die die allgemeine Schranke übertreffen. Unser Algorithmus hat außerdem das Potential, auf viele weitere Probleme der algebraischen Geometrie angewendet zu werden. Neben diesen algorithmischen Methoden beschreiben wir eine Konstruktion von Hyperflächen vom Grad d im P^n mit vielen A_j-Singularitäten, j>=2. Diese Beispiele, deren Existenz wir mit Hilfe der Theorie der Dessins d'Enfants beweisen, übertreffen die bekannten unteren Schranken in den meisten Fällen und ergeben insbesondere neue asymptotische untere Schranken für j>=2, n>=3. 3. Teil: Visualisierung. Wir beschließen unsere Arbeit mit einer Anwendung unserer neuen Visualisierungs-Software surfex, die die Stärken mehrerer existierender Programme bündelt, auf die Konstruktion affiner Gleichungen aller 45 topologischen Typen reeller kubischer Flächen.

HCHO, H 2 O 2 und CH 3 OOH in der Troposphäre

In dieser Arbeit wurden im Rahmen der UTOPIHAN- und HOHPEX04-Projekte Peroxid- und Formaldehydmessungen in der Troposphäre durchgeführt und wissenschaftlich interpretiert. Die Messungen während UTOPIHAN fanden dabei an Bord eines für Forschungszwecke umgerüsteten Flugzeuges (Learjet 35A) im Wesentlichen in der freien, insbesondere in der oberen Troposphäre über Europa statt. Die Messungen während HOHPEX04 waren hingegen als Bodenmessungen an der sich abwechselnd in der bodennahen Grenzschicht und in von dieser Schicht entkoppelten Luftmassen liegenden Bergstation Hohenpeißenberg (bayerisches Voralpenland) konzipiert. Um eine quantitative Auswertbarkeit der Messungen sicherzustellen, wurden die verwendeten, auf chemischer Derivatisierung und fluorimetrischer Detektion basierenden Messgeräte AL 2001CA (Peroxide) und AL 4021 (Formaldehyd) (AEROLASER) genau charakterisiert. Dabei wurde speziell die bekannte Ozoninterferenz beider Geräte in einer großen Zahl von Laborexperimenten mit unterschiedlichen Randbedingungen bezüglich Wasserdampf- und Kohlenwasserstoffgehalt der Luft untersucht. Für beide Verbindungen wurden Höhen- sowie Breitenprofile erstellt und mit Ergebnissen eines 3D-Chemie-Transport-Modells (CTM) sowie früherer Studien verglichen. In einem weiteren Kapitel werden Ergebnisse einer quantitativen Studie zum Einfluss hochreichender Konvektion auf das HCHO-Budget in der mittleren und oberen Troposphäre präsentiert. Diese Studie kommt zu dem Schluss, dass der rasche Aufwärtstransport von Vorläufergasen von HCHO und HOx wie Methanol, Aceton und sogar gut löslicher Spurengase wie CH3OOH beziehungsweise H2O2 aus der Grenzschicht einen signifikanten, auf Grund der längeren Lebensdauer von NOx über mehrere Tage andauernden und damit großräumigen Einfluss auf die Budgets von HCHO, HOx und auch O3 in der oberen Troposphäre haben kann. Die Befunde der Studie legen desweiteren nahe, dass fehlerhafte Modellvorhersagen für die NO-Mischungsverhältnisse in der Tropopausenregion, die zum Beispiel mit Mängeln des Modells bezüglich der Höhe der Konvektion und des Stratosphären-Troposphären-Austauschs zu tun haben, hauptverantwortlich sind für gefundene Differenzen zwischen Messdaten und dem verwendeten 3D-Chemie-Transport-Modell. Um die Signifikanz der Aussagen zu erhöhen, wurde eine Sensitivitätsstudie durchgeführt, in der die Konzentration einiger chemischer Verbindungen sowie die Photolyseraten variiert wurden. Eine weitere Studie zum Einfluss verschiedener Parameter auf das CH3OOH/H2O2-Verhältnis kommt zu dem Schluss, dass dieses Verhältnis keinen idealen Indikator für Wolkenprozessierung von Luftmassen darstellt, während eine signifikant positive Abweichung vom H2O2/H2O-Verhältnis in der oberen Troposphäre ein guter Indikator für rasch aufwärts transportierte Luftmassen sein kann. Im Rahmen dieser Studie werden auch Höhen- und Breitenprofile des CH3OOH/H2O2-Verhältnisses diskutiert. In einer letzten Untersuchung zu HCHO-Messungen am Observatorium Hohenpeißenberg im Sommer 2004 werden für die in zwei Windrichtungssektoren eingeteilten Daten Korrelationen anderer Spurengase wie O3, PAN, CO, NOy und Isopren mit HCHO interpretiert. In diesem Zusammenhang wird auch versucht, den beobachteten Tagesgang von HCHO zu erklären.

Anwendung der geoelektrischen 3D-Tomographie für die Analyse thermisch induzierter Strömungen im Labor

Sowohl in der Natur als auch in der Industrie existieren thermisch induzierte Strömungen. Von Interesse für diese Forschungsarbeit sind dabei die Konvektionen im Erdmantel sowie in den Glasschmelzwannen. Der dort stattfindende Materialtransport resultiert aus Unterschieden in der Dichte, der Temperatur und der chemischen Konzentration innerhalb des konvektierenden Materials. Um das Verständnis für die ablaufenden Prozesse zu verbessern, werden von zahlreichen Forschergruppen numerische Modellierungen durchgeführt. Die Verifikation der dafür verwendeten Algorithmen erfolgt meist über die Analyse von Laborexperimenten. Im Vordergrund dieser Forschungsarbeit steht die Entwicklung einer Methode zur Bestimmung der dreidimensionalen Temperaturverteilung für die Untersuchung von thermisch induzierten Strömungen in einem Versuchsbecken. Eine direkte Temperaturmessung im Inneren des Versuchsmaterials bzw. der Glasschmelze beeinflusst allerdings das Strömungsverhalten. Deshalb wird die geodynamisch störungsfrei arbeitende Impedanztomographie verwendet. Die Grundlage dieser Methode bildet der erweiterte Arrhenius-Zusammenhang zwischen Temperatur und spezifischer elektrischer Leitfähigkeit. Während der Laborexperimente wird ein zähflüssiges Polyethylenglykol-Wasser-Gemisch in einem Becken von unten her erhitzt. Die auf diese Weise generierten Strömungen stellen unter Berücksichtigung der Skalierung ein Analogon sowohl zu dem Erdmantel als auch zu den Schmelzwannen dar. Über mehrere Elektroden, die an den Beckenwänden installiert sind, erfolgen die geoelektrischen Messungen. Nach der sich anschließenden dreidimensionalen Inversion der elektrischen Widerstände liegt das Modell mit der Verteilung der spezifischen elektrischen Leitfähigkeit im Inneren des Versuchsbeckens vor. Diese wird mittels der erweiterten Arrhenius-Formel in eine Temperaturverteilung umgerechnet. Zum Nachweis der Eignung dieser Methode für die nichtinvasive Bestimmung der dreidimensionalen Temperaturverteilung wurden mittels mehrerer Thermoelemente an den Beckenwänden zusätzlich direkte Temperaturmessungen durchgeführt und die Werte miteinander verglichen. Im Wesentlichen sind die Innentemperaturen gut rekonstruierbar, wobei die erreichte Messgenauigkeit von der räumlichen und zeitlichen Auflösung der Gleichstromgeoelektrik abhängt.