An efficient asymptotically correct error estimator for collocation solutions to singular index-1 DAEs

Abstract

A computationally efficient a posteriori error estimator is introduced and analyzed for collocation solutions to linear index-1 DAEs (differential-algebraic equations) with properly stated leading term exhibiting a singularity of the first kind. The procedure is based on a modified defect correction principle, extending an established technique from the context of ordinary differential equations to the differentialalgebraic case. Using recent convergence results for stiffly accurate collocation methods, we prove that the resulting error estimate is asymptotically correct. Numerical examples demonstrate the performance of this approach. To keep the presentation reasonably self-contained, some arguments from the literature on DAEs concerning the decoupling of the problem and its discretization, which is essential for our analysis, are also briefly reviewed. The appendix contains a remark about the interrelation between collocation and implicit Runge-Kutta methods for the DAE case.