Numerical Solution of Differential-Algebraic Systems Arising in Circuit Simulation

In this thesis, quasi-linear Differential-Algebraic Equations (DAE) as they arise in circuit simulation are examined. The circuit is described with help of a netlist, e.g. a list which contains the information about all devices in the circuit and the way these devices are connected to each other. With the help of this netlist a DAE is created by using either the classical or the charge oriented Modified Nodal Analysis. One focus of this thesis is the graph theoretical structure of the circuit, which is reflected in the properties of the circuit DAE. For example, it is known that it is possible under certain conditions to determine the differentiation index of such a DAE based on the graph of the circuit. Moreover, it is possible to determine the derivatives of equations which yield hidden constraints. In the course of this work, two methods have been developed that use this structural information to reduce the differentiation index of a circuit DAE if the DAE has differentiation index 2. To this end the graph theoretical foundations are laid out and applied to the existing results. It is shown that because of the graph theoretical results both index reduction methods work without time consuming rank decisions. Therefore these methods are applicable to large circuits. Moreover, the second method which was developed for DAEs arising from the charge oriented Modified Nodal Analysis allows for a physical interpretation. Hence this method not only can be realized as a numerical method that alters the DAE in order to reduce the differentiation index, but also as a method that alters the circuit itself. The second approach has the advantage that existing software for circuit simulation only has to be adapted slightly to use the index reduction method. Since the charge oriented Modified Nodal Analysis is preferred in industrial circuit simulation, the academical circuit simulator Psim by Robert Melville has been adapted to use this variant to produce the circuit DAE. Following this step the differentiation index of the resulting DAE is determined by graph theoretical algorithms and the index reduction method for DAEs from charge oriented Modified Nodal Analysis as proposed in this thesis is applied. The method has been tested on various examples. A comparison with the numerical results obtained without reducing the differentiation index shows the robustness and efficiency of the method.

[1]  Anders Barrlund,et al.  Constrained least squares methods for linear timevarying DAE systems , 1991 .

[2]  F. Krogh,et al.  Solving Ordinary Differential Equations , 2019, Programming for Computations - Python.

[3]  J. Lambert.  Numerical Methods for Ordinary Differential Equations , 1991 .

[4]  C. Pantelides.  The consistent intialization of differential-algebraic systems , 1988 .

[5]  Kishore Singhal,et al.  Computer Methods for Circuit Analysis and Design , 1983 .

[6]  Helmut Werner,et al.  Gewöhnliche Differentialgleichungen , 1986, Hochschultext.

[7]  William J. Cook,et al.  Combinatorial optimization , 1997 .

[8]  James M. Ortega,et al.  Iterative solution of nonlinear equations in several variables , 2014, Computer science and applied mathematics.

[9]  Volker Mehrmann,et al.  A New Software Package for Linear Differential-Algebraic Equations , 1997, SIAM J. Sci. Comput..

[10]  Ernst Hairer,et al.  On the Instability of the BDF Formulas , 1983 .

[11]  Alexander A. Stepanov,et al.  C++ Standard Template Library , 2000 .

[12]  E. Griepentrog,et al.  Differential-algebraic equations and their numerical treatment , 1986 .

[13]  C. W. Gear,et al.  Differential-algebraic equations index transformations , 1988 .

[14]  R. März.  Differential algebraic systems anew , 2002 .

[15]  J. Butcher.  Implicit Runge-Kutta processes , 1964 .

[16]  V. Mehrmann,et al.  Mathematik Weak formulations of linear differential-algebraic systems , 2006 .

[17]  R. März.  Differential Algebraic Systems with Properly Stated Leading Term and MNA Equations , 2003 .

[18]  C. Tischendorf,et al.  Mathematical Problems in Circuit Simulation , 2001 .

[19]  Josef Stoer,et al.  Numerische Mathematik 2 , 1990 .

[20]  Stephen L. Campbell,et al.  Least squares completions for nonlinear differential algebraic equations , 1993 .

[21]  René Lamour,et al.  The computation of consistent initial values for nonlinear index-2 differential–algebraic equations , 2004, Numerical Algorithms.

[22]  Volker Mehrmann,et al.  Regular solutions of nonlinear differential-algebraic equations and their numerical determination , 1998 .

[23]  C. W. Gear,et al.  Differential-Algebraic Equations , 1984 .

[24]  Ulrich Nowak,et al.  A Family of Newton Codes for Systems of Highly Nonlinear Equations. , 1992 .

