Improving the accuracy of BDF methods for index 3 dierential algebraic equations

• Methods for solving index 3 DAEs based on BDFs suffer a loss of accuracy when there is a change of step size or a change of order of the method. A layer of nonuniform convergence is observed in these cases, and O(1) errors may appear in the algebraic variables. From the viewpoint of error control, it is beneficial to allow smooth changes of step size, and since most codes based on BDFs are of variable order, it is also of interest to avoid the inaccuracies caused by a change of order of the method. In the case of BDFs applied to index 3 DAEs in semi-explicit form, we present algorithms that correct to O(h) the inaccurate approximations to the algebraic variables when there are changes of step size in the backward Euler method. These algorithms can be included in an existing code at a very small cost. We have also described how to obtain formulas that correct the O(1) errors in the algebraic variables appearing after a change of order.

Subject headings

NATURVETENSKAP  -- Matematik -- Beräkningsmatematik (hsv//swe)

NATURAL SCIENCES  -- Mathematics -- Computational Mathematics (hsv//eng)

Keyword

BDF

differential-algebraic equation

index 3

Publication and Content Type

art (subject category)

ref (subject category)