On modelling and differential/algebraic systems

Abstract

It is of a vital interest that software for model development and simulation supports models given as differential/algebraic equa tions (DAE), since many models for physical systems are naturally given as sets of ordinary differential and algebraic equations and since, as motivated in the report, the equation form is the only reasonable for model libraries. Unfortunately, most simulation packages of today do not allow models given as DAE systems, but require assignment statements for derivatives and algebraic variables.

There are numerical solvers for DAE systems, but they are not as reliable and robust as numerical solvers for ordinary differential equations on state space form. Unfortunately, DAE systems may exhibit bad properties, which make them difficult to solve numerically today. The purpose of the paper is to give an over view of important properties of DAE systems, indicate possible ways to handle difficulties and to motivate others to do research in this area.

By means of examples it is illustrated that even small, simple DAE systems can exhibit bad properties. The examples are by no means pathological. They arise naturally during model development.

There are many open questions to answer. The questions cover a broad range from theoretical to practical problems. To en courage research in this area, interesting issues are listed.

Keywords

modelling simulation CSSL,differential/algebraic systems

Get full access to this article

View all access options for this article.

References

Acsl (1986): Advanced Continuous Simulation Language (ACSL)-Reference Manual, Fourth Edition , Mitchell and Gauthier Associates, Concord, Mass. 01742, USA.

Åstrom, K.J. (1983): "Computer Aided Modelling, Analysis and Design of Control Systems - A perspective," IEEE Control Systems Magazine, Vol. 3, No. 2, May 1983, pp. 4-16.

Åstrom, K.J. (1985): "Computer Aided Tools for Control System Design -A perspective," In Janshidi and Herget (1985 ).

Barrlund, A. (1985): "Numerisk behandling av differentialalgebraiska system (Numeric solution of differential/algebraic systems)," UMNAD-16.85, Institute of Information Processing, University of Umea, Umea, Sweden.

Cellier, F.A. (1979): Combined Continuous/Discrete System Simulation by Use of Digital Computers: Techniques and Tools, Ph.D-thesis, Diss ETH No 6483, The Swiss Federal Institute of Technology, ETH, Zurich.

Chambers, T.H.E. and M.J. Whitmarsh-Everiss (1984): "A Methodology for the design of Plant Analysers," PKR/SE/201. Central Electricity Generating Board, U.K.

Coddington, E.A. and N. Levinson (1955): The Theory of Ordinary Differential Equations , McGraw-Hill Book Company, Inc., New York.

Elmqvist, H. (1975): "SIMNON - User's Manual," Rept. TFRT-3091. Department of Automatic Control, Lund Institute of Technology , Sweden.

Elmqvist, H. (1978): "A Structures Model Language for Large Continuous Systems," CODEN: LUTED2/TFRT1015, Department of Automatic Control, Lund Institute of Technology, Lund, Sweden.

10.

Gantmacher, F.R. (1959): The Theory of Matrices, Vol. 1-II, Chelsea Publishing Company, New York, N.Y., USA.

11.

Gear, C.W. and L.R. Petzold (1984): "ODE Methods for the Solution of Differential /Algebraic Systems," SIAM J. Numerical Analysis, Vol. 21, No. 4, Aug. 1984, pp. 716-728.

12.

Jamshidi M. and C.J. Herget (Eds.) (1985): Computer-Aided Control Systems Engineering, North-Holland.

13.

Marz, R. (1982a): "Multistep Methods for Initial Value Problems in Implicit Differential-Algebraic Equations," Preprint 22, Sektion Mathematik, Humbolt-University, Berlin.

14.

Marz, R. (1982b): "On Difference and Shooting Methods for Boundary Value Problems in Differential-Algebraic Equations," Preprint 24, Sektion Mathematik, Humbolt-University, Berlin.

15.

Petzold, L.R. (1982a): "Differential-Algebraic Systems are not ODE'S ," SIAM Journal on Scientific and Statistical Computing, Vol. 3, No. 3, Sept. 1982, pp. 367-384.

16.

Petzold, L.R. (19826): "A Description of DASSL: A Differential/Algebraic System Solver," Rept. SAND 82-8637, Sandia National Laboratories or Proceedings of IMACS World Congress, Montreal , Canada, 1982.

17.

Soderlind, G. (1980): "DASP3 - A Program for the Numerical Integration of Partitioned Stiff ODEs and Differential-Algebraic Systems," Rept. TRITA-NA-8008. Department of Numerical Analysis and Computing Science, The Royal Institute of Technology, Stockholm, Sweden.

18.

Tarjan, R.E. (1972): "Depth First Search and Linear Graph Algorithms ," SIAM J. Comp., 1, 146-160.

19.

Wiberg, T. (1977): "Permutation of an Unsymmetric Matrix to Block Triangular Form," Diss., Department of Information Processing, University of Umeå, Umeå, Sweden.