Search

Electric circuits are present in a number of applications, e.g. in home computers, television, credit cards, electric power networks, etc. The development of integrated circuit requires numerical simulation. Modern modeling techniques like the Modified Nodal Analysis (MNA) lead to differential algebraic equations (DAEs). Properties like the stability of solutions of such systems depend strongly on the DAE index. The paper deals with lumped circuits containing voltage sources, current sources as well as general nonlinear but time-invariant capacitances, inductances and resistances. We present network-topological criteria for the index of the DAEs obtained by the classical and the charge oriented MNA. Furthermore, the index is shown to be limited to 2 for our model-class.

The development of integrated circuits requires powerful numerical simulation programs. Of course, there is no method that treats all the different kinds of circuits successfully. The numerical simulation tools provide reliable results only if the circuit model meets the assumptions that guarantee the successful application of the integration software. Because of the large dimension of many circuits (about $10^7$ circuit elements) it is often difficult to find the circuit configurations that lead to numerical difficulties. In this paper, we analyze electric circuits with respect to their structural properties in order to give circuit designers some help for fixing modelling problems if the numerical simulation fails. We consider one of the most frequently used modelling technique, the modified nodal analysis (MNA), and discuss the index of the differential algebraic equations (DAEs) obtained by this kind of modelling.

Circuit simulation is a standard task for the computer-aided design of electronic circuits. The transient analysis is well understood and realized in powerful simulation packages for conventional circuits. But further developments in the production engineering lead to new classes of circuits which may cause difficulties for the numerical integration. The dimension of circuit models can be quite large (105 equations). The complexity of the models demands a higher abstraction level. In this paper, we analyze electric circuits with respect to their structural properties. We discuss the relevant subspaces of the resulting differential algebraic equations (DAEs) and present algorithms for calculating the index as well as consistent initial values.

For index-1 DAEs with properly stated leading term, we characterize dissipative and contractive flows and study how the qualitative properties of the DAE solutions are reflected by numerical approximations. The best situation occurs when the discretization and the decoupling procedure commute. It turns out out that this is the case if the relevant part of the inherent regular ODE has a constant state space. Different kinds of reformulations are studied to obtain numerically well formulated systems. Those reformulations might be expensive, hence, in order to avoid them, criteria ensuring the given DAE to be numerically equivalent to a numerically well formulated representation are proved.

differential algebraic systems (ADASs), i.e., differential algebraic systems with operators acting in real Hilbert spaces are introduced for a systematical treatment of coupled systems of PDEs, DAEs and integral equations. Using the finite-dimensional decoupling theory for DAEs as motivation, this paper will examine what one appropriate analogue is for infinite-dimensional systems. This leads to an index definition for ADASs. Thereby, instead of the inherent regular ODE one obtains an explicit (abstract) differential equation. In particular, when discussing PDAEs, the inherent regular differential equation is actually a parabolic PDE. The decoupling procedure provides, additionally, appropriate initial and boundary conditions for unique solvability of the coupled systems. The concept to handle ADASs is explained in different case studies.

For linear index-2 DAEs with properly stated leading term we characterize contractive and dissipative flows. We study under which conditions the qualitative properties of the DAE solutions are reflected by the numerical approximations. This is the case if the discretisation and the decoupling processes commute. Commutativity is achieved if two subspaces associated with the index-2 DAE are constant; in this case we say that the index-2 DAE is numerically well formulated. If both subspaces are time dependent, the problem should be reformulated; in order to avoid numerically equivalent reformulations, a criterion is given. If only one of these subspaces is constant, the problem is in some sense close to a numerically well formulated one and thus depending on the problem context, no reformulations are needed. Contractive and dissipative flows induced by the DAE are characterized and results on qualitative properties of the numerical solution are given.

Several qualitative properties of equilibria in electrical circuits are analyzed in this paper. Specifically, non-singularity, hyperbolicity, and asymptotic stability are addressed in terms of the circuit topology, which is captured through the use of Modified Nodal Analysis (MNA) models. The differential-algebraic or semistate nature of these models drives the analysis of the spectrum to a matrix pencil setting, and puts the results beyond the ones already known for state-space models, unfeasible in many actual problems. The topological conditions arising in this qualitative study are proved independent of those supporting the index, and therefore they apply to both index-1 and index-2 configurations. The analysis combines results coming from graph theory, matrix analysis, matrix pencil theory, and Lyapunov theory for DAEs. The study is restricted to problems with independent sources; qualitative properties of circuits including controlled sources are the focus of future research.

