Distributional solution theory for linear DAEs

A solution theory for switched linear differential–algebraic equations (DAEs) is developed. To allow for non–smooth coordinate transformation, the coefficients matrices may have distributional entries. Since also distributional solutions are considered it is necessary to define a suitable multiplication for distribution. This is achieved by restricting the space of distributions to the smaller space of piecewise–smooth distributions. Solution formulae for two special DAEs, distributional ordinary differential equations (ODEs) and pure distributional DAEs, are given. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Constraint preserving integrators for general nonlinear higher index DAEs

Summary. In the last few years there has been considerable research on numerical methods for differential algebraic equations (DAEs) where is identically singular. The index provides one measure of the singularity of a DAE. Most of the numerical analysis literature on DAEs to date has dealt with DAEs with indices no larger than three. Even in this case, the systems were often assumed to have a special structure. Recently a numerical method was proposed that could, in principle, be used to integrate general unstructured higher index solvable DAEs. However, that method did not preserve constraints. This paper will discuss a modification of that approach which can be used to design constraint preserving integrators for general nonlinear higher index DAEs. Received August 25, 1993 / Revised version received April 7, 1994     

Robust Controllability of Linear Differential-Algebraic Equations with Unstructured Uncertainty

We consider the linear stationary systems of ordinary differential equations (ODEs) that are unsolvedwith respect to the derivative of the unknown vector-function and degenerate identically in the domain of definition. These systems are usually called differential-algebraic equations (DAEs). The measure of how a system of DAEs is unsolved with respect to the derivative is an integer which is called the index of the system of DAEs. The analysis is carried out under the assumption of existence of a structural form with separated differential and algebraic subsystems. We investigate the robust controllability of these systems (controllability in the conditions of uncertainty). The sufficient conditions for the robust complete and R-controllability of a system of DAEs with the indices 1 and 2 are obtained.     

On differential–algebraic equations in infinite dimensions

We consider a class of differential–algebraic equations (DAEs) with index zero in an infinite dimensional Hilbert space. We define a space of consistent initial values, which lead to classical continuously differential solutions for the associated DAE. Moreover, we show that for arbitrary initial values we obtain mild solutions for the associated problem. We discuss the asymptotic behaviour of solutions for both problems. In particular, we provide a characterisation for exponential stability and exponential dichotomies in terms of the spectrum of the associated operator pencil.     

Geometric Properties of Runge-Kutta Discretizations for Nonautonomous Index 2 Differential Algebraic Systems

We analyze Runge-Kutta discretizations applied to nonautonomous index 2 differential algebraic equations (DAEs) in semi-explicit form. It is shown that for half-explicit and projected Runge-Kutta methods there is an attractive invariant manifold for the discrete system which is close to the invariant manifold of the DAE. The proof combines reduction techniques to autonomou index 2 differential algebraic equations with some invariant manifold results of Schropp [9]. The results support the favourable behavior of these Runge-Kutta methods applied to index 2 DAEs for t = 0.     

拟线性微分代数方程的标准形

讨论了拟线性微分代数方程在一类特殊的奇点-拟障碍点附近的标准形.通过矩阵广义逆理论,拟线性微分代数方程可化为半显式形式.然后运用标准形理论,在微分同胚变换下,给出了拟线性微分代数方程在拟障碍点附近的标准形.在此基础上进一步讨论了这类标准形的去奇异化性质.     

几类微分-代数方程的神经网络求解法

在非线性科学中，寻求微分方程的近似解析解一直是重要的研究课题和研究热点.利用人工神经网络原理，结合最优化方法，研究了几类微分-代数方程的近似解析解，包括指标1，2，3型Hessenberg方程及指标3型Euler-Lagrange方程，得到了方程近似解析解的表达式.通过与精确解或Runge-Kutta(龙格-库塔)数值计算结果对比，表明神经网络方法的结果有很高的精度.     

Operator differential-algebraic equations with noise arising in fluid dynamics

We study linear semi-explicit stochastic operator differential algebraic equations (DAEs) for which the constraint equation is given in an explicit form. In particular, this includes the Stokes equations arising in fluid dynamics. We combine a white noise polynomial chaos expansion approach to include stochastic perturbations with deterministic regularization techniques. With this, we are able to include Gaussian noise and stochastic convolution terms as perturbations in the differential as well as in the constraint equation. By the application of the polynomial chaos expansion method, we reduce the stochastic operator DAE to an infinite system of deterministic operator DAEs for the stochastic coefficients. Since the obtained system is very sensitive to perturbations in the constraint equation, we analyze a regularized version of the system. This then allows to prove the existence and uniqueness of the solution of the initial stochastic operator DAE in a certain weighted space of stochastic processes.     

Numerical integration of rigid body dynamics in terms of quaternions

Unit–quaternions (or Euler parameter) are known to be well–suited for the singularity–free parametrization of finite rotations. Despite of this advantage, unit quaternions were rarely used to formulate the equations of motion (exceptions are the works by Nikravesh [1] and Haug [2]). This might be related to the fact, that the unit–quaternions are redundant, which requires the use of algebraic constraints in the equations of motion. Nowadays robust energy consistent integrators are available for the numerical solution of these differential–algebraic equations (DAEs). In the present work a mechanical integrator for the quaternions will be derived. This will be done by a size–reduction from the director formulation of the equations of motion, which also has the form of DAEs. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On Singular Values Related to DAEs in Kronecker Canonical Form

