Abstract

We implement a relatively new analytic iterative technique to get approximate solutions of differential algebraic equations system based on generalized Taylor series formula. The solution methodology is based on generating the residual power series expansion solution in the form of a rapidly convergent series with easily computable components. The residual power series method (RPSM) can be used as an alternative scheme to obtain analytical approximate solution of different types of differential algebraic equations system applied in mathematics. Simulations and test problems were analyzed to demonstrate the procedure and confirm the performance of the proposed method, as well as to show its potentiality, generality, viability, and simplicity. The results reveal that the proposed method is very effective, straightforward, and convenient for solving different forms of such systems.

1. Introduction

Differential algebraic equations (DAEs) have gained in the recent years considerable importance and popularity partly due to its powerful potential applications. The applications of the DAEs can be found in many fields of scientific and engineering including electrical networks, circuit analysis, control theory, robotics, power systems, and chemical process simulations [15], which can be successfully modeled by coupled differential and algebraic equation. In this regard, the DAEs are type of differential equations, in which the unknown functions are satisfying additional algebraic equations, whereas the derivatives are not in general expressed explicitly and typically derivatives of some of the dependent variables may not appear in the equations at all. In other words, the DAEs consist of a set of differential equations together with additional algebraic constraints.

Many mathematicians, engineers, and physicists have devoted considerable effort to studying the nature of the DAEs involving explicit and approximate solutions where a complicated system of DAEs can be described quite precisely with mathematical expressions. Generally, many explicit solutions have been found to the linear DAEs, but there exists no method that yields an explicit solution for the nonlinear DAEs. Most often, these systems are typically modeled as complex nonlinear DAEs. Thus, a purely mathematical analysis is often not possible to obtain a closed-form solution. Approximations based on numerical methods are therefore needed.

In science literature, a lot of studies about the DAEs have appeared according to [69] and references therein. In this point, it is usually difficult to solve these kinds of DAEs analytically. Therefore, there are many powerful numerical methods in literature that can be used to approximate solutions to the DAEs system. To mention but a few, the backward differentiation formula (BDF) has been the first numerical method that was used to obtain the solutions for a class of DAEs [10]. In [11], the implicit Runge-Kutta method has been applied to solve the DAEs numerically. Çelik et al., in [12], have presented approximate solutions for chemical DAEs using the Padé series method. Also, Çelik and Yeloglu [13] have used the Chebyshev series approximation for solving other kinds of DAEs. In contrast, Hosseini has employed the Adomian decomposition method (ADM) for solving linear and nonlinear forms of the DAEs [14, 15]. In addition, the variational iteration method (VIM) [16] and the homotopy perturbation method (HPM) [17] have been introduced for dealing with different kinds of DAEs and their engineering applications. Lately, the authors [18] have developed a new modified method to further investigate for solving linear and nonlinear system of DAEs of higher index without index reduction based on the HPM. However, none of previous studies propose a methodical way to solve systems of DAEs. Moreover, previous studies require more effort to achieve the results and usually they are suited for linear form of DAEs system. The applications of other versions for linear and nonlinear problems can be found in [1926] and references therein.

Furthermore, the RPSM has been developed as an efficient numerical as well as analytical method to determine the coefficients of power series solutions for a class of fuzzy differential equation by Abu Arqub [27]. Besides, the RPSM has been successfully applied to get numerical solutions for many other problems, for instance, generalized Lane-Emden equation which is a highly nonlinear singular differential equation [28], regular initial value problems [29], composite and noncomposite fractional differential equations [30], and boundary value problems of fractional order [31]. This method is effective to construct power series expansion solution for strongly linear and nonlinear equations without linearization, perturbation, or discretization [32], which computes the coefficients of the power series solutions by chain of linear equations of one variable. The RPSM is an alternative procedure for obtaining analytical Taylor series solution for system of DAEs. Consequently, using the residual error concept, we get a series solution in practice as well a truncated series solution. For linear problems, the exact solution can be obtained by few terms of the RPS solution.

The RPSM is an approximate analytical method for solving large amount of problems of ordinary, partial, integral, and integrodifferential equations. The biggest advantage over the existing standard analytical techniques is that it overcomes the difficulty arising in calculating complicated terms. However, the method provides the solution in terms of a rapidly convergent series with easily computable components. The analytical approximate solution should be constructed in the form of a polynomial which does not exhibit the real behaviors of the problem but gives a good approximation to the true solution in the given interval. Further, this approach is different from the traditional higher order Taylor series method. The Taylor series method is computationally expensive for large orders and is suited for the linear problems, whilst the RPS techniques need small computational requirements together with high precision.

In this paper, we consider the following system of differential algebraic equations: subject to the initial conditions where , , , and are real finite constants with , are nonlinear continuous functions of , , , , which satisfy all necessary requirements of the existence of a unique solution, and are unknown functions of independent variable to be determined. Throughout this paper, we assume that , are analytical functions on the given interval.

