Reversible methods of Runge-Kutta type for Index-2 DAEs

Summary. A new interpretation of Runge-Kutta methods for differential algebraic equations (DAEs) of index 2 is presented, where a step of the method is described in terms of a smooth map (smooth also with respect to the stepsize). This leads to a better understanding of the convergence behavior of Runge-Kutta methods that are not stiffly accurate. In particular, our new framework allows for the unified study of two order-improving techniques for symmetric Runge-Kutta methods (namely post-projection and symmetric projection) specially suited for solving reversible index-2 DAEs.Mathematics Subject Classification (1991): 65L05, 65L06     

Numerical Integration of Differential Algebraic Systems and Invariant Manifolds

The dynamics of a differential algebraic equation takes place on a lower dimensional manifold in phase space. Applying a numerical integration scheme, it is natural to ask if and how this geometric property is preserved by the discrete dynamical system. In the index-1 case answers to this question are obtained from the singularly perturbed case treated by Nipp and Stoffer, Numer. Math. 70 (1995), 245–257, for Runge-Kutta methods and in K. Nipp and D. Stoffer, Numer. Math. 74 (1996), 305–323, for linear multistep methods. As main result of this paper it is shown that also for Runge-Kutta methods and linear multistep methods applied to a index-2 problem of Hessenberg form there is a (attractive) invariant manifold for the discrete dynamical system and this manifold is close to the manifold of the differential algebraic equation.     

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.     

Efficient Runge-Kutta integrators for index-2 differential algebraic equations

In seeking suitable Runge-Kutta methods for differential algebraic equations, we consider singly-implicit methods to which are appended diagonally-implicit stages. Methods of this type are either similar to those of Butcher and Cash or else allow for the importation of a final derivative from a previous step. For these two classes, with up to three additional diagonally-implicit stages, we derive methods that satisfy appropriate order conditions for index-2 DAEs.

     

Higher order fully discrete scheme combined withH 1-Galerkin mixed finite element method for semilinear reaction-diffusion equations

We first apply a first order splitting to a semilinear reaction-diffusion equation and then discretize the resulting system by anH ¹-Galerkin mixed finite element method in space. This semidiscrete method yields a system of differential algebraic equations (DAEs) ofindex one. Apriori error estimates for semidiscrete scheme are derived for both differential as well as algebraic components. For fully discretization, an implicit Runge-Kutta (IRK) methods is applied to the temporal direction and the error estimates are discussed for both components. Finally, we conclude the paper with a numerical example.     

MULTISTEP DISCRETIZATION OF INDEX 3 DAES

1. IntroductionIn this paper, we will consider the multistep discrezations of the differential--algebraicequations (DAEs) in Hessenberg formwhere F e AN M M R", K E AN M L - AM, G E AN - RL, the initial value(yo, ic, no) at xo are assumed to be consistent, i.e., they satisfyWe supposes, F, G and K are sufficiently differentiable, and thatin a neighbourhood of the solution. Such problems often appear in the simulation ofmechanical systems with constraints and the singularly perturbed…     

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.     

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)     

Completions of nonlinear DAE flows based on index reduction techniques and their stabilization

Differential algebraic equations (DAEs) define a differential equation on a manifold. A number of ways have been developed to numerically solve some classes of DAEs. Motivated by problems in control theory, numerical simulation, and the use of general purpose modeling environments, recent research has considered the embedding of the DAE solutions of a general DAE into the solutions of an ODE where the added dynamics have special properties. This paper both provides new results on the linear time-varying case and considers the important nonlinear case.     

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     

Runge-Kutta Discretizations of Singularly Perturbed Gradient Equations

We analyze Runge-Kutta discretizations applied to singularly perturbed gradient systems. It is shown in which sense the discrete dynamics preserve the geometric properties and the longtime behavior of the underlying ordinary differential equation. If the continuous system has an attractive invariant manifold then numerical trajectories started in some neighbourhood (the size of which is independent of the step-size and the stiffness parameter) approach an equilibrium in a nearby manifold. The proof combines invariant manifold techniques developed by Nipp and Stoffer for singularly perturbed systems with some recent results of the second author on the global behavior of discretized gradient systems. The results support the favorable behavior of ODE methods for stiff minimization problems.     

Runge-Kutta schemes that maintain invariant solutions for differential algebraic systems

Implicit Runge-Kutta (IRK) methods and projected IRK methods for the solution of semiexplicit index-2 systems of differential algebraic systems (DAEs) have been proposed by several authors. In this paper we prove that if a method satisfiesBA+A ^t B–bb ^t =0, it conserves quadratic invariants of DAEs.     

