A Structure Preserving Index Reduction Method For MNA

Transient analysis in industrial chip design leads to very large systems of differential-algebraic equations (DAEs). The numerical solution of these DAEs strongly depends on the so called index of the DAE. In general, the higher the index of the DAE is, the more sensitive the numerical solution will be to errors in the computation. So, it is advisable to use mathematical models with small index or to reduce the index. This paper presents an index reduction method that uses information based on the topology of the circuit. In addition, we show that the presented method retains structural properties of the DAE. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The index of general nonlinear DAEs

Summary. In the last few years there has been considerable research on differential algebraic equations (DAEs) where is identically singular. Much of the mathematical effort has focused on computing a solution that is assumed to exist. More recently there has been some discussion of solvability of DAEs. There has historically been some imprecision in the use of the two key concepts of solvability and index for DAEs. The index is also important in control and systems theory but with different terminology. The consideration of increasingly complex nonlinear DAEs makes a clear and correct development necessary. This paper will try to clarify several points concerning the index. After establishing some new and more precise terminology that we need, some inaccuracies in the literature will be corrected. The two types of indices most frequently used, the differentiation index and the perturbation index, are defined with respect to solutions of unperturbed problems. Examples are given to show that these indices can be very different for the same problem. We define new "maximum indices," which are the maxima of earlier indices in a neighborhood of the solution over a set of perturbations and show that these indices are simply related to each other. These indices are also related to an index defined in terms of Jacobians. Received November 15, 1993 / Revised version received December 23, 1994     

Half-explicit Runge-Kutta methods for semi-explicit differential-algebraic equations of index 1

Summary For the numerical solution of non-stiff semi-explicit differentialalgebraic equations (DAEs) of index 1 half-explicit Runge-Kutta methods (HERK) are considered that combine an explicit Runge-Kutta method for the differential part with a simplified Newton method for the (approximate) solution of the algebraic part of the DAE. Two principles for the choice of the initial guesses and the number of Newton steps at each stage are given that allow to construct HERK of the same order as the underlying explicit Runge-Kutta method. Numerical tests illustrate the efficiency of these methods.     

Index Reduction and Discontinuity Handling Using Substitute Equations

Several techniques exist for index reduction and consistent initialization of higher index DAEs. Many such techniques change the original set of equations by differentiation, substitution, and/or introduction of new variables. This paper introduces substitute equations as a new language element. By means of a substitute equation, the value of a continuous variable or its time derivative is specified by an expression. This expression is evaluated each time that the variable or its time derivative, respectively, is referenced in the model. The advantage of substitute equations is that they enable index reduction and consistent initialization of higher index DAEs without changing the original equations; no existing variables are removed and no new variables are introduced. Substitute equations can also be used to enable the use of general purpose numerical solvers for equations where one or more of the unknowns are discontinuous, and they can be used to prevent functions to be called outside of their domain.     

Element-based index reduction in electrical circuit simulation

The numerical simulation of very large scale integrated circuits is an important tool in the development of new industrial circuits. In the course of the last years, this topic has received increasing attention. Common modeling approaches for circuits lead to differential-algebraic systems (DAEs). In circuit simulation, these DAEs are known to have index 2, given some topological properties of the network. This higher index leads to several undesirable effects in the numerical solution of the DAEs. Recent approaches try to lower the index of DAEs to improve the numerical behaviour. These methods usually involve costly algebraic transformations of the equations. Especially, for large scale circuit equations, these transformations become too costly to be efficient. We will present methods that change the topology of the network itself, while replacing certain elements in oder to obtain a network that leads to a DAE of index 1, while not altering the analytical solution of the DAE. This procedure can be performed prior to the actual numerical simulation. The decreasing of the index usually leads to significantly improved numerical behaviour. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Index of an Infinite Dimensional Implicit System

The idea of the index of a differential algebraic equation (DAE) (or implicit differential equation) has played a fundamental role in both the analysis of DAEs and the development of numerical algorithms for DAEs. DAEs frequently arise as partial discretizations of partial differential equations (PDEs). In order to relate properties of the PDE to those of the resulting DAE it is necessary to have a concept of the index of a possibly constrained PDE. Using the finite dimensional theory as motivation, this paper will examine what one appropriate analogue is for infinite dimensional systems. A general definition approach will be given motivated by the desire to consider numerical methods. Specific examples illustrating several kinds of behavior will be considered in some detail. It is seen that our definition differs from purely algebraic definitions. Numerical solutions, and simulation difficulties, can be misinterpreted if this index information is missing.     

Non-stiff integrators for differential–algebraic systems of index 2

