Research

# Analysis of the new homotopy perturbation method for linear and nonlinear problems

Ali Demir1*, Sertaç Erman1, Berrak Özgür1 and Esra Korkmaz2

Author Affiliations

1 Department of Mathematics, Kocaeli University Umuttepe, Izmit, Kocaeli, 41380, Turkey

2 Ardahan University, Merkez, Ardahan, 75000, Turkey

For all author emails, please log on.

Boundary Value Problems 2013, 2013:61  doi:10.1186/1687-2770-2013-61

The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2013/1/61

 Received: 14 December 2012 Accepted: 6 March 2013 Published: 26 March 2013

© 2013 Demir et al.; licensee Springer

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

### Abstract

In this article, a new homotopy technique is presented for the mathematical analysis of finding the solution of a first-order inhomogeneous partial differential equation (PDE) u x ( x , y ) + a ( x , y ) u y ( x , y ) + b ( x , y ) g ( u ) = f ( x , y ) . The homotopy perturbation method (HPM) and the decomposition of a source function are used together to develop this new technique. The homotopy constructed in this technique is based on the decomposition of a source function. Various decompositions of source functions lead to various homotopies. Using the fact that the decomposition of a source function affects the convergence of a solution leads us to development of a new method for the decomposition of a source function to accelerate the convergence of a solution. The purpose of this study is to show that constructing the proper homotopy by decomposing f ( x , y ) in a correct way determines the solution with less computational work than using the existing approach while supplying quantitatively reliable results. Moreover, this method can be generalized to all inhomogeneous PDE problems.

### 1 Introduction

Many problems in the fields of physics, engineering and biology are modeled by linear and nonlinear PDEs. In recent years, the homotopy perturbation method (HPM) has been employed to solve these PDEs [1-6]. HPM has been used extensively to solve nonlinear boundary and initial value problems. Therefore HPM is of great interest to many researchers and scientists. HPM, first presented by Ji Huan He [7-9], is a powerful mathematical tool to investigate a wide variety of problems arising in different fields. It is obtained by successfully coupling homotopy theory in topology with perturbation theory. In HPM, a complicated problem under study is continuously deformed into a simple problem which is easy to solve to obtain an analytic or approximate solution [10].

In this work, a new homotopy perturbation technique is proposed to find the solution u ( x , y ) based on the decomposition of a right-hand side function f ( x , y ) which leads to the construction of new homotopies. It is shown that decomposing f ( x , y ) in a correct way, requiring less computational work than the existing approach, helps us to determine the solution. The results reveal that the proposed method is very effective and simple.

The paper is organized as follows. In Section 2, an analysis of the new homotopy perturbation method is given. The construction of a new homotopy based on the decomposition of a source function for a first-order inhomogeneous PDE is covered in Section 3. Some examples are given to illustrate the construction of a new homotopy in Section 4. Finally, some concluding remarks are given in Section 5.

### 2 New homotopy perturbation method

To illustrate the basic ideas of the homotopy analysis method, we consider the following nonlinear differential equation:

A ( u ) f ( r ) = 0 , r Ω (1)

with the boundary conditions

B ( u , u n ) = 0 , r Γ , (2)

where A is a general differential operator, B is a boundary operator, f ( r ) is a known analytical function, and Γ is the boundary of domain Ω.

The operator A in (1) can be rewritten as a sum of L and N, where L and N are linear and nonlinear parts of A, respectively, as follows:

L ( u ) + N ( u ) f ( r ) = 0 .

By the homotopy technique, we construct the homotopy

H ( v , p ) = ( 1 p ) ( L ( v ) L ( u 0 ) ) + p ( A ( v ) f ( r ) ) = 0 , (3)

which is equivalent to

H ( v , p ) = L ( v ) L ( u 0 ) + p L ( u 0 ) + p ( N ( v ) f ( r ) ) = 0 , (4)

where

p is an embedding parameter, u 0 is an initial approximation of (1), which satisfies the boundary conditions. As p changes from zero to unity, v ( r , p ) changes from u 0 to u ( r ) . In this technique, the convergence of a solution depends on the choice of u 0 , that is, we can have different approximate solutions for different u 0 .

