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.

[1]  C. W. Gear,et al.  Simultaneous Numerical Solution of Differential-Algebraic Equations , 1971 .

[2]  J. Baumgarte.  Stabilization of constraints and integrals of motion in dynamical systems , 1972 .

[3]  E. Haug,et al.  Generalized Coordinate Partitioning for Dimension Reduction in Analysis of Constrained Dynamic Systems , 1982 .

[4]  L. Petzold.  Differential/Algebraic Equations are not ODE's , 1982 .

[5]  W. Rheinboldt.  Differential-algebraic systems as differential equations on manifolds , 1984 .

[6]  U. Ascher,et al.  Stability of collocation at Gaussian points , 1984 .

[7]  C. W. Gear,et al.  Maintianing solution invariants in the numerical solution of ODEs , 1986 .

[8]  G. J. Cooper.  Stability of Runge-Kutta Methods for Trajectory Problems , 1987 .

[9]  C. W. Gear,et al.  Differential algebraic equations, indices, and integral algebraic equations , 1990 .

[10]  B. Leimkuhler,et al.  Stabilization and projection methods for multibody dynamics , 1990 .

[11]  U. Ascher,et al.  Projected implicit Runge-Kutta methods for differential-algebraic equations , 1990 .

[12]  B. Leimkuhler,et al.  Numerical solution of differential-algebraic equations for constrained mechanical motion , 1991 .

[13]  F. Potra,et al.  On the Numerical Solution of Euler-Lagrange Equations∗ , 1991 .

[14]  M. C. Recchioni,et al.  A Quadratically Convergent Method for Linear Programming , 1991 .

[15]  Hans Bock,et al.  Recent Advances in the Numerical Integration of Multibody Systems , 1993 .

[16]  Uri M. Ascher,et al.  Stability of Computational Methods for Constrained Dynamics Systems , 1991, SIAM J. Sci. Comput..

[17]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[18]  E. Eich.  Convergence results for a coordinate projection method applied to mechanical systems with algebraic constraints , 1993 .

[19]  J. M. Sanz-Serna,et al.  Numerical Hamiltonian Problems , 1994 .

[20]  U. Ascher,et al.  Stabilization of DAEs and invariant manifolds , 1994 .

[21]  R. Russell,et al.  Unitary integrators and applications to continuous orthonormalization techniques , 1994 .

[22]  B. Leimkuhler,et al.  Symplectic Numerical Integrators in Constrained Hamiltonian Systems , 1994 .

[23]  W. Rheinboldt,et al.  On the numerical solution of the Euler-Lagrange equations , 1995 .

[24]  U. Ascher,et al.  Stabilization of Constrained Mechanical Systems with DAEs and Invariant Manifolds , 1995 .

[25]  E. Hairer,et al.  Accurate long-term integration of dynamical systems , 1995 .

[26]  Hong Sheng Chin,et al.  Stabilization methods for simulations of constrained multibody dynamics , 1995 .

[27]  S. Reich.  Enhancing energy conserving methods , 1996 .

[28]  Ernst Hairer,et al.  Reversible Long-Term Integration with Variable Stepsizes , 1997, SIAM J. Sci. Comput..

[29]  Antonella Zanna,et al.  Numerical solution of isospectral flows , 1997, Math. Comput..

[30]  D. Higham.  Time-stepping and preserving orthonormality , 1997 .

[31]  Uri M. Ascher,et al.  Sequential Regularization Methods for Nonlinear Higher-Index DAEs , 1997, SIAM J. Sci. Comput..

[32]  L. Shampine.  Conservation laws and the numerical solution of ODEs, II , 1999 .

[33]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[34]  M. A. Akanbi,et al.  Numerical solution of initial value problems in differential - algebraic equations , 2005 .