In chapter 2, we study two Kaup-Newell-type matrix spectral problems, derive their soliton hierarchies within the zero curvature formulation, and furnish their bi-Hamiltonian structures by the trace identity to show that they are integrable in the Liouville sense. In chapter 5, we obtain the Riemann theta function representation of solutions for the first hierarchy of generalized Kaup-Newell systems. In chapter 3, using Hirota bilinear forms, we discuss positive quadratic polynomial solutions to generalized bilinear equations, which generate lump or lump-type solutions to nonlinear evolution equations, and propose an algorithm for computing higher-order lump or lump-type solutions. In chapter 4, we study mixed exponential and trigonometric wave solutions (called complexitons) to general bilinear equations, and propose two methods to find complexitons to generalized bilinear equations. We also succeed in proving that by choosing suitable complex coefficients in soliton solutions, multi-complexitons are actually real wave solutions from complex soliton solutions and establish the linear superposition principle for complexion solutions. In each chapter, we present computational examples.

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).

Algebro-Differentialgleichungen (engl. differential-algebraic equations - DAEs) sind implizite singuläre gewöhnliche Differentialgleichungen, die restringierte dynamische Prozesse beschreiben. Sie unterscheiden sich von expliziten gewöhnlichen Differentialgleichungen dahingehend, dass Anfangswerte nicht beliebig vorgegeben werden können. Weiterhin sind in einer DAE neben Integrations- auch Differentiationsaufgaben involviert. Der Differentiationsindex gibt an, wieviele Differentiationen zur Lösung notwendig sind. Seit den 1980er Jahren wird vorwiegend an der Charakterisierung und Klassifizierung regulärer DAEs und der Konstruktion nebst Fundierung von Integrationsmethoden gearbeitet. I. Higueras, R. März und C. Tischendorf haben gezeigt, dass man lineare DAEs mit properem Hauptterm, A(t)(D(t)x(t))'' + B(t)x(t) = q(t), die regulär mit Traktabilitätsindex 2 sind, zuverlässig numerisch integrieren kann - im Unterschied zu linearen DAEs in Standardform. In Publikationen von R. Riaza und R. März wird die Klassifizierungen kritischer Punkten von linearen DAEs an die Verletzung bestimmter Rangbedingungen von Matrixfunktionen im Rahmen des Traktabilitätsindexes geknüpft. Im wesentlichen heißt ein kritischer Punkt harmlos, wenn der durch die inhärente Differentialgleichung beschriebene Fluß nicht tangiert ist. Gegenstand der vorliegenden Arbeit sind lineare quasi-proper formulierte DAEs. Es werden Index 2 DAEs mit harmlosen kritischen Punkten charakterisiert. Unter Verwendung von quasi-zulässigen Projektorfunktionen können neben DAEs, die fast überall gleiche charakteristische Werte haben, nun erstmalig auch solche mit Indexwechseln behandelt werden. Der Hauptteil der Arbeit besteht im Nachweis von Durchführbarkeit, Konvergenz und nur schwacher Instabilität von numerischen Integrationsmethoden (BDF, IRK(DAE)) für lineare Index 2 DAEs mit harmlosen kritischen Punkten, sowie in der Entwicklung von Fehlerschätzern und Schrittweitensteuerung. / Differential-algebraic equations (DAEs) are implicit singular ordinary differential equations, which describe dynamical processes that are restricted by some constraints. In contrast to explicit regular ordinary differential equations, for a DAE not any value can be imposed as an initial condition. Furthermore, DAEs involve not only integration problems but also differentiation problems. The differentiation index of a DAE indicates the number of differentiations required in order to solve a DAE. Since the 1980th, research focuses primarily on the characterization and classification of regular problem classes and the construction and foundation of integration methods for simulation software. I. Higueras, R. Maerz, and C. Tischendorf have shown that one can reliably integrate a general linear DAE with a properly stated leading term, A(t)(D(t)x(t))'' + B(t)x(t) = q(t), which is regular with tractability index 2 - in contrast to linear standard form DAEs. The first classification of critical points of linear DAEs has been published by R. Riaza and R. Maerz. Based on the tractability index, critical points are classified according to failures of certain rank conditions of matrix functions. Essentially, a critical point is said to be harmless, if the flow described by the inherent differential equation is not affected. The subject of this work are quasi-proper linear DAEs. Index-2 DAEs with harmless critical points are characterized. Under the application of quasi-admissible projector functions. Besides DAEs which have almost everywhere the same characteristic values, DAEs with index changes can now be discussed for the first time. The main part of the work is to provide a proof of feasibility, convergence, and only weak instability of numerical integration methods (BDF, IRK (DAE)) for linear index-2 DAEs with harmless critical points, as well as the development and testing of error estimators and stepsize control.

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.