Let us decompose the source function as f ( r ) = f 1 ( r ) + f 2 ( r ) . If we take

L ( u 0 ) = f 1 ( r )

in (3), we obtain the following homotopy:

H ( v , p ) = ( 1 p ) ( L ( v ) f 1 ( r ) ) + p ( A ( v ) f ( r ) ) = 0 , (5)

which is equivalent to

H ( v , p ) = L ( v ) f 1 ( r ) + p ( N ( v ) f 2 ( r ) ) = 0 . (6)

Obviously, from (6) we have

As the embedding parameter p changes from zero to unity, v ( r , p ) changes from L 1 ( f 1 ( r ) ) to u ( r ) .

According to He’s homotopy perturbation method, we can first use the embedding parameter p as a small parameter and assume that the solution of (6) can be written as a power series in p as follows:

v = v 0 + p v 1 + p 2 v 2 + p 3 v 3 + . (7)

Setting p = 1 , we get the approximate solution of (1)

u = lim p 1 v = v 0 + v 1 + v 2 + v 3 + .

In the next section, it is shown that the decomposition of the right-hand side function has a great effect on the amount of calculation and the speed of convergence of the solution.

### 3 Construction of a new homotopy based on the decomposition of a source function

Let us consider the following boundary value problem with the following inhomogeneous PDE:

(8)

(9)

where a, b, g and f are continuous functions in some region of the plane and g ( 0 ) = 0 . By solving this boundary value problem by the homotopy perturbation method, we obtain an approximate or exact solution u ( x , y ) . Before proceeding further, let us introduce the integral operator S defined in the following form:

S ( ) = 0 x ( ) d x . (10)

Then the derivative of operator S with respect to y is defined as

S y ( ) = y 0 x ( ) d x .

Rewriting f ( x , y ) as a sum of two functions f ( x , y ) = f 1 ( x , y ) + f 2 ( x , y ) and then constructing a homotopy in equation (6), we have

H ( v , p ) = ( 1 p ) ( v x f 1 ) + p ( v x + a ( x , y ) v y + b ( x , y ) g ( v ) f ) , (11)

which is equivalent to

H ( v , p ) = v x ( x , y ) f 1 ( x , y ) + p ( a ( x , y ) v y ( x , y ) + b ( x , y ) g ( v ( x , y ) ) f 2 ( x , y ) ) = 0 . (12)

Substituting (7) into (12), and equating the coefficient of the terms by the same power in p, we have

(13)

If there exists the relationship

a ( x , y ) [ S y ( f 1 ( x , y ) ) + h y ( y ) ] + b ( x , y ) g ( S ( f 1 ( x , y ) ) + h ( y ) ) = f 2 ( x , y ) (14)

between f 1 ( x , y ) and f 2 ( x , y ) , then we have from (13)

( v 1 ) x = 0 , v 1 ( 0 , y ) = 0 v 1 = 0

and

Hence the approximate or exact solution of problem (8)-(9) is obtained as

u ( x , y ) = v 0 ( x , y ) ,

and the source function is obtained from (14) in the following form:

f ( x , y ) = f 1 + a ( x , y ) [ S y ( f 1 ) + h y ( y ) ] + b ( x , y ) g ( S ( f 1 ) + h ( y ) ) . (15)

Therefore equation (15) is a necessary condition to accelerate the convergence. If equation (15) has a solution, then we obtain the approximate or exact solution of problem (8)-(9) in two steps.

However, it is not always possible to decompose the source function f ( x , y ) in such a way that the functions f 1 ( x , y ) and f 2 ( x , y ) have the relationship in (14). If we have such a case, then we are looking for an arbitrary P ( x , y ) such that the functions f 1 ( x , y ) and f 2 ( x , y ) have the following relationship:

a ( x , y ) [ S y ( f 1 ) + h y ( y ) ] + b ( x , y ) g ( S ( f 1 ) + h ( y ) ) = f 2 + P ( x , y ) . (16)

