Journal of Applied Mathematics

Volume 2012 (2012), Article ID 387365, 9 pages

doi:10.1155/2012/387365

Research Article

Removal of Noise Oscillation Term Appearing in the Nonlinear Equation Solution

Yasir Khan,¹ Hector Vázquez-Leal,² and Luis Hernandez-Martínez³

¹Department of Mathematics, Zhejiang University, Hangzhou 310027, China

²School of Electronic Instrumentation, University of Veracruz, Xalapa, VER, Mexico

³Electronics Department, National Institute for Astrophysics, Optics and Electronics, 72000 Tonantzintla, PUE, Mexico

Received 27 April 2012; Accepted 4 July 2012

Academic Editor: J. Biazar

Copyright © 2012 Yasir Khan et al. 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.

This paper suggests a novel modified Laplace method for removal of noise oscillation term appearing in the nonlinear equation solutions. The modified method overcomes the noise oscillation during the iteration procedure by suitable choice of an initial solution. Several examples are tested, and the obtained results suggest that this newly developed technique could lead to a promising tool and powerful improvement for many applications in differential and integral equations.

1. Introduction

Nonlinear science emerged in its present form following a series of decisive analytic, numerical, and experimental development took place in close interaction in the last three decades. And the application of many nonlinear science partial differential equations has attracted a great deal of attention in recent years. For example, partial differential equations are increasing used to model many problems in biology, chemistry, economic, engineering, physics, and other areas of applications. Many analytical and numerical techniques have been developed by various scientists to cope with the nonlinearity of such problems. In reality, the Laplace transform is one of only few methods that can be useful to linear systems with periodic or discontinuous driving inputs. In spite of its great usefulness in solving linear problems, however, the Laplace transform is totally incapable of handling nonlinear equations because of the difficulties that are caused by the nonlinear terms. Various ways have been proposed recently to deal with these nonlinearities such as the Adomian decomposition method [1–3], the Laplace decomposition method [4, 5], the homotopy perturbation method [6–14], variational approach [15, 16], Hamiltonian approach [17, 18], Boubaker polynomials expansion scheme [19, 20], and recently the homotopy perturbation transform method was proposed by Khan and Wu [21].

In this paper, we proposed a new method to solve any kind of nonlinear equations. This method is called modified homotopy perturbation transform method (MHPTM). It is important to note that the modified technique works effectively independent of other phenomena in some cases, or it may be combined with the powerful phenomena of the self-canceling “noise terms’’. We aim to extend the works of [21] and make further progress beyond the achievements made so far in this regard. The proposed modification will accelerate the rapid convergence of the series solution if compared with the standard homotopy perturbation transform method (HPTM), and therefore provides a major progress. The modified technique has been shown to be computationally efficient in several examples that are important to researchers in applied fields. In addition, the modified technique may give the exact solution for nonlinear equations without any need of the so-called Adomian polynomials. Although the modified technique needs only a slight variation from the standard homotopy perturbation transform method [21], the results are promising in that it minimizes the size of calculations needed. While this slight variation is rather simple, it does demonstrate the reliability and the power of the proposed modification. The modified homotopy perturbation transform method is much easier to implement as compared with the Adomian decomposition method and the Laplace decomposition method, where huge complexities are involved. Several examples are given to reconfirm the efficiency and accuracy of the proposed modified homotopy perturbation transform method.

2. Description of the Method

To illustrate the basic idea of this method, we consider a general nonlinear nonhomogeneous partial differential equation with initial condition of the following form:

where is the second order linear differential operator , the linear differential operator of less order than , representing the general nonlinear differential operator, and as the source term.

According to the homotopy perturbation transform method [21], we apply the Laplace transform (denoted throughout this paper by ) on both sides of (2.1) as follows:

Using the differentiation property of Laplace transform, we have

Operating with Laplace inverse on both sides of (2.4) gives

where represents the term arising from the source term and prescribed initial condition. The homotopy perturbation transform method [21] admits a solution in the form

