Bei der Modellierung elektrischer Schaltungen ergeben sich Algebro-Differentialgleichungen (DAEs) mit proper formuliertem Hauptterm. Diese Gleichungen müssen z.B. bei der transienten Schaltungssimulation numerisch gelöst werden. Bei den klassischen Ansätzen der Linearen Mehrschrittverfahren oder der Runge-Kutta Verfahren ergeben sich Nachteile, die durch Verwendung von Allgemeinen Linearen Verfahren vermieden werden können. Sowohl Lineare Mehrschrittverfahren als auch Runge-Kutta Verfahren sind als Spezialfälle in dieser allgemeineren Klasse enthalten. Darüberhinaus sind aber neue Verfahren mit verbesserten Eigenschaften möglich. In dieser Arbeit werden DAEs der Schaltungssimulation eingehend studiert und Allgemeine Lineare Verfahren für solche Gleichungen untersucht. Die Verfahrenskonstruktion und Implementierungsfragen werden ausführlich diskutiert. Diese Arbeit erscheint im Logos Verlag Berlin (www.logos-verlag.de, ISBN 3-8325-1353-1). / Modelling electrical circuits leads to differential algebraic equations (DAEs) having a properly stated leading term. These equations need to be solved numerically, e.g. in case of a transient analysis of the given circuit. Classical methods such as linear multistep methods or Runge-Kutta schemes suffer from disadvantages that can be overcome by studying general linear schemes. Both Runge-Kutta methods and linear multistep schemes belong to this class as special cases, but there is plenty of room for new methods with improved properties. This work presents both a detailed study of DAEs in the framework of integrated circuit design and a thorough analysis of general linear methods for these kind of equations. The construction and implementation of general linear methods for DAEs is discussed in detail. This work is published by Logos Verlag Berlin (www.logos-verlag.de, ISBN 3-8325-1353-1).

Ziel dieser Dissertation ist die Untersuchung der asymptotischen Stabilität numerischer Verfahren für Index-2-Algebro-Differentialgleichungen. Es werden Anfangswertaufgaben für quasilineare Algebro-Differentialgleichungen (ADGln). Die meisten anwendungsrelevanten Aufgaben können damit behandelt werden. Zuerst werden einige Stabilitätsbegrife und Aussagen vorgestellt, die das Fundament für den Rest der Arbeit darstellen. Dies erstreckt sich sowohl auf den kontinuierlichen als auch auf den diskreten Fall. Insbesondere werden Kontraktivitätskonzepte eingeführt und Beziehungen zwischen der Kontraktivität der ADGl und derer der Anwendung eines numerischen Verfahrens. Die eingeführte Kontraktivitätsbegriffe erweitern oder verallgemeinern die bereits bekannten Konzepte. Als wichtigste Aussage in dem Kontraktivitätskontext geht ein Theorem hervor, das allgemeine Bedingungen aufstellt, damit die Anwendung eines IRK(DAE)-Verfahrens auf eine ADGl stabil ist. Bekannte Aussagen für gewöhnliche und Algebro-Differntialgleichungen können als Sonderfälle dieses Ergebnisses gesehen werden. Im weiteren Verlauf der Arbeit wird anhand von neuartigen Index-2-Entkopplungs- und Indexreduktionstechniken die Stabilität von Diskretisierungsverfahren untersucht. Die durchgeführte Analyse erbringt neue Ergebnisse, die eine Verbesserung des Kenntnissstandes in diesem Gebiet darstellen. Die erzielte Aussagen stellen hinreichende Bedingungen, damit ein BDF- oder IRK-Verfahren für eine ADGl das gleiche Stabilitätsverhalten wie für eine gewöhnliche Differentialgleichung besitzt. Diese Ergebnisse werden durch numerishce Beispiele veranschaulicht. Weiterhin stellt man fest, dass eine der gefundenen Voraussetzungen für die Kontraktivität der Anwendung eines algebraisch stabilen IRK(DAE)-Verfahrens, auf eine ebenfalls kontraktive ADGl, genügt. Dieses Ergebnis wurde durch die Anwendung der im ersten Teil dieser Arbeit erzielten Kontraktivitätsaussagen ermöglicht. Die Konsequenzen der soeben genannten Aussage für bestimmte Modelle der Schaltkreissimulation werden ebenfalls erläutert. Aus der oben genannten Analyse, ebenso wie aus der Fachliteratur, geht hervor, dass bei manchen ADGl-Aufgaben die Diskretisierungsverfahren Stabilitätsprobleme aufweisen. Um solche Probleme zu behandeln sind bereits einige Ansätze bekannt. Im letzten Teil der Arbeit werden zwei repräsentativen Ansätze betrachtet und ihre Aussichtschancen für Index-2-Aufgaben anhand eines kritischen Beispieles evaluiert. Des Weiteren wird eine Verallgemeinerung für vollimplizite lineare ADGln des Gear-Gupta-Leimkuhler-Ansatzes (GGL) vorgeschlagen. Der Rest der Arbeit beschäftigt sich mit der Stabilitätsuntersuchung der GGL-Formulierung und der auf sie angewandten numerischen Verfahren. Dafür werden Aussagen dieser Arbeit eingesetzt und man kommt zu der Schlussfolgerung, dass sowohl für die IRK(DAE)- als auch für die BDF-Verfahren die Integration der GGL-Formulierung, natürlich unter bestimmten Voraussetzungen, stabil ist. Dieses Ergebniss wird durch ein numerisches Beispiel belegt. Dabei handelt es um eine Gleichung, die mit einer direkten Anwendung eines Verfahrens Instabilitäten aufweist. Jedoch ist die Integration der entsprechenden GGL stabil. / The purpose of the present PhD work is the asymptotic stability investigation of numerical methods for index 2 differential algebraic equations. Initial value problems are considered for quasi linear differential algebraic equations (DAEs) that cover the most important applications. First some stability concepts and related results are presented, which represent the basis for further investigations. This background concerns both, the continuous and the discreet case. Especially contractivity concepts are introduced and the relationship between the asymptotic stability of the DAE and the numerical method applied to it is established. The new contractivity concepts extend or generalize the already known concepts. The most important result in this context is a theorem that establishes general conditions under which the application of an algebraic stable IRK(DAE) method to a DAE is contractive. Well-known assertions for ordinary and differential algebraic equations can be considered as special cases of this general result. Later on the stability of numerical discretizations applied to index-2 DAEs is investigated. This is made possible by the introduction of new decopling and index reduction techniques. The analysis makes new insights in the asymptotic of numerical methods for DAEs possible. The obtained results state sufficient conditions in order that a BDF or an IRK(DAE) method applying to DAEs shows the same asymptotic stability properties as for ODEs. These results are illustrated by some numerical examples. Moreover, it can be realized that one of the found conditions is sufficient in order to show contractivity of the application of an algebraic stable IRK(DAE) method, supposed the DAE is contractive. This assertion is possible based on the general theorem mentioned in the paragraph above. Further some consequences of the mentioned results for electric network models are shown. According to both, the above mentioned analysis and the specialized literature of this field, the application of numerical methods to some special DAEs shows asymptotic stability problems. A few approaches are known to manage such difficult equations. Two exponents of these techniques are considered and their chances of success for index-2 DAEs are evaluated with the application to a critical example. A generalization of the Gear-Gupta-Leimkuhler (GGL) approach is proposed for full implicit linear DAEs. This generalization is investigated in detail in the rest of the paper, concerning both the analytical and the numerical asymptotic stability of the GGL equation and the numerical methods applied to it correspondingly. The result is, that, if some conditions are fulfilled, IRK(DAE) and BDF methods for the GGL equation will produce stable solutions. This result is illustrated by a numerical example. The application of the methods directly to the considered DAE produces unstable solutions. However, the integration of the corresponding GGL formulation is stable. The obtained result opens new possibility for the numerical treatment of instabilities by differential algebraic equations.

