Mathematical Problems in Engineering

Volume 2015, Article ID 717404, 7 pages

http://dx.doi.org/10.1155/2015/717404

## Power Series Extender Method for the Solution of Nonlinear Differential Equations

^{1}Electronic Instrumentation and Atmospheric Sciences School, Universidad Veracruzana, Cto. Gonzalo Aguirre Beltrán, S/N, 91000 Xalapa, VER, Mexico^{2}National Institute for Astrophysics, Optics and Electronics, Luis Enrique Erro No. 1, Santa Maria, 72840 Tonantzintla, PUE, Mexico

Received 1 October 2014; Accepted 1 December 2014

Academic Editor: Salvatore Alfonzetti

Copyright © 2015 Hector Vazquez-Leal and Arturo Sarmiento-Reyes. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We propose a power series extender method to obtain approximate solutions of nonlinear differential equations. In order to assess the benefits of this proposal, three nonlinear problems of different kind are solved and compared against the power series solution obtained using an approximative method. The problems are homogeneous Lane-Emden equation of index, governing equation of a burning iron particle, and an explicit differential-algebraic equation related to battery model simulations. The results show that PSEM generates highly accurate handy approximations requiring only a few steps. The main advantage of PSEM is to extend the domain of convergence of the power series solutions of approximative methods as Taylor series method, homotopy perturbation method, homotopy analysis method, variational iteration method, differential transform method, and Adomian decomposition method, among many others. From the application of PSEM, it results in handy easy computable expressions that extend the domain of convergence of high order power series solutions.

#### 1. Introduction

Solving nonlinear differential equations is an important task in sciences because many physical phenomena are modelled using such equations. Recently, the generalized homotopy method (GHM) [1] was proposed as a generalization of the homotopy perturbation method (HPM). The application of such method is based on a power series matching that enables GHM to obtain complex and rich expression impossible to obtain using HPM. Therefore, using as a guide the main idea of power series matching, behind GHM, we propose a power series method extender (PSEM). The key steps of this proposal are as follows.(1)First, we apply an approximative method to obtain the coefficients of a power series solution. The approximative method can be one from the literature like Taylor series method (TSM) [2, 3], power series method (PSM) [4], homotopy perturbation method (HPM) [5–11], perturbation method (PM), homotopy analysis method (HAM), variational iteration method (VIM), differential transform method (DTM), and Adomian decomposition method (ADM), among many others.(2)In the same fashion as GHM [1], we propose a trial function (TF) that can potentially describe the qualitative behaviour of the nonlinear problem.(3)Next, we apply the Taylor series method to a trial function (TF).(4)Then, we equate the coefficients of the power series obtained in steps 1 and 3 to obtain a nonlinear algebraic equation system in terms of the coefficients of the TF, which can be solved using symbolic or numerical methods.(5)Finally, the approximate solution is obtained by substituting the calculated coefficients from step 4 into the proposed TF.

In order to study the potential of the proposed technique, three nonlinear problems will be solved and compared versus numerical methods: homogeneous Lane-Emden equation of index [12], governing equation of burning iron particle [13], and an explicit differential-algebraic equation related to battery model simulations [14].

This paper is organized as follows. In Section 2, we introduce the basic concept of PSEM method. In Section 3, we show the approximation of three nonlinear differential equations related to different phenomenons in physics. Numerical simulations and a discussion about the results are provided in Section 4. Finally, a brief conclusion is given in Section 5.

#### 2. Basic Concept of PSEM Method

In a broad sense a nonlinear differential equation can be expressed as having as boundary condition where and are a linear operator and a nonlinear operator, respectively, is a known analytic function, is a boundary operator, is the boundary of domain , and denotes differentiation along the normal drawn outwards from [8].

Next, we express the solution of (1) as a power series where () are the coefficients of the power series.

It is important to notice that (3) can be obtained by some approximative method from literature as HPM, HAM, VIM, DTM, ADM, TSM, and PSM, among others.

Now, we propose that the solution for (1) can be written as a finite sum of functions in the general form [1] or where are constants to be determined by PSEM, are arbitrary trial functions, and and are the orders of approximations (4) and (5), respectively. We will denominate (4) and (5) as the trial function (TF).

Next, we calculate the Taylor series of (4) or (5), resulting in the power series respectively, where Taylor coefficients are expressed in terms of parameters .

