Abstract
Numerical solution of differential-algebraic equations with Hessenberg index-3 is considered by variational iteration method. We applied this method to two examples, and solutions have been compared with those obtained by exact solutions.
1. Introduction
Many important mathematical models can be expressed in terms of differential-algebraic equations (DAEs). Many physical problems are most easily initially modeled as a system of differential-algebraic equations (DAEs) [1]. Some numerical methods have been developed, using both BDF [1–3] and implicit Runge-Kutta methods [1], Padé and Chebysev approximations method [4–6]. These methods are only directly suitable for low-index problems and often require that the problem, have special structure. Although many important applications can be solved by these methods, there is a need for more general approaches. There are many new publication in the field of analytical sueveys such as [7–10]. The variational iteration method (VIM) was developed by He in [11]. The method is used by many researchers in a variety of scientific fields. The method has been proved by many authors [12–16] to be reliable and efficient for a variety of scientific applications, linear and nonlinear as well.
The most general form of a DAE is given by where may be singular. The rank and structure of this Jacobian matrix may depend, in general, on the solution , and for simplicity we will always assume that it is independent of . The important special case is of a semiexplicit DAE or an ODE with constraints:This is a special case of (1.1). The index is 1 if is nonsingular, because then one differentiation of (1.2b) yields in principle. For the semi-explicit index-1 DAE we can distinguish between differential variables and algebraic variables [1]. The algebraic variables may be less smooth than the differential variables by one derivative. In the general case, each component of may contain a mix of differential and algebraic components, which makes the numerical solution of such high-index problems much harder and riskier.
2. Special Differential-Algebraic Equations (DAEs) Forms
Most of the higher-index problems encountered in practice can be expressed as a combination of more restrictive structures of ODEs coupled with constraints. In such systems the algebraic and differential variables are explicitly identified for higher-index DAEs as well, and the algebraic variables may all be eliminated using the same number of differentiations. These are called Hessenberg forms of the DAE and are given below. In this paper, the variational iteration method has been proposed for solving differential-algebraic equations with Hessenberg index-3.
2.1. Hessenberg Index-1
One has Here the Jacobian matrix function is assumed to be nonsingular for all . This is also often referred to as a semi-explicit index-1 system. Semi-explicit index-1 DAEs are very closely related to implicit ODEs.
2.2. Hessenberg Index-2
One hasHere the product of Jacobians is nonsingular for all . Note the absence of the algebraic variables from the constraints (2.2b). This is a pure index-2 DAE, and all algebraic variables play the role of index-2 variables.
2.3. Hessenberg Index-3
One has Here the product of three matrix functions is nonsingular.
The index of a Hessenberg DAE is found, as in the general case, by differentiation. However, here only algebraic constraints must be differentiated.
3. He’s Variational Iteration Method (VIM)
Consider the differential equation where and are linear and nonlinear operators, respectively, and is the source inhomogeneous term. In [11], He proposed the variational iteration method where a correction functional for (3.1) can be written as where is a general Lagrange’s multiplier, which can be identified optimally via the variational theory and as a restricted variation which means . It is to be noted that the Lagrange multiplier can be a constant or a function.
The variational iteration method should be employed by following two essential steps. It is required first to determine the Lagrange multiplier that can be identified optimally via integration by parts and by using a restricted variation. Having determined, an iteration formula, without restricted variation, should be used for the determination of the successive approximations , , of the solution . The zeroth approximation can be any selective function. However, using the initial values , , and are preferably used for the selective zeroth approximation as will be seen later. Consequently, the solution is given by
3.1. First-Order ODEs
We first start our analysis by studying the first-order linear ODE of a standard form The VIM admits the use of the correction functional for this equation by where is the lagrange multiplier, that in this method may be a constant or a function, and is a restricted value where .
Taking the variation of both sides of (3.5) with respect to the independent variable we have that gives obtained upon using and . Integrating the integral of (3.6) by parts we obtain or equivalently The extremum condition of requires that . This means that the left hand side of (3.9) is 0, and as a result the right hand side should be 0 as well. This yields the stationary conditions This in turn gives Substituting this value of the lagrange multiplier into the functional (3.5) gives the iteration formula obtained upon deleting the restriction on that was used for the determination of . Considering the given condition , we can select the zeroth approximation . Using the selection into (3.12) we obtain the following successive approximations: Recall that that may give the exact solution if a closed form solution exists, or we can use the approximation for numerical purposes.
4. Applications
Example 4.1. We first considered the following differential-algebraic equations with Hessenberg index-3 form:
with initial conditions
The exact solutions are
where , represent the differential variables and represents the algebraic variables. After three times of differentiation of (4.1) we have the following ODE system:
Differential-algebraic equation (DAE) is a Hessenberg index-3 form.
To solve system (4.4), we can construct following correction functionals:
where , , and, are general Lagrange multipliers and , denote restricted variations, that is, .
Making the above correct functional stationary,
Its stationary conditions can be obtained as follows:
The Lagrange multipliers can be identified as follows:
and the following multipliers can be obtained as
Beginning with , , by the iteration formula (4.9), we have
Example 4.2. One has
with initial conditions
The exact solutions are , , :
where , represent the differential variables and represents the algebraic variables. After three times of differentiation of (4.11) we have the following ODE system:
Differential-algebraic equation (DAE) is a Hessenberg index-3 form.
To solve system (4.14), we can construct the following correction functionals:
By using the basic definition of the variational iteration method can obtain that
5. Conclusion
The method has been proposed for solving differential-algebraic equations with Hessenberg index-3. Results show the advantages of the method. Tables 1–4 and Figures 1–4 show that the numerical solution approximates the exact solution very well in accordance with the above method.
