On the Computation of Impasse Points of Quasi-Linear Differential-Algebraic Equations

Abstract

Abstract : We present computational algorithms for the calculation of impasse points and higher order singularities in quasilinear differential-algebraic equations. Our method combines a reduction step transforming the DAE into a singular ODE with an augmentation procedure inspired by numerical bifurcation theory. Singularities are characterized by the vanishing of a scalar quantity that may be monitored along any trajectory. Two numerical examples with physical relevance are given.