Abstract
Automated modeling of multi-physics dynamical systems often results in large-scale high-index differential-algebraic equations (DAEs). Since direct numerical simulation of such systems leads to instabilities and possibly non-convergence of numerical methods, a regularization or remodeling is required. In many simulation environments, a structural analysis based on the sparsity pattern of the system is used to determine the index and an index-reduced system model. Here, usually the Pantelides algorithm in combination with the Dummy Derivative Method is used. We present a new approach for the regularization of DAEs that is based on the Signature method (\(\varSigma \)-method).
Similar content being viewed by others
Notes
The term structurally singular is also used slightly differently in the literature: on the one hand, it is used for systems that do not allow an assignment of the highest occurring derivative of each variables to a specific equation (i.e., systems that do not have a transversal with finite value) [18], and on the other hand it is used for systems with singular \(\varSigma \)-Jacobian [16].
References
Altmeyer, R., Steinbrecher, A.: Regularization and Numerical Simulation of Dynamical Systems Modeled with Modelica. Institut für Mathematik, TU Berlin, Preprint 29–2013 (2013)
Barton, P.I., Martinson, W.S., Reißig, G.: Differential-algebraic equations of index 1 may have an arbitrarily high structural index. SIAM J. Sci. Comput. 21, 1987–1990 (2000)
Brenan, K., Campbell, S., Petzold, L.: Numerical Solution of Initial-Value Problems in Differential Algebraic Equations, Classics in Applied Mathematics, vol. 14. SIAM, Philadelphia (1996)
Campbell, S., Gear, C.: The index of general nonlinear DAEs. Numerische Mathematik 72(2), 173–196 (1995)
Fu, Z., Eglese, R., Wright, M.: A branch-and-bound algorithm forfinding all optimal solutions of the assignment problem. Asia Pac. J. Oper. Res. 24(06), 831–839 (2007)
Fukuda, K., Matsui, T.: Finding all the perfect matchings in bipartite graphs. Appl. Math. Lett. 7(1), 15–18 (1994)
Fritzson, P.: Principles of Object-Oriented Modeling and Simulation with Modelica 2.1. Wiley-IEEE Press, New York (2004)
Hairer, E., Wanner, G.: Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems, 2nd edn. Springer-Verlag, Berlin (1996)
Kuhn, H.W.: The Hungarian Method for the assignment problem. Naval Res. Logistics Q. 2, 83–97 (1955)
Kunkel, P., Mehrmann, V.: Analysis of over- and underdetermined nonlinear differential-algebraic systems with application to nonlinear control problems. Math. Control Sign. Syst. 14, 233–256 (2001)
Kunkel, P., Mehrmann, V.: Index reduction for differential-algebraic equations by minimal extension. Zeitschrift für Angewandte Mathematik und Mechanik 84(9), 579–597 (2004)
Kunkel, P., Mehrmann, V.: Differ. Algebraic Equ. Anal. Numer. Solut. EMS Publishing House, Zürich (2006)
Leitold, A., Hangos, K.M.: Structural solvability analysis of dynamic process models. Comput. Chem. Eng. 25(11), 1633–1646 (2001)
Mattsson, S., Söderlind, G.: Index reduction in differential-algebraic equations using dummy derivatives. SIAM J. Sci. Stat. Comput. 14, 677–692 (1993)
Murota, K.: Matrices and Matroids for Systems Analysis. Algorithms and Combinatorics, 20th edn. Springer, New York (2000)
Nedialkov, N., Pryce, J.: Solving differential-algebraic equations by Taylor series (i): computing Taylor coefficients. BIT Numer. Math. 45(3), 561–591 (2005)
Nedialkov, N., Pryce, J., Tan, G.: DAESA: a Matlab tool for structural analysis of DAEs: Software. Technical report CAS-12-01-NN, Department of Computing and Software, McMaster University, Hamilton (2012)
Nilsson, H.: Type-based structural analysis for modular systems of equations. In: Proceedings of the 2nd International Workshop on Equation-Based Object-Oriented Languages and Tools 029, 71–81 (2008)
Pantelides, C.: The consistent initialization of differential-algebraic systems. SIAM J. Sci. Stat. Comput. 9, 213–231 (1988)
Pryce, J.: A simple structural analysis method for DAEs. BIT Numer. Math. 41, 364–394 (2001)
Scholz, L., Steinbrecher, A.: A combined structural-algebraic approach for the regularization of coupled systems of DAEs. Institut für Mathematik, TU Berlin, Preprint 30–2013 (2013)
Scholz, L., Steinbrecher, A.: Efficient numerical integration of dynamical systems based on structural-algebraic regularization avoiding state selection. In: Proceedings of the 10th International Modelica Conference (2014)
Scholz, L., Steinbrecher, A.: Structural-algebraic regularization for coupled systems of DAEs. Submitted to BIT Numerical Mathematics (2014)
Soares, R.P., Secchi, A.R.: Direct initialisation and solution of high-index DAE systems. Comput. Aided Chem. Eng. 20, 157–162 (2005)
Zeng, Y., Wu, X., Cao, J.: An improved KM algorithm for computing structural index of DAE system. In: 12th International Symposium on Distributed Computing and Applications to Business, Engineering and Science (DCABES), IEEE, pp. 95–99 (2013)
Acknowledgments
We would like to thank two anonymous referees for their careful reading and for thoughtful suggestions for the improvement of the paper.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Hans-Petter Langtangen.
Research supported by European Research Council, through ERC Advanced Grant MODSIMCONMP.
Rights and permissions
About this article
Cite this article
Scholz, L., Steinbrecher, A. Regularization of DAEs based on the Signature method. Bit Numer Math 56, 319–340 (2016). https://doi.org/10.1007/s10543-015-0565-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10543-015-0565-x
Keywords
- Differential-algebraic equation
- Regularization
- Structural analysis
- \(\varSigma \)-method
- Index reduction