For several reasons, we consider (1) directly rather than to try to reformulate it as an initial value problem of the form Some often, via simulation of physical phenomena as a matter of fact, the model takes form of system of DAEs depicting a collection of relationships between the variables with interest and some of their derivatives. These relationships may even be generated automatically by modeling or simulation program. The variables usually have a physical significance; changing the model to (2) may produce less meaningful variables. In the case of computer-generated or nonlinear models, especially complex (nonlinear) models, it may be time consuming or impossible to obtain an explicit model. Parameters are present in many applications where any change in the parameter values can alter the relationships between variables, which require different explicit models with solution manifolds of different dimensions. On the other hand, the index of DAEs systems is a measure of the degree of singularity of the system and also widely regarded as an indication of certain difficulties for numerical methods. However, if the original DAEs can be solved directly, then it becomes easier to explore the effect of modeling changes and parameter variation and also to interface the modeling software directly with design software. These advantages enable researchers to focus their attention on the physical phenomena of interest. The desirability of working directly with DAEs has been recognized for over twenty years by scientists and engineers in several areas. Depending on the area, DAEs have been called singular, implicit, differential-algebraic, descriptor, generalized state space, noncanonic, noncausal, degenerate, semistate, constrained, reduced order model, and nonstandard systems.

The basic motivation of this paper is to apply the RPSM to develop an approach for obtaining the representation of exact and approximate solutions to strongly linear and nonlinear system of DAEs. This approach is simple and needs less effort to achieve the results. It does not require any converting while switching from first to higher order in which the solutions and all its derivatives are applicable for each arbitrary point in the given interval. Thus, the RPSM can be applied directly to the given problems by choosing an appropriate value for the initial guess approximation.

The remainder of this paper is organized as follows. In Section 2, basic idea of the residual power series method (RPSM) together with analysis of the method is presented. In the meantime, the RPSM is extended to provide symbolic approximate series solutions for system of differential algebraic equations (1) and (2). Based on the above, numerical examples are given to illustrate the capability of the proposed method in Section 3. Results reveal that only few terms are required to deduce the approximate solutions which are found to be accurate and efficient. Finally, some conclusions are summarized in the last section.

2. Adaptation of Residual Power Series Method (RPSM)

In this section, we present a brief description and preliminary results of the standard RPSM that will be used in the remainder of this paper. In the meantime, we employ our technique for finding out a series solution of system (1) and (2) by formulating and analyzing the method in solving such systems. Afterward, a convergence theorem is presented in order to capture the behavior of the solution.

The RPSM consists of expressing the solutions of system (1) and (2) as a power series expansion about the initial point . To achieve our goal, we suppose that these solutions take the following form: where are the terms of approximations such that .

For , substituting the initial guesses , which are known from initial conditions (2), into , , leads to the approximate solutions for the DAEs system such that , , whereas , for , can be calculated by the following th-truncated series:

Regarding applying the RPSM, the th-residual functions and the th residual functions are given, respectively, by

Obviously, for each , which are infinitely differentiable functions at . Furthermore, , , and these relations are fundamental rules in the RPSM and its applications. Particularly, , , .

Now, substituting the th-truncated series into (6) yields

Consequently, based on , , setting and in (8) leads to which contained equations in variables that can be solved directly by using the following Mathcad commands:

Thus, by using (5), the first approximate solutions of system (1) and (2) can be written as follows:

Next, in order to obtain the second approximate solutions, we set such that . Then, we differentiate both sides of (8) with respect to and substitute to obtain

According to the fact that , , the values of can be found easily by using the following Mathcad commands:

Hence, the second approximate solutions of system (1) and (2) can be written as follows:

By the same technique, the process can be repeated to generate a sequence of approximate solutions for the system (1) and (2). Moreover, higher accuracy can be achieved by evaluating more components of the solution.

It will be convenient to have a notation for the error in the approximation . Accordingly, we will let denote the difference between and its th Taylor polynomial; that is, where the functions are called the th remainder for the Taylor series of . In fact, it often happens that the remainder become smaller and smaller, approaching zero, as gets large.

Next, we present a convergence theorem of the RPSM to capture the behavior of these solutions. Afterwards, we introduce the error functions to study the accuracy and efficiency of the method. Actually, continuous approximations to the solution will be obtained. Taylor’s theorem allows us to represent fairly general functions exactly in terms of polynomials with a known, specified, and bounded error. The next theorem will guarantee convergence to the exact analytic solution of (1) and (2).

Theorem 1. Suppose that , , is the exact solution for system (1) and (2). Then, the approximate solution obtained by the RPSM is in fact the Taylor expansion of .

Proof. The proof of Theorem 1 is similar to proof of Theorem 1 in [28].

Corollary 2. Let , be a polynomial for some ; then the RPSM will obtain the exact solution.

To show the accuracy of the present method for our problems, we report four types of error. The first one is the residual error, , and is defined as while the exact, , relative, , and consecutive,, errors are defined, respectively, by for , where , are the th-order approximation of obtained by the RPSM and are the exact solution. An excellent account of the study of error analysis, which includes its definitions, varieties, applications, and method of derivations, can be found in [34].

3. Simulations and Test Problems