The nonlinear term is decomposed as

where are He’s polynomials [22, 23] of which can be calculated by the following formula:

Using (2.7) and (2.6) in (2.5), we get

Comparing the coefficient of like powers of , the following approximations are obtained:

Clearly, the general recursive relation is given by

The initial solution is important, and the choice of (2.10) as the initial solution always leads to noise oscillation/term during the iteration procedure. It is important to note that these terms may appear for inhomogeneous problems, whereas homogeneous problems do not generate noise terms. A necessary condition for the generation of the noise terms for inhomogeneous problems is that the zeroth component must contain the exact solution among other terms. The noise terms are defined as the identical terms, with opposite signs, that may appear in various components . It was formally shown that by canceling the noise terms that appear in and from , even though contains further terms, the remaining noncancelled terms of may give the exact solution of the inhomogeneous problem. In order to overcome this shortcoming, we decompose into two parts as

Instead of the iteration procedure, (2.10), (2.11), and (2.12), we suggest the following modification:

The solution through the modified homotopy perturbation transform method highly depends upon the choice of and . We will show how to suitably choose and through the use of examples.

3. Nonlinear Klein-Gordan Equations

In order to elucidate the solution procedure of the modified method and how to choose a suitable initial solution, we consider first the nonlinear Klein-Gordon equations.

Example 3.1. Consider the following nonlinear Klein-Gordon equation:

Taking Laplace transform of both sides of (3.1) gives

The initial conditions (3.2) imply

The inverse Laplace transform applied to (3.4) results in

Applying the homotopy perturbation transform method [21], we get

where are He’s polynomials [22, 23] that represent the nonlinear terms. In view of (3.6), recursive relation is given as follows:

It is important to recall here that the noise terms appear between the two adjacent components, for example, in and . The noise terms between the adjacent components can be cancelled using our modified version. The exact solution is

According to the modified homotopy perturbation transform method, we first set

As suggested before, we split into two parts as

According to the iteration algorithm equation (2.14), we obtain

One iteration leads to the exact solution . It is obvious from this example that “noise terms” appeared in the standard homotopy perturbation transform method can be completely eliminated in our solution procedure.

Example 3.2. Consider the following nonlinear Klein-Gordon equation:

By applying the aforesaid method subject to the initial conditions, we have

The inverse of Laplace transform implies that

In view of the homotopy perturbation transform method [21], we have

where the nonlinear operator is decomposed as in (2.8) in terms of the He’s polynomials. Matching both sides of (3.15), the components of can be defined as follows:

The solution procedure can continue, and we find the remaining components are all zero. We, therefore, obtain a closed-form solution which is .

4. Conclusion

This paper suggests an effective modification of the Laplace transform method. A suitable choice of and results in an extremely simple solution procedure. The modified version is also valid for other nonlinear equations, and this paper can be used as a standard paradigm for other applications.

Acknowledgment

The authors wish to express their cordial thanks to the anonymous referees for useful suggestions and comments.

References

G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, Kluwer Academic, Boston, Mass, USA, 1994.

E. Hetmaniok, D. Słota, R. Witua, and A. Zielonka, “Comparison of the Adomian decomposition method and the variational iteration method in solving the moving boundary problem,” Computers and Mathematics with Applications, vol. 61, no. 8, pp. 1931–1934, 2011.

J. Biazar, “Solution of the epidemic model by Adomian decomposition method,” Applied Mathematics and Computation, vol. 173, no. 2, pp. 1101–1106, 2006.

S. A. Khuri, “A laplace decomposition algorithm applied to a class of nonlinear differential equations,” Journal of Applied Mathematics, vol. 1, no. 4, pp. 141–155, 2001.

Y. Khan, “An effective modification of the laplace decomposition method for nonlinear equations,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 10, pp. 1373–1376, 2009.

J. H. He, “Homotopy perturbation method for solving boundary value problems,” Physics Letters A, vol. 350, no. 1-2, pp. 87–88, 2006.

