| |
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 = 0of 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.
|