The Open UniversitySkip to content
 

Multiple shooting using a dichotomically stable integrator for solving differential-algebraic equations

England, Roland; Lamour, René and López-Estrada, Jesús (2002). Multiple shooting using a dichotomically stable integrator for solving differential-algebraic equations. Applied Numerical Mathematics, 42(1-3) pp. 117–131.

DOI (Digital Object Identifier) Link: http://dx.doi.org/10.1016/S0168-9274(01)00145-3
Google Scholar: Look up in Google Scholar

Abstract

In previous work by the first author, it has been established that a dichotomically stable discretization is needed when solving a stiff boundary-value problem in ordinary differential equations (ODEs), when sharp boundary layers may occur at each end of the interval. A dichotomically stable implicit Runge–Kutta method, using the 3-stage, fourth-order, Lobatto IIIA formulae, has been implemented in a variable step-size initial-value integrator, which could be used in a multiple-shooting approach.

In the case of index-one differential–algebraic equations (DAEs) the use of the Lobatto IIIA formulae has an advantage, over a comparable Gaussian method, that the order is the same for both differential and algebraic variables, and there is no need to treat them separately.

The ODE integrator (SYMIRK [R. England, R.M.M. Mattheij, in: Lecture Notes in Math., Vol. 1230, Springer, 1986, pp. 221–234]) has been adapted for the solution of index-one DAEs, and the resulting integrator (SYMDAE) has been inserted into the multiple-shooting code (MSHDAE) previously developed by R. Lamour for differential–algebraic boundary-value problems. The standard version of MSHDAE uses a BDF integrator, which is not dichotomically stable, and for some stiff test problems this fails to integrate across the interval of interest, while the dichotomically stable integrator SYMDAE encounters no difficulty. Indeed, for such problems, the modified version of MSHDAE produces an accurate solution, and within limits imposed by computer word length, the efficiency of the solution process improves with increasing stiffness. For some nonstiff problems, the solution is also entirely satisfactory.

Item Type: Journal Article
Copyright Holders: 2001 IMACS
ISSN: 1873-5460
Extra Information: Numerical Solution of Differential and Differential-Algebraic Equations, 4-9 September 2000, Halle, Germany
Keywords: differential–algebraic equations; boundary-value problems; dichotomic stability; multiple shooting
Academic Unit/Department: Mathematics, Computing and Technology > Mathematics and Statistics
Item ID: 2403
Depositing User: Roland England
Date Deposited: 12 Jun 2006
Last Modified: 22 May 2014 10:40
URI: http://oro.open.ac.uk/id/eprint/2403
Share this page:

Actions (login may be required)

View Item
Report issue / request change

Policies | Disclaimer

© The Open University   + 44 (0)870 333 4340   general-enquiries@open.ac.uk