A model for a linear electric circuit containing semiconductors is presented. The modified nodal analysis leads to a differential algebraic equation (DAE) describing the electric network. The non-linear behaviour of the semiconductors is modelled by the drift diffusion equations. Coupling relations are defined and a sensitivity analysis concept that generalises the DAE index for finite systems to infinite ones is presented and applied to the resulting partial differential algebraic equation (PDAE). It is shown that the coupled system is of index 1 if the voltages applied to the semiconductors are low and the network without semiconductors is of index 1.

For a refined network analysis, we are interested in circuit simulation including distributed models of semiconductors. We construct a mathematical model for nonlinear electric networks containing semiconductors described by the drift-diffusion equations. The focus lies on the coupling of the network DAEs and the semiconductor PDEs. Furthermore, we study the behavior of the coupled systems with respect to time dependent perturbations using an index concept for abstract DAEs. We present a network topological criterion that guarantees index-1 systems.

In electric circuit simulation the charge oriented modified nodal analysis may lead to highly nonlinear DAEs with low smoothness properties. They may have index 2 but they do not belong to the class of Hessenberg form systems that are well understood. In the present paper, on the background of a detailed analysis of the resulting structure, it is shown that charge oriented modified nodal analysis yields the same index as the classical modified nodal analysis. Moreover, for index 2 DAEs in the charge oriented case, a further careful analysis with respect to solvability, linearization and numerical integration is given.

In this work we are interested in the numerical solution of a coupled model of differential algebraic equations (DAEs) and partial differential equations (PDEs). The DAEs describe the behavior of an electrical circuit that contains semiconductor devices and the partial differential equations constitute drift-diffusion equations modeling the semiconductor devices in the circuit. After space discretization using a finite element method, the coupled system results in a differential-algebraic system with a properly stated leading term. We investigate the structure and the properties of this DAE system. In particular, we develop structural criteria for the DAE index. This is of basic interest since DAE properties like stability, existence and uniqueness of solutions depend strongly on its index.

Die Anwenung von Operator-Splitting-Methoden auf gewöhnliche Differentialgleichungen ist gut etabliert. Für Differential-algebraische Gleichungen und partielle Differential-algebraische Gleichungen unterliegt sie jedoch vielen Einschränkungen aufgrund des Vorhandenseins von Nebenbedingungen. Die räumliche Diskretisierung reduziert PDAEs und lenkt unseren Fokus auf das Konzept der DAEs. Um eine reibungslose Übertragung des Operator-Splittings von ODEs auf DAEs durchzuführen, ist es wichtig, eine geeignete entkoppelte Struktur für das gewünschte Differential-algebraische System zu haben. In dieser Arbeit betrachten wir ein Modell, das partielle Differentialgleichungen für elektromagnetische Bauelemente - modelliert durch die Maxwell-Gleichungen - mit Differential-algebraischen Gleichungen koppelt, die die elementaren Schaltungselemente beschreiben. Nach der räumlichen Diskretisierung der klassischen Formulierung der Maxwell-Gleichungen mit Hilfe der finiten Integrationstechnik formulieren wir das resultierende gekoppelte System als Differential-algebraische Gleichung. Um eine geeignete Entkopplung zu bekommen, verwenden wir den zweigorientierten Loop-Cutset-Ansatz für die Schaltungsmodellierung. Daraus folgt, dass wir in der Lage sind, eine geeignete Operatorzerlegung so zu konstruieren, dass wir eine natürliche topologisch entkoppelte Port-Hamiltonsche DAE-Struktur erhalten. Wir schlagen einen Operator-Splitting-Ansatz für die Schaltungs-DAEs und gekoppelten Feld-Schaltungs-DAEs in entkoppelter Form vor und analysieren seine numerischen Eigenschaften. Darüber hinaus nutzen wir das Hamiltonsche Verhalten der inhärenten gewöhnlichen Differentialgleichung durch die Verwendung expliziter und energieerhaltender Zeitintegrations-methoden. Schließlich führen wir numerische Tests, um das mathematische Modell zu illustrieren und die Konvergenzergebnisse für das vorgeschlagene DAE-Operator-Splitting zu demonstrieren., The application of operator splitting methods to ordinary differential equations (ODEs) is well established. However, for differential-algebraic equations (DAEs) and partial differential-algebraic equations (PDAEs), it is subjected to many restrictions due to the presence of constraints. In constrained dynamical problems as electrical circuits or energy transport networks, DAEs arise. In order to perform a smooth transfer of the operator splitting from ODEs to DAEs, it is important to have a suitable decoupled structure for the desired differential-algebraic system. In this thesis, we consider a model which couples partial differential equations for electro- magnetic devices -modeled by Maxwell’s equations- with differential-algebraic equations describing the basic circuit elements. After spatially discretizing the classical formulation of Maxwell’s equations using the finite integration technique, we formulate the resulting coupled system as a differential-algebraic equation. To perform an appropriate decoupling, we use the branch oriented loop-cutset approach for circuit modeling. It follows that we are able to construct a suitable operator decomposition such that we obtain a natural topologically decoupled port-Hamiltonian DAE structure. We propose an operator splitting approach for the decoupled circuit and coupled field- circuit DAEs using the Lie-Trotter and Strang splitting algorithms and analyze its numerical properties. Furthermore, we exploit the Hamiltonian behavior of the system’s inherent ordinary differential equation by the utilization of explicit and energy-preserving time integration methods. Based on the convergence analysis of the ODE operator splitting method, we derive convergence results for the proposed approach that depends on the index of the system and thus on its topological structure. Finally, we perform numerical tests, to underline the mathematical model and to demonstrate the convergence results for the proposed DAE operator splitting.

