Differential algebraic equations DAEs

Most of the research on the analysis and numerical simulation of nonlinear DAEs has focused on systems in the fully implicit form of Equation (A.9). However, the generality of the form of the system in Equation (A.9) does not allow the development of explicit controller synthesis results. Also, the majority of chemical process applications (see examples throughout this book), as well as other engineering applications, are modeled by DAEs in a semi-explicit form, such that there is a distinct separation of the differential and algebraic equations  

The common underlying principle in the approaches for characterizing the solvability of a DAE system is to obtain, either explicitly, or implicitly, a local representation of an equivalent ODE system, for which available results on existence and uniqueness of solutions are applicable. The derivation of the underlying ODE system involves the repeated differentiation of the algebraic constraints of the DAE, and it is this differentiation process that leads to the concept of a DAE index that is widely used in the literature. For the semi-explicit DAE systems (A. 10) that are of interest to us here, the index has the following definition. 

Definition A.5. The index yj of the DAE system in Equation (A.10) with specified smooth inputs u(t) is the minimum number of times the algebraic equations or their subset have to be differentiated to obtain a set of differential equations for z, i.e., in order to be able to solve z = iF(x, z,t) for z. 

The index yj provides a measure of the singularity of the algebraic equations and the resulting differences from ODE systems. More specifically consider the DAE system of Equation (A.10) in the case in which the matrix L(x) 

Definition A.6. A differential algebraic equation system (A. 10) is said to be regular, if 

---

Off-line analysis, controller design, and optimization are now performed in the area of dynamics. The largest dynamic simulation has been about 100,000 differential algebraic equations (DAEs) for analysis of control systems. Simulations formulated with process models having over 10,000 DAEs are considered frequently. Also, detailed training simulators have models with over 10,000 DAEs. On-line model predictive control (MPC) and nonlinear MPC using first-principle models are seeing a number of industrial applications, particularly in polymeric reactions and processes. At this point, systems with over 100 DAEs have been implemented for on-line dynamic optimization and control. 

The application of simultaneous optimization to reactor-based flowsheets leads us to consider the more general problem of differentiable/algebraic optimization problems. Again, the optimization problem needs to be reconsidered and reformulated to allow the application of efficient nonlinear programming algorithms. As with flowsheet optimization, older conventional approaches require the repeated execution of the differential/algebraic equation (DAE) model. Instead, we briefly describe these conventional methods and then consider the application and advantages of a simultaneous approach. Here, similar benefits are realized with these problems as with flowsheet optimization. 

An alternative to the standard-form representation is the differential-algebraic equation (DAE) representation, which is stated in a general form as g(r, y, y). The lower portion of Fig. 7.3 illustrates how the heat equation is cast into the DAE form. The boundary conditions can now appear as algebraic constraints (i.e., they have no time derivatives). For a problems as simple as the heat equations, this residual representation of the boundary conditions is not necessary. However, recall that implicit boundary-condition specification is an important aspect of solving boundary-layer equations. 

There are many high-quality, well-documented, software packages available to solve stiff problems in this form. However, one often encounters chemically reacting flow problems that are not easily posed as standard-form ODEs. In these cases problems can often be posed easily in a more general form, called differential-algebraic equations (DAE),... 

Taken together, the system of equations represents a set of stiff ordinary differential equations, which can be solved numerically. Because more than one dependent-variable derivative can appear in a single equation (e.g., the momentum equation has velocity and pressure derivatives), it is usually more convenient to use differential-algebraic equation (DAE) software (e.g., Dassl) for the solution rather than standard-form ODE software. 

A brief explanation of differential-algebraic equations (DAE) facilitates a further mathematical discussion of the stagnation-flow equations. In general, DAEs are stated as a vector residual equation, where w is the dependent-variable vector and the prime denotes a time derivative. For the discussion here, it is convenient to consider a restricted class of DAEs called semi-explicit nonlinear DAEs, which are represented as... 

