On the computation of impasse points of quasilinear differential-algebraic equations

Abstract

We present computational algorithms for the calculation of impasse points and higher-order singularities in quasi-linear 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.