To give a clear overview of the content of this work, we consider some examples to demonstrate the performance and efficiency of the present technique. The proposed method provides an analytical approximate solution in terms of an infinite power series. However, there is a practical need to evaluate this solution and to obtain numeric values from the infinite power series. The consequent series truncation and the practical procedure are conducted to accomplish this task and transform the otherwise analytical results into an exact solution, which is evaluated to a finite degree of accuracy. In most real life situations, the differential equation that models the problem is too complicated to solve exactly, and there is a practical need to approximate the solution.

Example 1. Consider the following linear system of DAEs [33]: subject to the initial conditions The exact solution is , , .
To apply the RPS approach in solving system (18) and (19), we start with selecting initial guesses of the approximations such that , , and ; then the th-truncated series solutions , , for system (18) and (19) are given by
Accordingly, the unknown coefficients , , , can be found by constructing the following th residual functions , , such that
Now, in order to find the 1st-approximate solutions, put through (21) to get that
Using the facts that , the values of , , are given by , , and . Based on this, the first approximations , , of the RPS solution for (18) and (19) can be written by ,  , and . In contrast, the second approximations , , of the RPS solutions have the forms ,   , and , in which the coefficients , , and can be found by putting through (21) and differentiating both sides of (21) with respect to such that
Consequently, , , and as soon as . Therefore, the second approximations for (18) and (19) are ,  , and . On the other hand, the fourth approximations of the RPS solution for (18) and (19) according to of the th residual function (21) are given by
Likewise, the 10th truncated series , , of the RPS solution are given by
Thus, the approximate solutions , , of system (18) and (19) have general forms that coincide with the exact solutions such that
To show the accuracy of the method, numerical results at some selected grid points together with comparison between the absolute errors of RPSM for 8th-order approximations and the ADM [33] are given in Tables 1 and 2. From the tables, it can be seen that the present method provides us with an accurate approximate solution to system (18) and (19). Indeed, the results reported in these tables confirm the effectiveness of the RPS method.

Example 2. Consider the following nonlinear system of DAEs [18]: for , subject to initial conditions The exact solution is , , and .
As we mentioned earlier, if we select the first terms as initial guesses of the approximations as , , and , then the th-truncated series solutions , , for (27) and (28) have the following form: in which the values of the coefficients , , , can be found by constructing the th residual functions , , as follows:
Using the facts and through (30), we can obtain ,   and as soon as , , , , , and . Furthermore, based on , , , the 10th truncated series , , of the RPS solution for system (27) and (28) are given by
Therefore, the approximate solutions of system (27) and (28) can be expressed as which are compatible with the exact solutions.
In order to show how the th values affect the approximate solutions, we calculate the maximum error functions of , ,  , for various values of in with step size of 0.1 that are listed in Table 3 by using the th order RPS approximate solution. In Table 4, the residual error functions , , , have been listed for , , to demonstrate the rapid convergence of the present method by increasing the order of RPS approximation. As a result, it is clear from these tables that the approximate solutions are found to be in good agreement with the exact solutions for all values of in .

Example 3. Consider the following nonlinear system of DAEs: subject to the initial conditions The exact solution is , , .
To apply the RPS approach for solving system (33) and (34), if we start with selecting the initial guesses of the approximations such that , , and , then the th-truncated series solutions , , can be given by
Accordingly, the unknown coefficients , , , can be found by constructing the following th residual functions , , such that
Now, in order to find the 1st-approximate solutions, put through (36) to get that
Use the fact that , to get , , and . Based upon this, the first approximations , , of the RPS solution can be written by ,  , and . Furthermore, the second approximations , , of the RPS solution for system (33) and (34) are given by where the coefficients , , and will be found by putting through (36) and using the facts that such that , , and . Therefore, the second approximations for the system are ,  , and . By continuing with the similar fashion, the 10th-order approximations , , of the RPS solution lead to the following results:
Correspondingly, the general forms of the solutions ,  , and for (33) and (34) are given by which are coinciding with the exact solutions.
To illustrate the convergence of the approximate solutions , , to the exact solutions with respect to the th-order of solutions, we present numerical results of Example 3 graphically by Figures 1, 2, and 3, which show the exact solution , , and some iterated approximations , , , respectively. These graphs reveal that the proposed method is an effective and convenient method for solving nonlinear DAEs systems with less computational and iteration steps.

4. Conclusion

In this paper, a new efficient and accurate numerical technique is proposed and applied to handle the linear and nonlinear system of DAEs based on the use of the residual power series method (RPSM). The present technique is performed based mainly on the generation of residual error function and then applying the generalized Taylor series formula. The numerical solutions obtained are an infinite series for appropriate initial conditions. In addition, the RPSM does not require linearization, perturbation, or discretization of the variables, as well as avoiding the round-off errors. The results reveal that the RPSM is a powerful and straightforward tool and is convenient to handle a various range of engineering problems. The approximate solutions in this analysis are found in the closed form of a convergent series.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors would like to express their thanks to the unknown referees for their careful reading and helpful comments.