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.

Finite-element discretization of 3D energy-transport equations for semiconductors

In this thesis a mathematical model was derived that describes the charge and energy transport in semiconductor devices like transistors. Moreover, numerical simulations of these physical processes are performed. In order to accomplish this, methods of theoretical physics, functional analysis, numerical mathematics and computer programming are applied. After an introduction to the status quo of semiconductor device simulation methods and a brief review of historical facts up to now, the attention is shifted to the construction of a model, which serves as the basis of the subsequent derivations in the thesis. Thereby the starting point is an important equation of the theory of dilute gases. From this equation the model equations are derived and specified by means of a series expansion method. This is done in a multi-stage derivation process, which is mainly taken from a scientific paper and which does not constitute the focus of this thesis. In the following phase we specify the mathematical setting and make precise the model assumptions. Thereby we make use of methods of functional analysis. Since the equations we deal with are coupled, we are concerned with a nonstandard problem. In contrary, the theory of scalar elliptic equations is established meanwhile. Subsequently, we are preoccupied with the numerical discretization of the equations. A special finite-element method is used for the discretization. This special approach has to be done in order to make the numerical results appropriate for practical application. By a series of transformations from the discrete model we derive a system of algebraic equations that are eligible for numerical evaluation. Using self-made computer programs we solve the equations to get approximate solutions. These programs are based on new and specialized iteration procedures that are developed and thoroughly tested within the frame of this research work. Due to their importance and their novel status, they are explained and demonstrated in detail. We compare these new iterations with a standard method that is complemented by a feature to fit in the current context. A further innovation is the computation of solutions in three-dimensional domains, which are still rare. Special attention is paid to applicability of the 3D simulation tools. The programs are designed to have justifiable working complexity. The simulation results of some models of contemporary semiconductor devices are shown and detailed comments on the results are given. Eventually, we make a prospect on future development and enhancements of the models and of the algorithms that we used.

Polymer properties measured by micromechanical cantilever technique

Surface stress changes induced by specific adsorption of molecules were investigated using a micromechanical cantilever sensor (MCS) device. 16 MCS are grouped within four separate wells. Each well can be addressed independently by different liquid enabling functionalization of MCS separately by flowing different solutions through each well and performing sensing and reference experiments simultaneously. In addition, each well contains a fixed reference mirror, which allows measuring the absolute bending of MCS. The effect of the flow rate on the MCS bending change was found to be dependent on the absolute bending value of MCS. In addition, the signal from the reference mirror can be used to follow refractive index changes upon mixing different solutions. Finite element simulation of solution exchange in wells was compared with experiment results. Both revealed that one solution can be exchanged by another one after a total volume of 200 µl has flown through. Using MCS, the adsorption of thiolated deoxyribonucleic acid (DNA) molecules and 6-mercapto-1-hexanol (MCH) on gold surfaces, and the DNA hybridization were performed. The nanomechanical response is in agreement with data reported by Fritz et al.1 Thus, the multiwell device is readily applicable for sensing of multiple chemical and biological recognition events in a single step. In this context controlled release and uptake of drugs are currently widely discussed. As a model system, we have used polystyrene (PS) spheres with diameters in the order of µm. The swelling behavior of individual PS spheres in toluene vapor was studied via mass loading by means of micromechanical cantilever sensors. For 4–8% cross-linked PS a mass increase of 180% in saturated toluene vapor was measured. In addition, the diameter change in saturated toluene vapor was measured and the corresponding volume increase of 200% was calculated. The mass of the swollen PS sphere decreases with increasing exposure time to ultraviolet (UV) light. The swelling response is significantly different between the first and the second exposure to toluene vapor. This is attributed to the formation of a cross-linked shell at the surface of the PS spheres. Shape persistent parts were observed for locally UV irradiated PS spheres. These PS spheres were found to be fluorescent and cracks occur after exposure in toluene liquid. The diffusion time of dye molecules in PS spheres increases with increasing chemical cross-linking density. This concept of locally dissolving non cross-linked PS from the sphere was applied to fabricate donut structures on surfaces. Arrays of PS spheres were fabricated using spin coating. The donut structure was produced simply after liquid solvent rinsing. The complete cross-linking of PS spheres was found after long exposure time to UV. We found that stabilizers play a major role in the formation of the donut nanostructures.

Wet bioadhesion in tree frogs

