Abstract
It is well-known that linear time-varying high-index DAEs can be very sensitive to parametric perturbations, [1, p. 31]. Stability is also affected, as is known from singular perturbation theory. In this note, we show that arbitrarily small and smooth perturbations can cause dramatic instabilities by introducing “small” perturbations of the matrix pencil's generalized eigenvalues atz=∞, leading to large positive finite eigenvalues. The smaller the perturbation, the larger is the instability of the perturbed problem, in contrast to the ODE case. Some high-index problems can thus be considered as marginally stable, with neighboring problems (usually of lower index) exhibiting severe instabilities.
Zusammenfassung
Es ist bekannt, daß Algebrodifferentialgleichungen von höherem Index sehr empfindlich gegenüber Parameteränderungen sein können [1, p. 31]. Ebenso ist aus der Theorie der singulär gestörten Probleme bekannt, daß auch die Stabilitätseigenschaften bei derartigen Störungen sich ändern können. In diesem Aufsatz zeigen wir, daß beliebig kleine, glatte Störungen Instabilitäten dadurch verursachen können, daß die verallgemeinerten Eigenwerte des zugehörigen Matrizenbüschels beiz=∞ zu großen positiven Eigenwerte verfälscht werden. Das Problem wird dadurch umso instabiler je kleiner die Störung des Ausgangsproblems ist, ganz im Gegensatz zu der Situation bei expliziten DGLn. Auf diese Weise können einige Probleme von höherem Index als ein stabiler Grenzfall von benachbarten hoch instabilen Problemen, gewöhnlich von niedrigerem Index, angesehen werden.
Similar content being viewed by others
References
Brenan, K. E., Campbell, S. L., Petzold, L. R.: Numerical solution of initial-value problems in differential-algebraic equations. New York: North-Holland 1989.
Clark, K. D.: The numerical solution of linear time varying singular systems of differential equations by difference methods. Ph.D. thesis, Dept. of Math, North Carolina State Univ., Raleigh, NC, USA, 1985.
Dahlquist, G., Edsberg, L., Sköllermo, G., Söderlind, G.: Are the numerical methods and software satisfactory for chemical kinetics? In: Hinze, J. (ed.) Numerical integration of differential equations and large linear systems. Berlin, Heidelberg: Springer 1982.
Duff, I. S., Gear, C. W.: Computing the structural index. SIAM J. Alg. Discrete Methods7, 594–603 (1986).
Hairer, E., Lubich, Ch., Roche, M.: The numerical solution of differential-algebraic systems by Runge-Kutta methods. Berlin, Heidelberg: Springer 1989 (Springer Lecture Notes in Mathematics, 1409).
Hairer, E., Wanner, G.: Solving ordinary differential equations II. Berlin, Heidelberg: Springer 1991.
Hartman, P.: Ordinary differential equations, 2nd ed. Boston: Birkäuser 1982.
O'Malley, R. E., Jr.: Introduction to singular perturbations. New York: Academic Press 1974.
O'Malley, R. E., Jr.: On nonlinear singularly perturbed initial value problems. SIAM Review30, 193–212 (1988).
Author information
Authors and Affiliations
Additional information
Dedicated to the memory of Prof. Dr. H. Wacker
Rights and permissions
About this article
Cite this article
Söderlind, G. Remarks on the stability of high-index DAEs with respect to parametric perturbations. Computing 49, 303–314 (1992). https://doi.org/10.1007/BF02238934
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF02238934