Zuweilen wird gefördertes Erdgas als eine Brückentechnologie noch eine Weile erhalten bleiben, aber unsere Gasnetzinfrastruktur hat auch in einer Ära post-fossiler Brennstoffe eine Zukunft, um Klima-neutral erzeugtes Methan, Ammoniak oder Wasserstoff zu transportieren. Damit die Dispatcher der Zukunft, in einer sich fortwährend dynamisierenden Marktsituation, mit sich beständig wechselnden Kleinstanbietern, auch weiterhin einen sicheren Gasnetzbetrieb ermöglichen und garantieren können, werden sie auf moderne, schnelle Simulations- sowie performante Optimierungstechnologie angewiesen sein. Der Schlüssel dazu liegt in einem besseren Verständnis zur numerischen Behandlung nicht differenzierbarer Funktionen und diese Arbeit möchte einen Beitrag hierzu leisten. Wir werden stückweise differenzierbare Funktionen in sog. Abs-Normalen Form betrachten. Durch einen Prozess, der Abs-Linearisierung genannt wird, können wir stückweise lineare Approximationsmodelle erster Ordnung, mittels Techniken der algorithmischen Differentiation erzeugen. Jene Modelle können über Matrizen und Vektoren mittels gängiger Software-Bibliotheken der numerischen linearen Algebra auf Computersystemen ausgedrückt, gespeichert und behandelt werden. Über die Generalisierung der Formel von Faà di Bruno können auch Splinefunktionen höherer Ordnung generiert werden, was wiederum zu Annäherungsmodellen mit besserer Güte führt. Darauf aufbauend lassen sich gemischte Taylor-Kollokationsmethoden, darunter die mit Ordnung zwei konvergente generalisierte Trapezmethode, zur Integration von Gasnetzen, in Form von nicht glatten Algebro-Differentialgleichungssystemen, definieren. Numerische Experimente demonstrieren das Potential. Da solche implizite Integratoren auch nicht lineare und in unserem Falle zugleich auch stückweise differenzierbare Gleichungssysteme erzeugen, die es als Unterproblem zu lösen gilt, werden wir uns auch die stückweise differenzierbare, sowie die stückweise lineare Newtonmethode betrachten. As of yet natural gas will remain as a bridging technology, but our gas grid infrastructure does have a future in a post-fossil fuel era for the transportation of carbon-free produced methane, ammonia or hydrogen. In order for future dispatchers to continue to enable and guarantee safe gas network operations in a continuously changing market situation with constantly switching micro-suppliers, they will be dependent on modern, fast simulation as well as high-performant optimization technology. The key to such a technology resides in a better understanding of the numerical treatment of non-differentiable functions and this work aims to contribute here. We will consider piecewise differentiable functions in so-called abs-normal form. Through a process called abs-linearization, we can generate piecewise linear approximation models of order one, using techniques of algorithmic differentiation. Those models can be expressed, stored and treated numerically as matrices and vectors via common software libraries of numerical linear algebra. Generalizing the Faà di Bruno's formula yields higher order spline functions, which in turn leads to even higher order approximation models. Based on this, mixed Taylor-Collocation methods, including the generalized trapezoidal method converging with an order of two, can be defined for the integration of gas networks represented in terms of non-smooth system of differential algebraic equations. Numerical experiments will demonstrate the potential. Since those implicit integrators do generate non-linear and, in our case, piecewise differentiable systems of equations as sub-problems, it will be necessary to consider the piecewise differentiable, as well as the piecewise linear Newton method in advance.