In this case, we get the approximate or exact solution of the problem in more than two steps. Moreover, we can get the solution in the form of series.

### 4 Analysis of the new homotopy perturbation method for linear problems

In this section we illustrate the theory of the new homotopy perturbation method given in Sections 2 and 3 for some special f ( x , y ) in the following linear problem:

(17)

(18)

where a and b are constants.

#### 4.1 Case 1: the source function f ( x , y ) is a polynomial

We consider that h ( y ) = 0 and f ( x , y ) is an nth-order polynomial in problem (17)-(18). Hence f ( x , y ) can be written in the following form:

f ( x , y ) = k = 0 n α + β = k A α β x α y β , (19)

where A α β are constants and α, β are natural numbers.

If the polynomial f ( x , y ) is decomposed as

f ( x , y ) = f 1 ( x , y ) + a S y ( f 1 ( x , y ) ) + b S ( f 1 ( x , y ) ) (20)

from (15), we can get the solution of the problem in two steps. We take f 1 ( x , y ) as an nth-order polynomial given as

f 1 ( x , y ) = k = 0 n α + β = k c α β x α y β , (21)

where c α β are constants. Then S ( f 1 ( x , y ) ) becomes

S ( f 1 ( x , y ) ) = k = 0 n α + β = k c α β α + 1 x α + 1 y β , (22)

and S y ( f 1 ( x , y ) ) becomes

S y ( f 1 ( x , y ) ) = k = 0 n α + β = k β > 1 β c α β α + 1 x α + 1 y β 1 . (23)

In order to determine the unknown coefficients c α β in terms of A α β , we substitute (21), (22), and (23) into (20)

f ( x , y ) = k = 0 n α + β = k c α β x α y β + k = 0 n α + β = k β > 1 β a c α β α + 1 x α + 1 y β 1 + k = 0 n α + β = k b c α β α + 1 x α + 1 y β , (24)

and then using (19), we get the following equations:

A 00 = c 00 , A 10 = c 10 + b c 00 , A 01 = c 01 , A 20 = c 20 + b 2 c 10 , A 02 = c 02 , A 11 = c 11 + 2 a c 02 + b c 01 ,

Then we obtain

c 00 = A 00 , c 10 = A 10 b A 00 , c 01 = A 01 , c 20 = A 20 b 2 ( A 10 b A 00 ) , c 02 = A 02 , c 11 = A 11 2 a A 02 b A 01 ,

As a result, constructing the following homotopy:

H ( v , p ) = ( 1 p ) ( v x f 1 ( x , y ) ) + p ( v x + a v y + b v f ( x , y ) ) = 0

leads us to reaching the solution in two steps.

#### 4.2 Case 2: the source function f ( x , y ) is a sum of two functions r ( x ) and t ( y )

Let us take f ( x , y ) = r ( x ) + t ( y ) in problem (17)-(18). Since f ( x , y ) is a continuous function, r ( x ) and t ( y ) are also continuous functions. Based on (15), f ( x , y ) = r ( x ) + t ( y ) is decomposed as follows:

f ( x , y ) = r 1 ( x ) + t 1 ( y ) + a S y ( r 1 ( x ) + t 1 ( y ) ) + a h y ( y ) + b S ( r 1 ( x ) + t 1 ( y ) ) + b h ( y ) . (25)

Then we can get the solution of the problem in two steps. Hence the function f ( x , y ) can be rewritten in the following form:

f ( x , y ) = r 1 ( x ) + t 1 ( y ) + a t 1 y ( y ) x + a h y ( y ) + b 0 x r 1 ( x ) d x + b t 1 ( y ) x + b h ( y ) . (26)

Since f ( x , y ) = r ( x ) + t ( y ) , we have

a t 1 y ( y ) x + b t 1 ( y ) x = 0 for  x 0 . (27)

The solutions of the above ordinary differential equation are

t 1 ( y ) = 0 (28)

and

t 1 ( y ) = e b a y . (29)

On the other hand, we have

t ( y ) = a h y ( y ) + b h ( y ) + t 1 ( y ) (30)

