SciCADE99
 

Abstract

    

Computation of Consistent Initial Values for Nonlinear Lower Index DAEs

René Lamour
lamour@mathematik.hu-berlin.de
Humboldt-University, Germany

The solution of initial value problems for differential-algebraic equations (DAE)

f(x',x,t) = 0
(1)
PP1(x(t0),t0)(x(t0)-a) = 0
of lower index ( up to 2 ) requires the computation of consistent initial values e.g. to start an integration numerically. Using a suitable index reduction method, the computation of consistent initial values in t0 leads to the system

f(y0,x0,t0) = 0
(2)
PP1(x0,t0)(x0 - a)+ PQ1(x0,t0)(y0 +A2-1f't(y0,x0,t0))+Qy0 = 0

of nonlinear algebraic equations for determining the unknowns y0,x0. For quasilinear DAEs of the form

A(x,t) x' + b(x,t) = 0
we can prove, under mild additional assumptions, that (2) has a nonsingular Jacobian matrix. Examples from electrical network simulation show the applicability of the method.

MINISIMPOSIUM SESSION: 13. Numerical methods for DAEs

Submitted: 18/Feb/99
[SciCADE99 | Abstracts | Sessions]