Die Hauptthemen dieser Arbeit sind einerseits eine tiefgehende Analyse von nichtlinearen differential-algebraischen Gleichungen (DAEs) vom Index 2, die aus der modifizierten Knotenanalyse (MNA) von elektrischen Schaltkreisen hervorgehen, und andererseits die Entwicklung von Konvergenzkriterien für Waveform Relaxationsmethoden zum Lösen gekoppelter Probleme. Ein Schwerpunkt in beiden genannten Themen ist die Beziehung zwischen der Topologie eines Schaltkreises und mathematischen Eigenschaften der zugehörigen DAE. Der Analyse-Teil umfasst eine detaillierte Beschreibung einer Normalform für Schaltkreis DAEs vom Index 2 und Abschätzungen, die für die Sensitivität des Schaltkreises bezüglich seiner Input-Quellen folgen. Es wird gezeigt, wie diese Abschätzungen wesentlich von der topologischen Position der Input-Quellen im Schaltkreis abhängen. Die zunehmend komplexen Schaltkreise in technologischen Geräten erfordern oftmals eine Modellierung als gekoppeltes System. Waveform relaxation (WR) empfiehlt sich zur Lösung solch gekoppelter Probleme, da sie auf die Subprobleme angepasste Lösungsmethoden und Schrittweiten ermöglicht. Es ist bekannt, dass WR zwar bei Anwendung auf gewöhnliche Differentialgleichungen konvergiert, falls diese eine Lipschitz-Bedingung erfüllen, selbiges jedoch bei DAEs nicht ohne Hinzunahme eines Kontraktivitätskriteriums sichergestellt werden kann. Wir beschreiben allgemeine Konvergenzkriterien für WR auf DAEs vom Index 2. Für den Fall von Schaltkreisen, die entweder mit anderen Schaltkreisen oder mit elektromagnetischen Feldern verkoppelt sind, leiten wir außerdem hinreichende topologische Konvergenzkriterien her, die anhand von Beispielen veranschaulicht werden. Weiterhin werden die Konvergenzraten des Jacobi WR Verfahrens und des Gauss-Seidel WR Verfahrens verglichen. Simulationen von einfachen Beispielsystemen zeigen drastische Unterschiede des WR-Konvergenzverhaltens, abhängig davon, ob die Konvergenzbedingungen erfüllt sind oder nicht. The main topics of this thesis are firstly a thorough analysis of nonlinear differential-algebraic equations (DAEs) of index 2 which arise from the modified nodal analysis (MNA) for electrical circuits and secondly the derivation of convergence criteria for waveform relaxation (WR) methods on coupled problems. In both topics, a particular focus is put on the relations between a circuit's topology and the mathematical properties of the corresponding DAE. The analysis encompasses a detailed description of a normal form for circuit DAEs of index 2 and consequences for the sensitivity of the circuit with respect to its input source terms. More precisely, we provide bounds which describe how strongly changes in the input sources of the circuit affect its behaviour. Crucial constants in these bounds are determined in terms of the topological position of the input sources in the circuit. The increasingly complex electrical circuits in technological devices often call for coupled systems modelling. Allowing for each subsystem to be solved by dedicated numerical solvers and time scales, WR is an adequate method in this setting. It is well-known that while WR converges on ordinary differential equations if a Lipschitz condition is satisfied, an additional convergence criterion is required to guarantee convergence on DAEs. We present general convergence criteria for WR on higher index DAEs. Furthermore, based on our results of the analysis part, we derive topological convergence criteria for coupled circuit/circuit problems and field/circuit problems. Examples illustrate how to practically check if the criteria are satisfied. If a sufficient convergence criterion holds, we specify at which rate of convergence the Jacobi and Gauss-Seidel WR methods converge. Simulations of simple benchmark systems illustrate the drastically different convergence behaviour of WR depending on whether or not the circuit topological convergence conditions are satisfied.