and

r ( x ) = b 0 x r 1 ( x ) d x + r 1 ( x ) . (31)

If either (28) or (29) satisfy equation (30), and equation (31) has a solution, then it is possible to decompose the function f ( x , y ) to accelerate the convergence of the solution. As a result, by constructing the following homotopy:

H ( v , p ) = ( 1 p ) ( v x r 1 ( x ) t 1 ( y ) ) + p ( v x + a v y + b v r ( x ) t ( y ) ) = 0 ,

we reach the solution in two steps.

### 5 Numerical examples

#### 5.1 Example 1

Consider the inhomogeneous linear boundary value problem (BVP) with constant coefficients

u x u y + u = e y + e x , u ( 0 , y ) = e y . (32)

By constructing the following homotopy:

H ( v , p ) = ( 1 p ) ( v x r 1 ( x ) t 1 ( y ) ) + p ( v x v y + v e x e y ) = 0 , (33)

which is equivalent to

H ( v , p ) = v x r 1 ( x ) t 1 ( y ) + p ( v y + v + r 1 ( x ) + t 1 ( y ) e x e y ) = 0 , (34)

we get the functions r 1 ( x ) , t 1 ( y ) as it is explained in Section 4.2. It is obvious that t 1 ( y ) = e y satisfies equation (30). Moreover, from equation (31), r 1 ( x ) can be found as follows:

r 1 ( x ) = e x + e x 2 . (35)

If we put the function r 1 ( x ) and t 1 ( y ) in the homotopy (33) or (34), we get

H ( v , p ) = v x e x + e x 2 e y + p ( v y + v + e x e x 2 ) = 0 , (36)

and we get the following:

Hence, the solution u ( x , y ) becomes

u ( x , y ) = v 0 = e x e x 2 + e y x + e y , (37)

which is the exact solution of the problem with minimum amount of calculation.

#### 5.2 Example 2

Consider the following inhomogeneous linear BVP with variable coefficients:

u x + y u y u = 2 y 2 + 2 x y 2 , u ( 0 , y ) = 0 . (38)

In this problem, the right-hand side function f ( x , y ) is a third-order polynomial. From (15), we have

f ( x , y ) = f 1 ( x , y ) + y S y ( f 1 ( x , y ) ) S ( f 1 ( x , y ) ) . (39)

Based on the source function f ( x , y ) = 2 y 2 + 2 x y 2 , if we take f 1 ( x , y ) as

f 1 ( x , y ) = a 1 x 3 + a 2 x 2 y + a 3 x y 2 + a 4 y 3 + a 5 x 2 + a 6 x y + a 7 y 2 + a 8 x + a 9 y + a 0 (40)

and substitute (40) into (39), we obtain

a 7 = 2 , a 1 = a 2 = a 3 = a 4 = a 5 = a 6 = a 8 = a 9 = a 0 = 0 .

Based on these results, we have f 1 ( x , y ) = 2 y 2 , f 2 ( x , y ) = 2 x y 2 and the homotopy is constructed in the following form:

H ( v , p ) = v x 2 y 2 + p ( y v y v 2 x y 2 ) = 0 , (41)

and we get the following:

Hence, we obtain the exact solution in the following form:

u ( x , y ) = v 0 = 2 x y 2 (42)

with short-length calculation.

#### 5.3 Example 3

Let us find the solution of following BVP:

u x + u y x u = y e x y + x , u ( 0 , y ) = 0 . (43)

As in the previous examples, using (15) we get

y e x y + x = f 1 + S y ( f 1 ) x S ( f 1 ) . (44)

Based on the source function f ( x , y ) = y e x y + x , we conclude that f 1 ( x , y ) is in the following form:

f 1 ( x , y ) = a 1 e x y + a 2 x e x y + a 3 x . (45)

Substituting (45) into (44), we obtain

a 1 = a 2 = 0 , a 3 = 1 .

Now, these results let us take f 1 ( x , y ) = y e x y , f 2 ( x , y ) = x and construct the homotopy in the following form:

H ( v , p ) = v x y e x y + p ( v y x v x ) = 0 , (46)

which leads to the following equations:

Hence, the exact solution becomes

u ( x , y ) = v 0 = e x y 1 , (47)

which is obtained by minimum amount of calculations.

#### 5.4 Example 4

Let us find the solution of the following first-order non-linear BVP:

u x + u y u 2 = e 2 x + 3 e x 1 , u ( 0 , y ) = 0 . (48)

Using (15) we get

e 2 x + 3 e x 1 = f 1 + S y ( f 1 ) [ S ( f 1 ) ] 2 . (49)

Based on the source function f ( x , y ) = e 2 x + 3 e x 1 , we find f 1 ( x , y ) is in the following form:

f 1 ( x , y ) = a 1 e x + a 2 e 2 x + a 3 . (50)

Substituting (50) into (49), we obtain

a 1 = 1 , a 2 = a 3 = 0 .

Now we can take f 1 ( x , y ) = e x , f 2 ( x , y ) = e 2 x + 2 e x 1 and construct the homotopy in the following form:

H ( v , p ) = v x e x + p ( v y v 2 + e 2 x 2 e x + 1 ) = 0 , (51)

which leads to the following equations:

Hence, we obtain the exact solution

u ( x , y ) = v 0 = e x 1 (52)

with short-length calculation.

### 6 Conclusion

In this paper, we employ a new homotopy perturbation method to obtain the solution of a first-order inhomogeneous PDE. In this method, each decomposition of the source function f ( x , y ) leads to a new homotopy. However, we develop a method to obtain the proper decomposition of f ( x , y ) which lets us obtain the solution with minimum computation and accelerate the convergence of the solution. This study shows that the decomposition of the source function has a great effect on the amount of computations and the acceleration of the convergence of the solution. Comparing to the standard one, decomposing the source function f ( x , y ) properly is a simple and very effective tool for calculating the exact or approximate solutions with less computational work.

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

The authors declare that the study was realized in collaboration with the same responsibility. All authors read, checked and approved the final manuscript.

### Acknowledgements

Dedicated to Professor Hari M Srivastava.

The research was supported by parts by the Scientific and Technical Research Council of Turkey (TUBITAK).

### References

1. Biazar, J, Azimi, F: He’s homotopy perturbation method for solving Helmoltz equation. Int. J. Contemp. Math. Sci.. 3, 739–744 (2008)

2. Biazar, J, Ghazvini, H: Convergence of the homotopy perturbation method for partial differential equations. Nonlinear Anal., Real World Appl.. 10, 2633–2640 (2009). Publisher Full Text

3. Sweilam, NH, Khader, MM: Exact solutions of some coupled nonlinear partial differential equations using the homotopy perturbation method. Comput. Math. Appl.. 58, 2134–2141 (2009). Publisher Full Text

4. Biazar, J, Ghazvini, H: Homotopy perturbation method for solving hyperbolic partial differential equations. Comput. Math. Appl.. 56, 453–458 (2008). Publisher Full Text

5. Junfeng, L: An analytical approach to the sine-Gordon equation using the modified homotopy perturbation method. Comput. Math. Appl.. 58, 2313–2319 (2009). Publisher Full Text

6. Neamaty, A, Darzi, R: Comparison between the variational iteration method and the homotopy perturbation method for the Sturm-Liouville differential equation. Bound. Value Probl.. 2010, Article ID 317369 (2010)

7. He, JH: Homotopy perturbation technique. Comput. Methods Appl. Mech. Eng.. 178, 257–262 (1999). Publisher Full Text

8. He, JH: A coupling method of homotopy technique and a perturbation technique for non-linear problems. Int. J. Non-Linear Mech.. 35, 37–43 (2000). PubMed Abstract | Publisher Full Text

9. He, JH: Homotopy perturbation method: a new nonlinear analytical technique. Appl. Math. Comput.. 135, 73–79 (2003). Publisher Full Text

10. He, JH: Comparison of homotopy method and homotopy analysis method. Appl. Math. Comput.. 156, 527–539 (2004). Publisher Full Text