L. Xu, “He's homotopy perturbation method for a boundary layer equation in unbounded domain,” Computers and Mathematics with Applications, vol. 54, no. 7-8, pp. 1067–1070, 2007.

H. Koçak, T. Öziş, and A. Yildirim, “Homotopy perturbation method for the nonlinear dispersive K(m, n, 1) equations with fractional time derivatives,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 20, no. 2, pp. 174–185, 2010.

A. R. Ghotbi, A. Barari, M. Omidvar, and G. Domairry, “Application of homotopy perturbation and variational iteration methods to SIR epidemic model,” Journal of Mechanics in Medicine and Biology, vol. 11, no. 1, pp. 149–161, 2011.

M. Turkyilmazoglu, “Convergence of the homotopy perturbation method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 12, pp. 9–14, 2011.

E. Hetmaniok, I. Nowak, D. Slota, and R. Witula, “Application of the homotopy perturbation method for the solution of inverse heat conduction problem,” International Communications in Heat and Mass Transfer, vol. 39, no. 1, pp. 30–35, 2011.

H. Koçak and A. Yıldırım, “An efficient algorithm for solving nonlinear age-structured population models by combining homotopy perturbation and Pade techniques,” International Journal of Computer Mathematics, vol. 88, no. 3, pp. 491–500, 2011.

S. S. Ganji, A. Barari, and D. D. Ganji, “Approximate analysis of two-massspring systems and buckling of a column,” Computers and Mathematics with Applications, vol. 61, no. 4, pp. 1088–1095, 2011.

J. Biazar and Behzad Ghanbari, “The homotopy perturbation method for solving neutral functional-differential equations with proportional delays,” Journal of King Saud University, vol. 24, no. 4, pp. 33–37, 2012.

Y. Khan, N. Faraz, and A. Yildirim, “New soliton solutions of the generalized Zakharov equations using He's variational approach,” Applied Mathematics Letters, vol. 24, no. 6, pp. 965–968, 2011.

L. Xu and N. Zhang, “A variational approach to analyzing catalytic reactions in short monoliths,” Computers and Mathematics with Applications, vol. 58, no. 11-12, pp. 2460–2463, 2009.

Y. Khan, Q. Wu, H. Askari, Z. Saadatnia, and M. Kalami-Yazdi, “Nonlinear vibration analysis of a rigid rod on a circular surface via hamiltonian approach,” Mathematical and Computational Applications, vol. 15, no. 5, pp. 974–977, 2010.

L. Xu, “A Hamiltonian approach for a plasma physics problem,” Computers and Mathematics with Applications, vol. 61, no. 8, pp. 1909–1911, 2011.

W. X. Yue, H. Koçak, D. H. Zhang, and A. Yildirim, “A second attempt to establish an analytical expression to steam-water dipole orientation parameter using the boubaker polynomials expansion scheme,” Journal of Structural Chemistry, vol. 52, no. 1, pp. 106–110, 2011.

H. Koçak, Z. Dahong, and A. Yildirim, “A range-free method to determine Antoine vapor-pressure heat transfer-related equation coefficients using the Boubaker polynomial expansion scheme,” Russian Journal of Physical Chemistry A, vol. 85, no. 5, pp. 900–902, 2011.

Y. Khan and Q. Wu, “Homotopy perturbation transform method for nonlinear equations using He's polynomials,” Computers and Mathematics with Applications, vol. 61, no. 8, pp. 1963–1967, 2011.

Y. Khan, M. Madani, A. Yıldırım, M. A. Abdou, and N. Faraz, “A new approach to Van der Pol's oscillator problem,” Zeitschrift für Naturforschung, vol. 66, pp. 620–624, 2011.

A. Ghorbani, “Beyond Adomian polynomials: He polynomials,” Chaos, Solitons and Fractals, vol. 39, no. 3, pp. 1486–1492, 2009.