Diese Arbeit beschäftigt sich mit verschiedenen mathematischen Fragestellungen hinsichtlich der Modellierung, Analysis und numerischen Simulation von Gasnetzen. Hierbei liegt der Fokus auf der mathematischen Handhabung von partiellen differential-algebraischen Gleichungen, die mit algebraischen Gleichungen gekoppelt sind. Diese bieten einen einfachen Zugang hinsichtlich der Modellierung von dynamischen Strukturen auf Netzen Somit sind sie insbesondere für Gasnetze geeignet, denen im Zuge der steigenden Bedeutung von erneuerbaren Energien ein gestiegenes Interesse seitens der Öffentlichkeit, Politik und Wissenschaft entgegen gebracht wird. Wir führen zunächst die gängigsten Elemente, die in Gasnetzen benötigt werden ein und formulieren zwei PDAE-Klassen für solche Netze: Eine für reine Rohrnetze, und eine, die zusätzliche Elemente wie Verdichter und Widerstände beinhaltet. Des Weiteren untersuchen wir die Sensitivität der Lösung der Rohrnetz-PDAE hinsichtlich Störungen. Dabei berücksichtigen wir Störungen, die nicht nur den dynamischen Teil der PDAE beeinflussen, sondern auch Störungen in den algebraischen Gleichungen und weisen Stabilitätseigenschaften für die Lösung der PDAE nach. Darüber hinaus beschäftigen wir uns mit einer neu entwickelten, an die Netztopologie angepassten Ortsdiskretisierung, welche die Stabilitätseigenschaften der PDAE auf DAE Systeme überträgt. Des Weiteren zeigen wir, wie sich die Gasnetz-DAE zu einer gewöhnlichen Differentialgleichung, welche die inhärente Dynamik der DAE widerspiegelt entkoppeln lässt. Dieses entkoppelte System kann darüber hinaus direkt aus den Topologie- und Elementinformationen des Netzes aufgestellt werden. Abschließend demonstrieren wir die Ergebnisse an Benchmark-Gasnetzen. Dabei vergleichen wir sowohl die entkoppelte Differentialgleichung mit dem ursprünglichen DAE System, zeigen aber auch, welche Vorteile die an die Netztopologie angepasste Ortsdiskretisierung gegenüber existierenden Verfahren besitzt. This thesis addresses several aspects regarding modelling, analysis and numerical simulation of gas networks. Hereby, our focus lies on (partial) differential-algebraic equations, thus systems of partial and ordinary differential equations which are coupled by algebraic equations. These coupled systems allow an easy approach towards the modelling of dynamic structures on networks. Therefore, they are well suited for gas networks, which have gained a rise of attention in society, politics and science due to the focus towards renewable energies. We give an introduction towards gas network modelling that includes the most common elements that also appear in real gas networks and present two PDAE systems: One for pipe networks and one that includes additional elements like resistors and compressors. Furthermore, we investigate the impact of perturbations onto the pipe network PDAE, where we explicitly allow perturbations to affect the system in the differential as well as in the algebraic components. We conclude that the solution of the PDAE possesses stability properties. In addition, this thesis introduces a new spatial discretisation that is adapted to the net- work topology. This topology-adapted semi-discretisation results in a DAE which possesses the same perturbation behaviour as the space continuous PDAE. Furthermore, we present a topology based decoupling procedure that allows to reformulate the DAE as an ordinary differential equation (ODE), which represents the inherent dynamics of the DAE system. This ODE, together with a decoupled set of algebraic equations, can be derived from the topology and element information directly. We conclude by demonstrating the established results for several benchmark networks. This includes a comparison of numerical solutions for the decoupled ODE and the DAE system. In addition we present the advantages of the topology-adapted spatial discretisation over existing well established methods.