Viele Tiere wie etwa Geckos oder Laubfrösche können mittels ihrer Haftscheiben an Oberflächen kleben. Diese Haftscheiben ermöglichen es den Tieren, sich während ihrerrnFortbewegung an Oberflächen anzuheften und wieder zu lösen unabhängig von denrnvorherrschenden Umweltbedingungen. Frösche besitzen mikro- und nanostrukturierternsowie charakteristisch geformte Haftscheiben an Finger- und Zehenenden. Ihre besonderernevolutionäre Errungenschaft, sich stark und zugleich reversibel in sowohl trockenen alsrnauch feuchten Umgebungen anzuhaften, hat die Wissenschaft zur Nachahmung und Untersuchungrndieser Strukturen inspiriert. Zum besseren Verständnis der Mechanismen vonrnAnhaftung und Loslösung bei Laubfröschen wurden weiche, elastische und mikrostrukturierternOberflächen hergestellt, indem PDMS (Polydimethylsiloxan) auf einer Siliziummaskernmit Hexagonstruktur aufgetragen und vernetzt wurde. Dadurch wurden Anordnungenrnvon hexagonalen Mikrosäulen mit spezifischen geometrischen Eigenschaften undrnunterschiedlichen Kontaktgeometrien (normale, flache Form, T-Form und konkave Formrnder Säulenenden) erhalten. Um den Einfluss der van-der-Waals, hydrodynamischen,rnKapillar-und Adhäsionskräfte zu verstehen, wurden verschiedene experimentelle Ansätzernverfolgt: Die auf eine einzelne Säule wirkenden Adhäsionskräfte wurden mittelsrnRasterkraftmikroskopie gemessen. Dazu wurden speziell hergestellte kolloidale Sensorenrnverwendet. Diese Experimente wurden sowohl mit als auch ohne Flüssigkeitsfilm auf derrnSäule durchgeführt. Die Ergebnisse zeigten den Beitrag von Kapillarkraft und direktenrnKontaktkräften zur Adhäsionskraft bei Vorliegen eines Flüssigkeitsfilms. Die Adhäsionrnfiel umso größer aus, je weniger Flüssigkeit zwischen Sensor und Säule vorhanden war.rnIm Falle einer trockenen Adhäsion zeigte die Säule mit T-Form die höchste Adhäsion. Darndie Haftscheiben der Laubfrösche weich sind, können sie dynamisch ihre Form ändern,rnwas zu einer Änderung der hydrodynamischen Kraft zwischen Scheibe und Oberflächernführt. Der Einfluss der Oberflächenverformbarkeit auf die hydrodynamische Kraft wurderndaher am Modellsystem einer Kugel untersucht, welche sich einer weichen und ebenenrnOberfläche annähert. Dieses System wurde sowohl theoretisch über die Simulation finiterrnElemente als auch experimentell über die Messung mit kolloidalen Sonden untersucht.rnSowohl experimentelle Ergebnisse als auch die Simulationen ergaben eine Abnahme derrnhydrodynamischen Kraft bei Annäherung des kolloidalen Sensors an eine weiche undrnelastische Oberfläche. Beim Entfernen der Sensors von der Oberfläche verstärkte sichrndie hydrodynamische Anziehungskraft. Die Kraft, die zur Trennung eines Partikels von einer Oberfläche in Flüssigkeit notwendig ist, ist für weiche und elastischen Oberflächenrngrößer als für harte Oberflächen. In Bezug zur Bioadhäsion bei Laubfröschen konnternfestgestellt somit festgestellt werden, dass sich der hydrodynamische Anteil zur feuchtenrnBioadhäsion aufgrund der weichen Oberfläche erhöht. Weiterhin wurde der Einflussrndes Aspektverhältnisses der Säulen auf die Reibungskraft mittels eines kolloidalen Sensorsrnuntersucht. Gestreckte Säulen zeigten dabei eine höhere Reibung im Vergleich zu.rnSäulen mit einem gestreckten Hexagon als Querschnitt.

Analysis and numerical solution of the Peterlin viscoelastic model

Liquids and gasses form a vital part of nature. Many of these are complex fluids with non-Newtonian behaviour. We introduce a mathematical model describing the unsteady motion of an incompressible polymeric fluid. Each polymer molecule is treated as two beads connected by a spring. For the nonlinear spring force it is not possible to obtain a closed system of equations, unless we approximate the force law. The Peterlin approximation replaces the length of the spring by the length of the average spring. Consequently, the macroscopic dumbbell-based model for dilute polymer solutions is obtained. The model consists of the conservation of mass and momentum and time evolution of the symmetric positive definite conformation tensor, where the diffusive effects are taken into account. In two space dimensions we prove global in time existence of weak solutions. Assuming more regular data we show higher regularity and consequently uniqueness of the weak solution. For the Oseen-type Peterlin model we propose a linear pressure-stabilized characteristics finite element scheme. We derive the corresponding error estimates and we prove, for linear finite elements, the optimal first order accuracy. Theoretical error of the pressure-stabilized characteristic finite element scheme is confirmed by a series of numerical experiments.