State event location in differential-algebraic models

An efficient discontinuity handling algorithm for initial value problems in differential-algebraic equations is presented. The algorithm supports flexible representation of state conditions in propositional logic, and guarantees the location of all state events in strict time order. The algorithm consists of two phases:(1) event detection and(2) consistent event location. In the event detection phase, the entire integration step is searched for the state event by solving the interpolation polynomials for the discontinuity functions generated by the BDF method. An efficient hierarchical polynomial root-finding procedure based upon interval arithmetic guarantees detection of the state event even if multiple state condition transitions exist in an integration step, in which case many existing algorithms may fail. As a second phase of the algorithm, a consistent even location calculation is developed that accurately locates the state event detected earlier while completely eliminating incorrect reactivation of the same state event immediately after the consistent initialization calculation that may follow. This numerical phenomenon has not been explained before and is termed discontinuity sticking. Results from various test problems are presented to demonstrate the correctness and efficiency of the algorithm.

[1]  D. Ellison,et al.  Efficient automatic integration of ordinary differential equations with discontinuities , 1981 .

[2]  Robert L. Grossman,et al.  Timed Automata , 1999, CAV.

[3]  Constantinos C. Pantelides,et al.  SpeedUp—recent advances in process simulation , 1988 .

[4]  L. Dekker,et al.  Simulation of systems '79 : proceedings of the 9th IMACS Congress, Sorrento, Italy, September 24-28, 1979 , 1980 .

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

[6]  Ramon E. Moore.  A Test for Existence of Solutions to Nonlinear Systems , 1977 .

[7]  Paul I. Barton,et al.  Modeling of combined discrete/continuous processes , 1994 .

[8]  Ole Østerby,et al.  Solving Ordinary Differential Equations with Discontinuities , 1984, TOMS.

[9]  Gholamreza Torkzadeh,et al.  State event in combined simulation , 1988, WSC '88.

[10]  Lawrence F. Shampine,et al.  Reliable solution of special event location problems for ODEs , 1991, TOMS.

[11]  Tuncer I. Ören,et al.  A robust procedure for discontinuity handling in continuous system simulation , 1985 .

[12]  Sven Erik Mattsson,et al.  On modelling and differential/algebraic systems , 1989, Simul..

[13]  François E. Cellier.  Combined continuous/discrete system simulation by use of digital computers , 1979 .

[14]  Gintaras V. Reklaitis,et al.  A simulator for batch and semi-conductors processes☆ , 1984 .

[15]  M. B. Carver.  Efficient integration over discontinuities in ordinary differential equation simulations , 1978 .

[16]  Constantinos C. Pantelides,et al.  The Dynamic Simulation of Transient Systems Described by Index Two Differential-algebraic Equations , 1988 .

[17]  Ramon E. Moore.  Methods and applications of interval analysis , 1979, SIAM studies in applied mathematics.

[18]  M. Berzins,et al.  Algorithms for the location of discontinuities in dynamic simulation problems , 1991 .

[19]  H. Elmqvist,et al.  Automated formula manipulation supports object-oriented continuous-system modeling , 1993, IEEE Control Systems.

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

[21]  Richard E. Nance,et al.  The time and state relationships in simulation modeling , 1981, CACM.

[22]  A. Neumaier.  Interval methods for systems of equations , 1990 .

[23]  ParkTaeshin,et al.  State event location in differential-algebraic models , 1996 .

[24]  E. Hansen,et al.  Bounding solutions of systems of equations using interval analysis , 1981 .