Dynamik auf Netzwerken ist ein mathematisches Feld, das in den letzten Jahrzehnten schnell gewachsen ist und Anwendungen in zahlreichen Disziplinen wie z.B. Physik, Biologie und Soziologie findet. Die Funktion vieler Netzwerke hängt von der Fähigkeit ab, die Elemente des Netzwerkes zu synchronisieren. Mit anderen Worten, die Existenz und die transversale Stabilität der synchronen Mannigfaltigkeit sind zentrale Eigenschaften. Erst seit einigen Jahren wird versucht, den verwickelten Zusammenhang zwischen der Kopplungsstruktur und den Stabilitätseigenschaften synchroner Zustände zu verstehen. Genau das ist das zentrale Thema dieser Arbeit. Zunächst präsentiere ich erste Ergebnisse zur Klassifizierung der Kanten eines gerichteten Netzwerks bezüglich ihrer Bedeutung für die Stabilität des synchronen Zustands. Folgend untersuche ich ein komplexes Verzweigungsszenario in einem gerichteten Ring von Stuart-Landau Oszillatoren und zeige, dass das Szenario persistent ist, wenn dem Netzwerk eine schwach gewichtete Kante hinzugefügt wird. Daraufhin untersuche ich synchrone Zustände in Ringen von Phasenoszillatoren die mit Zeitverzögerung gekoppelt sind. Ich bespreche die Koexistenz synchroner Lösungen und analysiere deren Stabilität und Verzweigungen. Weiter zeige ich, dass eine Zeitverschiebung genutzt werden kann, um Muster im Ring zu speichern und wiederzuerkennen. Diese Zeitverschiebung untersuche ich daraufhin für beliebige Kopplungsstrukturen. Ich zeige, dass invariante Mannigfaltigkeiten des Flusses sowie ihre Stabilität unter der Zeitverschiebung erhalten bleiben. Darüber hinaus bestimme ich die minimale Anzahl von Zeitverzögerungen, die gebraucht werden, um das System äquivalent zu beschreiben. Schließlich untersuche ich das auffällige Phänomen eines nichtstetigen Übergangs zu Synchronizität in Klassen großer Zufallsnetzwerke indem ich einen kürzlich eingeführten Zugang zur Beschreibung großer Zufallsnetzwerke auf den Fall zeitverzögerter Kopplungen verallgemeinere. Since a couple of decades, dynamics on networks is a rapidly growing branch of mathematics with applications in various disciplines such as physics, biology or sociology. The functioning of many networks heavily relies on the ability to synchronize the network’s nodes. More precisely, the existence and the transverse stability of the synchronous manifold are essential properties. It was only in the last few years that people tried to understand the entangled relation between the coupling structure of a network, given by a (di-)graph, and the stability properties of synchronous states. This is the central theme of this dissertation. I first present results towards a classification of the links in a directed, diffusive network according to their impact on the stability of synchronization. Then I investigate a complex bifurcation scenario observed in a directed ring of Stuart-Landau oscillators. I show that under the addition of a single weak link, this scenario is persistent. Subsequently, I investigate synchronous patterns in a directed ring of phase oscillators coupled with time delay. I discuss the coexistence of multiple of synchronous solutions and investigate their stability and bifurcations. I apply these results by showing that a certain time-shift transformation can be used in order to employ the ring as a pattern recognition device. Next, I investigate the same time-shift transformation for arbitrary coupling structures in a very general setting. I show that invariant manifolds of the flow together with their stability properties are conserved under the time-shift transformation. Furthermore, I determine the minimal number of delays needed to equivalently describe the system’s dynamics. Finally, I investigate a peculiar phenomenon of non-continuous transition to synchrony observed in certain classes of large random networks, generalizing a recently introduced approach for the description of large random networks to the case of delayed couplings.

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