Extending explicit and linearly implicit ODE solvers for index-1 DAEs

a b s t r a c t Nonlinear differential-algebraic equations (DAE) are typically solved using implicit stiff solvers based on backward difference formula or RADAU formula, requiring a Newton-Raphson approach for the nonlin- ear equations or using Rosenbrock methods specifically designed for DAEs. Consistent initial conditions are essential for determining numeric solutions for systems of DAEs. Very few systems of DAEs can be solved using explicit ODE solvers. This paper applies a single-step approach to system initialization and simulation allowing for systems of DAEs to be solved using explicit (and linearly implicit) ODE solvers without a priori knowledge of the exact initial conditions for the algebraic variables. Along with using a combined process for initialization and simulation, many physical systems represented through large systems of DAEs can be solved in a more robust and efficient manner without the need for nonlinear solvers. The proposed approach extends the usability of explicit and linearly implicit ODE solvers and removes the requirement of Newton-Raphson type iteration. Published by Elsevier Ltd.

[1]  E. Süli,et al.  Numerical Solution of Ordinary Differential Equations , 2021, Foundations of Space Dynamics.

[2]  J. Dormand,et al.  A family of embedded Runge-Kutta formulae , 1980 .

[3]  Wayne H. Enright,et al.  Interpolants for Runge-Kutta formulas , 1986, TOMS.

[4]  Venkatasailanathan Ramadesigan,et al.  A perturbation approach for consistent initialization of index-1 explicit differential-algebraic equations arising from battery model simulations , 2011, Comput. Chem. Eng..

[5]  Michael L. Michelsen,et al.  Application of semi-implicit Runge—Kutta methods for integration of ordinary and partial differential equations , 1977 .

[6]  Paul W. C. Northrop,et al.  Efficient Simulation and Reformulation of Lithium-Ion Battery Models for Enabling Electric Transportation , 2014 .

[7]  李青,et al.  Wolfram综合征二例 , 2005 .

[8]  Lawrence F. Shampine,et al.  Solving Index-1 DAEs in MATLAB and Simulink , 1999, SIAM Rev..

[9]  R. Sargent,et al.  The mathematical modelling of transient systems using differential-algebraic equations , 1988 .

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

[11]  Yaoyu Li,et al.  Consistent initialization of system of differential-algebraic equations for dynamic simulation of centrifugal chillers , 2012 .

[12]  L. Shampine,et al.  Numerical Solution of Ordinary Differential Equations. , 1995 .

[13]  François E. Cellier,et al.  Continuous System Simulation , 2006 .

[14]  A. Hindmarsh.  LSODE and LSODI, two new initial value ordinary differential equation solvers , 1980, SGNM.

[15]  Matthew B. Pinson,et al.  Theory of SEI Formation in Rechargeable Batteries: Capacity Fade, Accelerated Aging and Lifetime Prediction , 2012, 1210.3672.

[16]  Martin Berzins,et al.  Developing software for time-dependent problems using the method of lines and differential-algebraic integrators , 1989 .

[17]  Venkat R. Subramanian,et al.  Model-Based SEI Layer Growth and Capacity Fade Analysis for EV and PHEV Batteries and Drive Cycles , 2014 .

[18]  Venkatasailanathan Ramadesigan,et al.  Continuum Representation for Simulating Discrete Events of Battery Operation , 2010 .

[19]  Ignacio E. Grossmann,et al.  Computers and Chemical Engineering , 2014 .

[20]  Ralph E. White,et al.  Review of Models for Predicting the Cycling Performance of Lithium Ion Batteries , 2006 .

[21]  Alan H. Karp,et al.  A variable order Runge-Kutta method for initial value problems with rapidly varying right-hand sides , 1990, TOMS.

[22]  Linda R. Petzold,et al.  Using Krylov Methods in the Solution of Large-Scale Differential-Algebraic Systems , 1994, SIAM J. Sci. Comput..

[23]  Lawrence F. Shampine,et al.  The MATLAB ODE Suite , 1997, SIAM J. Sci. Comput..

[24]  Peter M. Dew,et al.  A Set of Library Routines for Solving Parabolic Equations in One Space Variable , 1981, TOMS.

[25]  René Lamour,et al.  Computation of consistent initial values for properly stated index 3 DAEs , 2009 .

[26]  L. Petzold.  A description of dassl: a differential/algebraic system solver , 1982 .

[27]  Venkatasailanathan Ramadesigan,et al.  Coordinate Transformation, Orthogonal Collocation, Model Reformulation and Simulation of Electrochemical-Thermal Behavior of Lithium-Ion Battery Stacks , 2011 .

[28]  Holger Boche,et al.  On inconsistent initial conditions for linear time-invariant differential-algebraic equations , 2002, 2002 IEEE International Symposium on Circuits and Systems. Proceedings (Cat. No.02CH37353).

[29]  Venkat R. Subramanian,et al.  Computational Methods in Chemical Engineering with Maple , 2010 .

[30]  J. Cash.  Modified extended backward differentiation formulae for the numerical solution of stiff initial value problems in ODEs and DAEs , 2000 .

[31]  C. W. Gear,et al.  Approximation methods for the consistent initialization of differential-algebraic equations , 1988 .

[32]  Linda R. Petzold,et al.  Consistent Initial Condition Calculation for Differential-Algebraic Systems , 1998, SIAM J. Sci. Comput..

[33]  J. R. Cash,et al.  Review Paper: Efficient numerical methods for the solution of stiff initial-value problems and differential algebraic equations , 2003, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[34]  Robert M. Corless,et al.  Initial value problems for ODEs in problem solving environments , 2000 .

[35]  P. I. Barton,et al.  On inconsistent initial conditions for linear time-invariant differential-algebraic equations , 2002 .

[36]  Lawrence F. Shampine,et al.  Solving 0 = F(t, y(t), y′(t)) in Matlab , 2002, J. Num. Math..

[37]  L. Petzold,et al.  Software and algorithms for sensitivity analysis of large-scale differential algebraic systems , 2000 .

[38]  Takashi Hikihara,et al.  Discontinuous Dynamics of Electric Power System With DC Transmission: A Study on DAE System , 2008, IEEE Transactions on Circuits and Systems I: Regular Papers.

[39]  Richard D. Braatz,et al.  Modeling and Simulation of Lithium-Ion Batteries from a Systems Engineering Perspective , 2010 .

[40]  C. Schneider.  Rosenbrock-type methods adapted to differential-algebraic systems , 1991 .

[41]  J. Verner.  Explicit Runge–Kutta Methods with Estimates of the Local Truncation Error , 1978 .

[42]  L. Petzold,et al.  DASPK: A new high order and adaptive time-integration technique with applications to mantle convection with strongly temperature-and pressure-dependent rheology , 1995 .

[43]  Ralph E. White,et al.  An initialization subroutine for DAEs solvers: DAEIS , 2001 .

[44]  E. Hairer,et al.  Stiff differential equations solved by Radau methods , 1999 .

[45]  Lorenz T. Biegler,et al.  A Successive Linear Programming Approach for Initialization and Reinitialization after Discontinuities of Differential-Algebraic Equations , 1998, SIAM J. Sci. Comput..

[46]  Walter M. Lioen,et al.  Specification of PSIDE , 1998 .

[47]  K. L. Praprost,et al.  A stability theory for constrained dynamic systems with applications to electric power systems , 1996 .

[48]  Linda R. Petzold,et al.  Differential-algebraic equations , 2008, Scholarpedia.

[49]  Gonzalez Garcia,et al.  A singular perturbation approach to modeling closed kinematic chains , 2000 .

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