Finally, we equate/match the coefficients of power series (6) or (7) with (3) to obtain the values of and substitute them into (4) or (5) to obtain the PSEM approximation.

It is important to notice that we can separately apply (4) or (5) to obtain an approximation of (1), where the selection of the TF depends of the behaviour of the problem under study. In addition, it is important to remark that if we choose the functions to be analytic, then (6) and (7) are convergent series [15–17].

#### 3. Case Studies

In the present section, we will solve three case studies to show the utility of the PSEM method to solve nonlinear differential equations. We know from literature that approximative methods as HPM, PM, HAM (using the as convergence control), VIM, DTM, ADM, TSM, and PSM are able to generate power series approximate solutions that in most of the cases are equivalent to the well-known TSM/PSM solutions. Therefore, the main difference in such cases is the difficulty of application of the specific approximative method to the particular case study. For its simplicity of application, we will use TSM [2] to exemplify the application of PSEM, although it is possible to use one of the available approximative methods from literature.

##### 3.1. Homogeneous Lane-Emden Equation of Index

The Lane-Emden singular equation describes a wide variety of problems in physics as some aspects related to the stellar structure, thermal history of spherical cloud of gas, isothermal gas spheres, and thermionic currents, among others [12]. The equation is expressed as where prime denotes differentiation with respect to and represents the index of the equation.

Using the initial conditions of (8) and considering as expansion point , it results the solution of (8) as where () are unknown constants to be determined by the Taylor series method.

In order to apply the TSM method [3], it is necessary to multiply (8) by to avoid the singularity, resulting in

Next, we derive successively (10) and resolve the system of equations obtained from the derivatives (), resulting in

Now, substituting the initial conditions of (8) into (11), it results in

From (5), we propose a specific solution for (8) with the following form: where () are constants to be determined and the order is chosen as 2 to obtain a handy approximation.

Next, Taylor series of (13) is

Then, we equate coefficients of -powers of (9) and (14) to obtain a system of equations, which can be solved symbolically, resulting in

Finally, the PSEM approximation is obtained by substituting (15) into (13), resulting in

##### 3.2. Combustion of a Single Iron Particle

The nonlinear energy equation for combustion of a single iron particle [13] is where prime denotes differentiation with respect to and is the dimensionless temperature.

The values of the parameters are obtained from [13]: , , , , and . Using the initial conditions of (17) and considering the expansion point yields to the following Taylor series: where () are unknown constants.

Now, we resolve (17) for and calculate the successive derivatives to obtain (), resulting in

Then, substituting the initial conditions of (17) into (19), we get

From (4), we propose a specific solution for (17) as follows:

Next, Taylor series of (21) is

Then, we equate coefficients of -powers of (18) and (22) to obtain a system of equations, which can be solved symbolically, resulting in where , , and .

Finally, substituting (23) into (21) yields to the PSEM approximation

For comparison purposes, we calculate the Taylor series by substituting and (20) into (18), yielding

##### 3.3. Differential-Algebraic Equation Related to Battery Model Simulation

Next, PSEM will be applied to a nonlinear equation [14] related to an oversimplified battery model, which is formulated as follows: where is the differential variable, is the algebraic variable, and primes denote derivative with respect to .

Using the initial conditions of (26) and considering the expansion point yields to the following Taylor series: where derivatives and () are unknown.

Firstly, we apply an implicit derivative with respect to of the second equation of (26), resulting in

As before, we calculate the successive derivatives of (28) and evaluate the resulting equations at , yielding where and .

From (4), we propose a specific solution for (26) as follows:

Finally, we calculate the fourth-order Taylor series of (30), equate the resulting coefficients of the same -powers with respect to (27), and solve the two systems of equations, resulting in the following PSEM approximation:

For comparison purposes, we calculate the Taylor series by substituting (29) into (27), yielding

#### 4. Numerical Simulation and Discussion

For all case studies, we used built-in numerical routines from Maple 15 for comparison purposes. The Fehlberg fourth-fifth order Runge-Kutta method with degree four interpolant (RKF45) was used [18, 19]. The command was setup with a tolerance of absolute error (A. E.) of .

We obtained a highly accurate rational approximate solution (16) for the singular second order Lane-Emden equation [12] (8) as depicted in Figures 1 and 2 for . Thus, the PSEM method may be useful for such kind of problems with singularities.