Zur numerischen L\"osung von Algebro-Differentialgleichungen (ADGln) m\"ussen konsistente Anfangswerte berechnet werden. Diese Arbeit befasst sich mit einem Ansatz zur Behandlung dieses Problems f\"ur Index-2 DAEs unter Verwendung von Projektoren auf die zur DAE zugeh\"origen Unterr\"aume. Die Arbeit hat zwei Schwerpunkte.\\ Zum einen werden neue Struktureigenschaften aus schwachen Voraussetzungen hergeleitet. Anschlie{\ss}end wird eine Vorgehensweise zur Auswahl von geeigneten Gleichungen einer Index-2 ADGln vorgeschlagen, deren Differentiation zu einer Indexreduktion f\"uhrt. Diese Indexreduktion liefert neue Existenz- und Eindeutigkeitsergebnisse f\"ur L\"osungen von Index-2 ADGln. Die Ergebnisse umfassen eine allgemeinere Aufgabenklasse als die bisherigen Resultate. Beruhend auf dieser Vorgehensweise wird ein stufenweiser Ansatz zur Berechnung konsistenter Anfangswerte hergeleitet. Auf diese Weise werden neue Einsichten hinsichtlich der Ausnutzung von Struktureigenschaften von Index-2 ADGln gewonnen. Insbesondere stellt sich heraus, dass im Vergleich zu Index-1 ADGln der zus\"atzliche Schritt oft in der L\"osung eines linearen Systems besteht. Die sich hieraus ergebenden numerischen Folgen werden f\"ur zwei in der Schaltungssimulation h\"aufig verwendete Verfahren, das implizite Eulerverfahren und die Trapezregel, erl\"autert. \\ Zum anderen wird die Anwendung der erhaltenen Ergebnisse auf die Gleichungen, die bei der Schaltungssimulation mittels modifizierter Knotenanalyse entstehen, ausgearbeitet. Abschlie{\ss}end wird eine kurze \"Ubersicht der durchgef\"uhrten Umsetzung gegeben.\\ / For solving DAEs numerically, consistent initial values have to be calculated. This thesis deals with an approach for handling this problem for index-2 DAEs by considering projectors onto the spaces related to the DAE. There are two major aspects in this work.\\ On the one hand, new structural properties are deduced from weak assumptions. Subsequently, a method is proposed to choose suitable equations of an index-2 DAE, whose differentiation leads to an index reduction. This index reduction yields new theoretical results for the existence and uniqueness of solutions of index-2 DAEs which apply to a wider class of applications than previous results. Based on this method, a step-by-step approach to compute consistent initial values is developed. In this way, we gain new insights about how to deal with structural properties of index-2 DAEs. In particular, it turns out that, in comparison to index-1 DAEs, the additional step that has to be undertaken in practice often consists in solving a linear system. The numerical consequences of this fact are exemplified for two methods commonly used in circuit simulation, the implicit Euler method and the trapezoidal rule.\\ On the other hand, the application of the obtained results to the equations arising in circuit simulation by means of the modified nodal analysis (MNA) is worked out. Finally, a short overview of the specifics of their realization is given.