Non-stiff differential-algebraic equations (DAEs) can be solved efficiently by partitioned methods that combine well-known non-stiff integrators from ODE theory with an implicit method to handle the algebraic part of the system. In the present paper we consider partitioned one-step and partitioned multi-step methods for index-2 DAEs in Hessenberg form and the application of these methods to constrained mechanical systems. The methods are presented from a unified point of view. The comparison of various classes of methods is completed by numerical tests for benchmark problems from the literature. This revised version was published online in June 2006 with corrections to the Cover Date.     

Continuous-time accelerated block successive overrelaxation methods for time-dependent Stokes equations

In order to solve the time-dependent Stokes equation, we follow the “Method of Lines” to obtain structured linear constant-coefficient differential–algebraic equations (DAEs). By taking advantage of the structure, we propose a class of waveform relaxation methods, called continuous-time accelerated block SOR (CABSOR) methods, for solving the obtained DAEs. The new methods are theoretically analyzed. The theory is applied to a two-dimensional time-dependent Stokes equation and verified by numerical experiments.     

Index minimization of differential-algebraic equations in hybrid analysis for circuit simulation

Modern modeling approaches for circuit analysis lead to differential-algebraic equations (DAEs). The index of a DAE is a measure of the degree of numerical difficulty. In general, the higher the index is, the more difficult it is to solve the DAE. The index of the DAE arising from the modified nodal analysis (MNA) is determined uniquely by the structure of the circuit. Instead, we consider a broader class of analysis method called the hybrid analysis. For linear time-invariant electric circuits, we devise a combinatorial algorithm for finding an optimal hybrid analysis in which the index of the DAE to be solved attains the minimum. The optimal hybrid analysis often results in a DAE with lower index than MNA.      

Optimal Control of Differential–Algebraic Equations of Higher Index, Part 1: First-Order Approximations

This paper deals with optimal control problems described by higher index DAEs. We introduce a class of these problems which can be transformed to index one control problems. For this class of higher index DAEs, we derive first-order approximations and adjoint equations for the functionals defining the problem. These adjoint equations are then used to state, in the accompanying paper, the necessary optimality conditions in the form of a weak maximum principle. The constructive way used to prove these optimality conditions leads to globally convergent algorithms for control problems with state constraints and defined by higher index DAEs.     

Boundary value methods for the solution of differential-algebraic equations

