Differential-algebraic system of equations

This article's lead section may be too long. Please read the length guidelines and help move details into the article's body. (March 2023)

In mathematics, a differential-algebraic system of equations (DAE) is a system of equations that either contains differential equations and algebraic equations, or is equivalent to such a system.

The set of the solutions of such a system is a differential algebraic variety, and corresponds to an ideal in a differential algebra of differential polynomials.

In the univariate case, a DAE in the variable t can be written as a single equation of the form

${\displaystyle F(x'(t),x(t),t)=0,}$

where ${\displaystyle x(t)}$ is a vector of unknown functions.

They are distinct from ordinary differential equation (ODE) in that a DAE is not completely solvable for the derivatives of all components of the function x because these may not all appear (i.e. some equations are algebraic); technically the distinction between an implicit ODE system [that may be rendered explicit] and a DAE system is that the Jacobian matrix ${\displaystyle {\frac {\partial F(u,v,t)}{\partial u}}}$ is a singular matrix for a DAE system.^[1] This distinction between ODEs and DAEs is made because DAEs have different characteristics and are generally more difficult to solve.^[2]

In practical terms, the distinction between DAEs and ODEs is often that the solution of a DAE system depends on the derivatives of the input signal and not just the signal itself as in the case of ODEs;^[3] this issue is commonly encountered in nonlinear systems with hysteresis,^[4] such as the Schmitt trigger.^[5]

This difference is more clearly visible if the system may be rewritten so that instead of x we consider a pair ${\displaystyle (x,y)}$ of vectors of dependent variables and the DAE has the form

${\displaystyle {\begin{aligned}{\dot {x}}(t)&=f(x(t),y(t),t),\\0&=g(x(t),y(t),t).\end{aligned}}}$

where ${\displaystyle x(t)\in \mathbb {R} ^{n}}$, ${\displaystyle y(t)\in \mathbb {R} ^{m}}$, ${\displaystyle f:\mathbb {R} ^{n+m+1}\to \mathbb {R} ^{n}}$ and ${\displaystyle g:\mathbb {R} ^{n+m+1}\to \mathbb {R} ^{m}.}$

A DAE system of this form is called semi-explicit.^[1] Every solution of the second half g of the equation defines a unique direction for x via the first half f of the equations, while the direction for y is arbitrary. But not every point (x,y,t) is a solution of g. The variables in x and the first half f of the equations get the attribute differential. The components of y and the second half g of the equations are called the algebraic variables or equations of the system. [The term algebraic in the context of DAEs only means free of derivatives and is not related to (abstract) algebra.]

The solution of a DAE consists of two parts, first the search for consistent initial values and second the computation of a trajectory. To find consistent initial values it is often necessary to consider the derivatives of some of the component functions of the DAE. The highest order of a derivative that is necessary for this process is called the differentiation index. The equations derived in computing the index and consistent initial values may also be of use in the computation of the trajectory. A semi-explicit DAE system can be converted to an implicit one by decreasing the differentiation index by one, and vice versa.^[6]