New Multivalue Methods for Differential Algebraic Equations

Multivalue methods are slightly different from the general linear methods John Butcher proposed over 30 years ago. Multivalue methods capable of solving differential algebraic equations have not been developed. In this paper, we have constructed three new multivalue methods for solving DAEs of index 1, 2 or 3, which include multistep methods and multistage methods as special cases. The concept of stiff accuracy will be introduced and convergence results will be given based on the stage order of the methods. These new methods have the diagonal implicit property and thus are cheap to implement and will have order 2 or more for both the differential and algebraic components. We have implemented these methods with fixed step size and they are shown to be very successful on a variety of problems. Some numerical experiments with these methods are presented.     

Half-explicit Runge-Kutta methods with explicit stages for differential-algebraic systems of index 2

Usually the straightforward generalization of explicit Runge-Kutta methods for ordinary differential equations to half-explicit methods for differential-algebraic systems of index 2 results in methods of orderq≤2. The construction of higher order methods is simplified substantially by a slight modification of the method combined with an improved strategy for the computation of the algebraic solution components. We give order conditions up to orderq=5 and study the convergence of these methods. Based on the fifth order method of Dormand and Prince the fifth order half-explicit Runge-Kutta method HEDOP5 is constructed that requires the solution of 6 systems of nonlinear equations per step of integration.     

Parareal for index two differential algebraic equations

Numerical Algorithms - This article proposes modifications of the Parareal algorithm for its application to higher index differential algebraic equations (DAEs). It is based on the idea of applying...     

求解延迟微分代数方程的多步Runge-Kutta方法的渐近稳定性

延迟微分代数方程(DDAEs)广泛出现于科学与工程应用领域．本文将多步Runge-Kutta方法应用于求解线性常系数延迟微分代数方程，讨论了该方法的渐近稳定性．数值试验表明该方法对求解DDAEs是有效的．     

An efficient asymptotically correct error estimator for collocation solutions to singular index-1 DAEs

A computationally efficient a posteriori error estimator is introduced and analyzed for collocation solutions to linear index-1 DAEs (differential-algebraic equations) with properly stated leading term exhibiting a singularity of the first kind. The procedure is based on a modified defect correction principle, extending an established technique from the context of ordinary differential equations to the differential-algebraic case. Using recent convergence results for stiffly accurate collocation methods, we prove that the resulting error estimate is asymptotically correct. Numerical examples demonstrate the performance of this approach. To keep the presentation reasonably self-contained, some arguments from the literature on DAEs concerning the decoupling of the problem and its discretization, which is essential for our analysis, are also briefly reviewed. The appendix contains a remark about the interrelation between collocation and implicit Runge-Kutta methods for the DAE case.     

Stabilization of invariants of discretized differential systems

Many problems of practical interest can be modeled by differential systems where the solution lies on an invariant manifold defined explicitly by algebraic equations. In computer simulations, it is often important to take into account the invariant's information, either in order to improve upon the stability of the discretization (which is especially important in cases of long time integration) or because a more precise conservation of the invariant is needed for the given application. In this paper we review and discuss methods for stabilizing such an invariant. Particular attention is paid to post-stabilization techniques, where the stabilization steps are applied to the discretized differential system. We summarize theoretical convergence results for these methods and describe the application of this technique to multibody systems with holonomic constraints. We then briefly consider collocation methods which automatically satisfy certain, relatively simple invariants. Finally, we consider an example of a very long time integration and the effect of the loss of symplecticity and time-reversibility by the stabilization techniques. This revised version was published online in June 2006 with corrections to the Cover Date.     

Convergence of multistep discretizations of DAEs

Standard ODE methods such as linear multistep methods encounter difficulties when applied to differential-algebraic equations (DAEs) of index greater than 1. In particular, previous results for index 2 DAEs have practically ruled out the use of all explicit methods and of implicit multistep methods other than backward difference formulas (BDFs) because of stability considerations. In this paper we embed known results for semi-explicit index 1 and 2 DAEs in a more comprehensive theory based on compound multistep and one-leg discretizations. This explains and characterizes the necessary requirements that a method must fulfill in order to be applicable to semi-explicit DAEs. Thus we conclude that the most useful discretizations are those that avoid discretization of the constraint. A freer use of e.g. explicit methods for the non-stiff differential part of the DAE is then possible.Dedicated to Germund Dahlquist on the occasion of his 70th birthdayThis author thanks the Centro de Estadística y Software Matemático de la Universidad Simón Bolivar (CESMa) for permitting her free use of its research facilities.Partial support by the Swedish Research Council for Engineering Sciences TFR under contract no. 222/91-405.