Dassl, solves stiff systems of differential-algebraic equations (DAE) using backward differentiation techniques [13,46]. The solution of coupled parabolic partial differential equations (PDE) by techniques like the method of lines is often formulated as a system of DAEs. It automatically controls integration errors and stability by varying time steps and method order. 

MESH) equations which are solved for the whole column, decanter included and taking into account the liquid-liquid phase split. Numerical treatment of the Differential Algebraic Equation (DAE) system and discrete events handling is performed with DISCo, a numerical package for hybrid systems with a DAE solver based on Gear s method. The column technical features and operating conditions are shown in Table 4. A sequence of two operational batch steps, namely... 

The resulting system is called a set of differential-algebraic equations (DAE) and their solution is now a specialised field with its own texts [130, 286] and there is a package program, DASSL [441], for their solution. This can be of use in the present context, for example with the method of lines, which indeed often results in a DAE system. This is gone into in some detail in Chap. 9, in the context of Rosenbrock methods. 

Method of Lines (MOL) and Differential Algebraic Equations (DAE) 165 Method of Wu and White... 

The dynamic models of chemical processes are represented by differential-algebraic equations (DAEs). Equation (2) and (3) define such a system. Equations (4), (5) and (6) are the path constraints on the state variables, control variables and algebraic variables respectively, while equation (7) represents the initial condition of the state variables. Obj is a scalar objective function at final time, tj. ... 

The mathematical model of a MAT reactor, considering a 12 lump model, has been discretized using a finite element method in the direction of gas flow. The resulting system of differential-algebraic equations (DAEs) has been solved by an appropriate computer code (DASSL). 

Several of these simple mass balances with basic rate expressions were solved analytically. In the case of multiple reactions with nonlinear rate expressions (i.e., not first-order reaction rates), the balances must be solved numerically. A high-quality ordinary differential equation (ODE) solver is indispensable for solving these problems. For a complex equation of state and nonconstant-volume case, a differential-algebraic equation (DAE) solver may be convenient. 

Some reactor models require a more general structure than the ODE, dxjdt = fix, t). The nonconstant density, nonconstant volume semibatch and CSTR reactors in Chapter 4 are more conveniently expressed as differential-algebraic equations (DAEs). To address these models, consider the more general form of implicit ODEs... 

The calculations for the nonconstant-density case may be greatly simplified by using a differential-algebraic equation (DAE) solver, All three cases enumerated above can be handled by modifying the residual equations provided to the DAE solver. We do not have to differentiate the equation of state or perform other algebraic manipulations that are required if one uses an ordinary differential equation (ODE) solver. 

In Equation-Oriented (EO) approach all the modelling equations are assembled in a large sparse system producing Non-linear Algebraic Equations (NAE) in steady state simulation, and stiff Differential Algebraic Equations (DAE) in dynamic simulation. Thus, the solution is obtained by solving simultaneously all the modelling equations. 

In Equation-Oriented (EO) approach the software architecture is close to a solver of equations. EO is more suited for dynamic simulation since this can be modelled by a system of differential-algebraic equations (DAE) of the form ... 

It should be noted that (MIP) problems, and their special cases, may be regarded as steady-state models. Hence, one important extension is the case of dynamic models, which in the case of discrete time models gives rise to multiperiod optimization problems, while for the case of continuous time it gives rise to optimal control problems that contain differential-algebraic equation (DAE) models. 

Large-order chemical engineering problems, when solved numerically, usually give rise to a set of coupled differential-algebraic equations (DAE). This type of coupling is more difficidt to deal with compared to coupled algebraic equations or coupled ODEs. The interested reader should refer to Brenan et al. (1989) for the exposition of methods for solving DAE. 

Primitive model is a single Differential-Algebraic Equation (DAE) derived from basic conservation laws in different domains describing the dynamic response of a constiment element (e.g., capacitors, inductors, resistors, and channels) in a complex system. 

The formulation of a mathematical model becomes mode dependent. In one mode it may be an explicit state space model. If storage elements become dependent the model turns into a set of differential-algebraic equations (DAEs). 

---