Summary Boundary value techniques for the solution of initial value problems of ODEs, despite their apparent higher cost, present some important advantages over initial value methods. Among them, there is the possibility to have greater accuracy, to control the global error, and to have an efficient parallel implementation.In this paper, the same techniques are applied to the solution of linear initial value problems of DAEs. We have considered three term numerical methods (Midpoint, Simpson, and an Adams type method) in order to obtain a block tridiagonal linear system as a discrete problem.Convergence results are stated in the case of constant coefficients, and numerical examples are given on linear time-varying problems.Work supported by the Ministero della Ricerca Scientifica, 40% project, and by the C.N.R. (contract of research # 92.00535.CT01)     

A new approach to the incorporation of servo constraints in multibody dynamics

A new index reduction approach is developed to solve the servo constraint problems [2] in the inverse dynamics simulation of underactuated mechanical systems. The servo constraint problem of underactuated systems is governed by differential algebraic equations (DAEs) with high index. The underlying equations of motion contain both holonomic constraints and servo constraints in which desired outputs (specified in time) are described in terms of state variables. The realization of servo constraints with the use of control forces can range from orthogonal to tangential [3]. Since the (differentiation) index of the DAEs is often higher than three for underactuated systems, in which the number of degrees of freedom is greater than the control outputs/inputs, we propose a new index reduction method [1] which makes possible the stable numerical integration of the DAEs. We apply the proposed method to differentially flat systems, such as cranes [1,4,5], and non-flat underactuated systems. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The computation of consistent initial values for nonlinear index-2 differential–algebraic equations

The computation of consistent initial values for differential–algebraic equations (DAEs) is essential for starting a numerical integration. Based on the tractability index concept a method is proposed to filter those equations of a system of index-2 DAEs, whose differentiation leads to an index reduction. The considered equation class covers Hessenberg-systems and the equations arising from the simulation of electrical networks by means of Modified Nodal Analysis (MNA). The index reduction provides a method for the computation of the consistent initial values. The realized algorithm is described and illustrated by examples.     

Beyond conventional Runge–Kutta methods in numerical integration of ODEs and DAEs by use of structures and local models

There are two parts in this paper. In the first part we consider an overdetermined system of differential-algebraic equations (DAEs). We are particularly concerned with Hamiltonian and Lagrangian systems with holonomic constraints. The main motivation is in finding methods based on Gauss coefficients, preserving not only the constraints, symmetry, symplecticness, and variational nature of trajectories of holonomically constrained Hamiltonian and Lagrangian systems, but also having optimal order of convergence. The new class of (s,s)$(s,s)$-Gauss–Lobatto specialized partitioned additive Runge–Kutta (SPARK) methods uses greatly the structure of the DAEs and possesses all desired properties. In the second part we propose a unified approach for the solution of ordinary differential equations (ODEs) mixing analytical solutions and numerical approximations. The basic idea is to consider local models which can be solved efficiently, for example analytically, and to incorporate their solution into a global procedure based on standard numerical integration methods for the correction. In order to preserve also symmetry we define the new class of symmetrized Runge–Kutta methods with local model (SRKLM).     

Variational splines for the numerical solution of singular differential equations

The numerical parametrization method (PM), originally created for optimal control problems, is specificated for classical calculus of variation problems that arise in connection with singular implicit (IDEs) and differential-algebraic equations (DAEs). The PM for IDEs is based on representation of the required solution as a spline with moving knots and on minimization of the discrepancy functional with respect to the spline parameters. Such splines are named variational splines. For DAEs only finite entering functions can be represented by splines, and the functional under minimization is the discrepancy of the algebraic subsystem. The first and the second derivatives of the functionals are calculated in two ways – for DAEs with the help of adjoint variables, and for IDE directly. The PM does not use the notion of differentiation index, and it is applicable to any singular equation having a solution. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Enumeration results for Runge–Kutta methods for index 2 DAEs

Recursive formulas are presented that give the number of order conditions for single-step Runge–Kutta methods for index 2 DAEs.     

A parametrization method for the numerical solution of singular differential equations

The paper explains the numerical parametrization method (PM), originally created for optimal control problems, for classical calculus of variation problems that arise in connection with singular implicit (IDEs) and differential-algebraic equations (DAEs) in frame of their regularization. The PM for IDEs is based on representation of the required solution as a spline with moving knots and on minimization of the discrepancy functional with respect to the spline parameters. Such splines are named variational splines. For DAEs only finite entering functions can be represented by splines, and the functional under minimization is the discrepancy of the algebraic subsystem. The first and the second derivatives of the functionals are calculated in two ways – for DAEs with the help of adjoint variables, and for IDE directly. The PM does not use the notion of differentiation index, and it is applicable to any singular equation having a solution.     

Some remarks on solvability and various indices for implicit differential equations

It has been shown [17,18,21] that the notion of index for DAEs (Differential Algebraic Equations), or more generally implicit differential equations, could be interpreted in the framework of the formal theory of PDEs. Such an approach has at least two decisive advantages: on the one hand, its definition is not restricted to a “state-space” formulation (order one systems), so that it may be computed on “natural” model equations coming from physics (which can be, for example, second or fourth order in mechanics, second order in electricity, etc.) and there is no need to destroy this natural way through a first order rewriting. On the other hand, this formal framework allows for a straightforward generalization of the index to the case of PDEs (either “ordinary” or “algebraic”). In the present work, we analyze several notions of index that appeared in the literature and give a simple interpretation of each of them in the same general framework and exhibit the links they have with each other, from the formal point of view. Namely, we shall revisit the notions of differential, perturbation, local, global indices and try to give some clarification on the solvability of DAEs, with examples on time-varying implicit linear DAEs. No algorithmic results will be given here (see [34,35] for computational issues) but it has to be said that the complexity of computing the index, whatever approach is taken, is that of differential elimination, which makes it a difficult problem. We show that in fact one essential concept for our approach is that of formal integrability for usual DAEs and that of involution for PDEs. We concentrate here on the first, for the sake of simplicity. Last, because of the huge amount of work on DAEs in the past two decades, we shall mainly mention the most recent results. This revised version was published online in June 2006 with corrections to the Cover Date.     

热弹性动力学耦合问题的微分代数方法

对一般的热机械问题提出了一种有效的数值方法,并对二维的热弹性问题进行了测试．该方法的基本思路是将描述热机械耦合问题的偏微分方程进行降阶,使之成为一组微分代数方程,应力应变关系被写成代数方程．所得到的微分代数系统采用全隐式的向后差分公式进行求解．对该方法进行了详细的说明．为了验证该方法的有效性,将其应用于一个动态非耦合的热弹性问题的求解和一个耦合的二维热弹性问题的求解．