The index and the structural properties of differential algebraic equations (DAEs) are often determined by rank considerations of the derivative array. Since the Kronecker canonical form is a well-understood standard form that permits deep insight into the properties of DAEs, in this contribution we undertake an analysis of the singular values of this specific derivative array. To this end, the special structure of the obtained block matrices is pointed out, such that some formulas for the computation and estimation of eigenvalues and singular values can be applied. Actually, we explore the relationship between the spectra of particular block tridiagonal matrices and some perturbed Jacobi matrices.

     

Frequency Domain Methods and Decoupling of Linear Infinite Dimensional Differential Algebraic Systems

We discuss the analysis of linear constant coefficient differential algebraic equations on infinite dimensional Hilbert spaces. We give solution concepts and discuss solvability criteria which are mainly based on Laplace transform. Furthermore, we investigate the decoupling of these systems motivated by the Kronecker normal form for the finite dimensional case. Applications are given by the analysis of mixed systems of ordinary differential, partial differential and differential algebraic equations.     

Mixed symbolic–numerical computations with general DAEs I: System properties

Differential-algebraic equations (DAEs) arise in many ways in many types of problems. In this expository paper we discuss a variety of situations where we have found mixed symbolic-numerical calculations to be essential. The paper is designed to both familiarize the reader with several fundamental DAE ideas and to present some applications. The situations discussed include the analysis of DAEs, the solution of DAEs, and applications which include DAEs. Both successes and challenges will be presented. This revised version was published online in June 2006 with corrections to the Cover Date.     

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

A network model for incompressible two-fluid flow and its numerical solution

In this paper we consider an incompressible version of the two-fluid network model proposed by Porsching (Nu. Methods Part. Diff. Eq., 1 , 295–313 [1985]). The system of equations governing the model is a mixed system of differential and algebraic equations (DAEs). These DAEs are then recast, through proper transformation, into a system of ordinary differential equations on a submanifold of ?ⁿ, for which uniqueness, existence, and stability theorems are proved. Numerical simulations are presented.     

The parameterization method in optimal control problems and differential–algebraic equations

The paper is devoted to the explanation of the numerical parameterization method (PM) for optimal control (OC) problems with intermediate phase constraint and to its circumstantiation for classical calculus of variation (CV) problems that arise in connection with singular ODEs or DAEs, especially in cases of their essential degeneracy. The PM is based on a finite parameterization of control functions and on derivation of the problem functional with respect to control parameters. The first and the second derivatives are calculated with the help of adjoint vector and matrix impulses. Results of the solution to one phase constrained OC and two degenerate CV problems, connected with singular DAEs nonreducible to the normal form, are presented.     

Homoclinic orbits and Hopf bifurcations in delay differential systems with T-B singularity

The paper carries the results on Takens-Bogdanov bifurcation obtained in [T. Faria, L.T. Magalhães, Normal forms for retarded functional differential equations and applications to Bogdanov-Takens singularity, J. Differential Equations 122 (1995) 201-224] for scalar delay differential equations over to the case of delay differential systems with parameters. Firstly, we give feasible algorithms for the determination of Takens-Bogdanov singularity and the generalized eigenspace associated with zero eigenvalue in Rn. Next, through center manifold reduction and normal form calculation, a concrete reduced form for the parameterized delay differential systems is obtained. Finally, we describe the bifurcation behavior of the parameterized delay differential systems with T-B singularity in detail and present an example to illustrate the results.     

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.     

Consistent initialization for DAEs in Hessenberg form

For nonlinear DAEs, we can hardly make a reasonable statement unless structural assumptions are given. Many results are restricted to explicit DAEs, often in Hessenberg form of order up to three. For the DAEs resulting from circuit simulation, different beneficial structures have been found and exploited for the computation of consistent initial values. In this paper, a class of DAEs in nonlinear Hessenberg form of arbitrary high order is defined and analyzed with regard to consistent initialization. For this class of DAEs, the hidden constraints can be systematically described and the consistent initialization can be determined step-by-step solving linear subproblems, an approach hitherto used for the DAEs resulting from circuit simulation. Finally, it is shown that the DAEs resulting from mechanical systems fulfill the defined structural assumptions. The algorithm is illustrated by several examples.     

Efficient Reduced Order State Space Model Computation for a Class of Second Order Index One Systems

Simulation, design optimization and controller design of modern machine tools heavily rely on adequat numerical models. In order to achieve results in shorter computation times, reduced order models (ROMs) are applied in either of these tasks. Most modern simulation tools expect these ROMs to come in standard state space form. Structural models of the machine tool are however of second order type. In case piezo actuators are used in the device they are even differential algebraic equations (DAEs) of index one due to the coupling to the equations describing the electric potentials. This contribution is dedicated especially to those systems. We combine the ideas for balanced truncation model order reduction of large and sparse index 1 DAEs with methods developed for the efficient numerical handling of second order systems. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)