Solution of nonlinear higher-index Hessenberg DAEs by Adomian polynomials and differential transform method

The solution of higher-index Hessenberg differential-algebraic equations (DAEs) is of great importance since this type of DAEs often arises in applications. Higher-index DAEs are known to be numerically and analytically difficult to solve. In this paper, we present a new analytical method for the solution of two classes of higher-index Hessenberg DAEs. The method is based on Adomian polynomials and the differential transform method (DTM). First, the DTM is applied to the DAE where the differential transforms of nonlinear terms are calculated using Adomian polynomials. Then, based on the index condition, the resulting recursion system is transformed into a nonsingular linear algebraic system. This system is then solved to obtain the coefficients of the power series solution. The main advantage of the proposed technique is that it does not require an index reduction nor a linearization. Two test problems are solved to demonstrate the effectiveness of the method. In addition, to extend the domain of convergence of the approximate series solution, we propose a post-treatment with Laplace-Padé resummation method.