A Successive Linear Programming Approach for Initialization and Reinitialization after Discontinuities of Differential-Algebraic Equations
- 1 January 1998
- journal article
- Published by Society for Industrial & Applied Mathematics (SIAM) in SIAM Journal on Scientific Computing
- Vol. 20 (2) , 447-467
- https://doi.org/10.1137/s1064827596307725
Abstract
Determination of consistent initial conditions is an important aspect of the solution of differential-algebraic equations (DAEs). Specification of inconsistent initial conditions, even if they are only slightly inconsistent, often leads to a failure in the initialization problem. In this paper, we present a successive linear programming (SLP) approach for the solution of the DAE derivative array equations for the initialization problem. The SLP formulation handles roundoff errors and inconsistent user specifications, among other things, and allows for reliable convergence strategies that incorporate variable bounds and trust region concepts. A new consistent set of initial conditions is obtained by minimizing the deviation of the variable values from the specified ones. For problems with discontinuities caused by a step change in the input functions, a new criterion is presented for identifying the subset of variables which are continuous across the discontinuity. The SLP formulation is then applied to determine a consistent set of initial conditions for further solution of the problem in the domain after the discontinuity. Numerous example problems are solved to illustrate these concepts.Keywords
This publication has 15 references indexed in Scilit:
- Nonsmooth dynamic simulation with linear programming based methodsComputers & Chemical Engineering, 1997
- A class of smoothing functions for nonlinear and mixed complementarity problemsComputational Optimization and Applications, 1996
- Reinitialization of DAE's after discontinuitiesComputers & Chemical Engineering, 1995
- Progress on a general numerical method for nonlinear higher index DAEs IICircuits, Systems, and Signal Processing, 1994
- Approximation Methods for the Consistent Initialization of Differential-Algebraic EquationsSIAM Journal on Numerical Analysis, 1991
- On methods for reducing the index of differential algebraic equationsComputers & Chemical Engineering, 1990
- A proposed numerical algorithm for solving nonlinear index problemsIndustrial & Engineering Chemistry Research, 1990
- Differential-Algebraic Equation Index TransformationsSIAM Journal on Scientific and Statistical Computing, 1988
- The Use of Linear Programming for the Solution of Sparse Sets of Nonlinear EquationsSIAM Journal on Scientific and Statistical Computing, 1987
- Sensitivity analysis of initial value problems with mixed odes and algebraic equationsComputers & Chemical Engineering, 1985