[25]  Karem A. Sakallah,et al.  A first-order charge conserving MOS capacitance model , 1990, IEEE Trans. Comput. Aided Des. Integr. Circuits Syst..

[26]  Volker Mehrmann,et al.  Stability properties of differential-algebraic equations and Spin-stabilized discretizations. , 2007 .

[27]  Béla Bollobás,et al.  Modern Graph Theory , 2002, Graduate Texts in Mathematics.

[28]  Caren Tischendorf,et al.  Structural analysis of electric circuits and consequences for MNA , 2000, Int. J. Circuit Theory Appl..

[29]  Roswitha März,et al.  Numerical methods for differential algebraic equations , 1992, Acta Numerica.

[30]  V. Mehrmann,et al.  Canonical forms for linear differential-algebraic equations with variable coefficients , 1994 .

[31]  Roswitha März,et al.  Characterizing differential algebraic equationswithout the use of derivative arrays , 2005 .

[32]  G. Dahlquist.  A special stability problem for linear multistep methods , 1963 .

[33]  V. Mehrmann,et al.  Characterization of Classes of Singular Linear Differential-Algebraic Equations , 2005 .

[34]  C. Tischendorf,et al.  Structural analysis of electric circuits and consequences for MNA , 2000 .

[35]  P.K. Chatterjee,et al.  An Investigation of the Charge Conservation Problem for MOSFET Circuit Simulation , 1983, IEEE Journal of Solid-State Circuits.

[36]  Ernst Hairer,et al.  The numerical solution of differential-algebraic systems by Runge-Kutta methods , 1989 .

[37]  E. Hairer,et al.  Solving Ordinary Differential Equations I , 1987 .

[38]  Linda R. Petzold,et al.  Numerical solution of initial-value problems in differential-algebraic equations , 1996, Classics in applied mathematics.

[39]  Volker Mehrmann,et al.  Analysis of Over- and Underdetermined Nonlinear Differential-Algebraic Systems with Application to Nonlinear Control Problems , 2001, Math. Control. Signals Syst..

[40]  Ricardo Riaza,et al.  On linear differential-algebraic equations with properly stated leading term , 2004 .

[41]  C. W. Gear,et al.  The index of general nonlinear DAEs , 1995 .

[42]  Michael Günther,et al.  CAD based electric circuit modeling in industry. Pt. 2: Impact of circuit configurations and parameters , 1997 .

[43]  H. Whitney.  On the Abstract Properties of Linear Dependence , 1935 .

[44]  R. März.  The index of linear differential algebraic equations with properly stated leading terms , 2002 .

[45]  V. Mehrmann,et al.  Index reduction for differential‐algebraic equations by minimal extension , 2004 .

[46]  R. Lamour.  Index determination and calculationof consistent initial values for DAEs , 2005 .

[47]  Sven Erik Mattsson,et al.  Index Reduction in Differential-Algebraic Equations Using Dummy Derivatives , 1993, SIAM J. Sci. Comput..

[48]  A. Schrijver.  A Course in Combinatorial Optimization , 1990 .

[49]  Stephen L. Campbell,et al.  Solvability of General Differential Algebraic Equations , 1995, SIAM J. Sci. Comput..

[50]  R.W. Dutton,et al.  A charge-oriented model for MOS transistor capacitances , 1978, IEEE Journal of Solid-State Circuits.

[51]  Bjarne Stroustrup,et al.  The C++ Programming Language: Special Edition , 2000 .

[52]  Christian Liebchen,et al.  Periodic Timetable Optimization in Public Transport , 2006, OR.

[53]  Volker Mehrmann,et al.  A New Class of Discretization Methods for the Solution of Linear Differential-Algebraic Equations with Variable Coefficients , 1996 .

[54]  D. E. Schwarz.  Topological analysis for consistent initialization in circuit simulation , 2005 .

[55]  Charles A. Desoer,et al.  Basic Circuit Theory , 1969 .

[56]  Stephen L. Campbell,et al.  A general form for solvable linear time varying singular systems of differential equations , 1987 .

[57]  Diana Estévez Schwarz.  A step-by-step approach to compute a consistent initialization for the MNA , 2002, Int. J. Circuit Theory Appl..

[58]  Paul I. Barton,et al.  Differential-Algebraic Equations of Index 1 May Have an Arbitrarily High Structural Index , 1999, SIAM J. Sci. Comput..