Mathematical Problems in Engineering

Volume 2013, Article ID 503137, 13 pages

http://dx.doi.org/10.1155/2013/503137

## Optimal Homotopy Asymptotic Method to Nonlinear Damped Generalized Regularized Long-Wave Equation

^{1}Department of Mathematics, Abdul Wali Khan University, Mardan, Pakistan^{2}Faculty of Numerical and Physical Sciences, Islamia College Peshawar (Charted University Peshawar), Khyber Pakhtoonkhwa, Pakistan

Received 4 April 2013; Revised 23 July 2013; Accepted 31 July 2013

Academic Editor: Ashraf M. Zenkour

Copyright © 2013 R. Nawaz 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

A new semianalytical technique optimal homology asymptotic method (OHAM) is introduced for deriving approximate solution of the homogeneous and nonhomogeneous nonlinear Damped Generalized Regularized Long-Wave (DGRLW) equation. We tested numerical examples designed to confine the features of the proposed scheme. We drew 3D and 2D images of the DGRLW equations and the results are compared with that of variational iteration method (VIM). Results reveal that OHAM is operative and very easy to use.

#### 1. Introduction

Partial differential equations used in modeling different problems in physics, biology, chemical reactions, and engineering sciences problems are frequently too difficult to be solved exactly, and even if an exact solution is possible, the required calculations may be too difficult.

The DGRLW equation is a partial differential equation that describes the amplitude of long-wave, which takes the following form: Here , is an integer, is known function, and is the amplitude of the long-wave at the position and time . For , (1) features a balance between nonlinear and dispersive effects but also takes into account mechanisms of dissipation. In the physical sense, (1) with the dissipative term is suggested if the good predictive power is preferred; such type of problem arises in the bore propagation as well as for water waves [1].

The Equal Width (EW) Wave, Regularized Long-Wave (RLW), and Generalized Regularized Long-Wave (GRLW) equations are special cases of the DGRLW equation [2–4].

The EW equation corresponds to , , whereas the RLW or Benjamin-Bona-Mahony equation corresponds to , , and . On the other hand, the GRLW equation corresponds to , , and .

Different methods have been used for numerical solutions of DGRLW equations. VIM was presented by Demir et al. for numerical solutions for the DGRLW equation [5]. Yousefi et al. used Bernstein Ritz-Galerkin Method for solving the DGRLW equation [6]. Achouri et al. worked on an article called *“A fully Galerkin method for the damped generalized regularized long-wave (DGRLW) equation,” *namely, in [7]. For the mathematical theory and physical significance of DGRLW equation see [8–16] and references therein.

Marinca and Herişanu proposed semianalytical technique OHAM for deriving approximate solution of nonlinear problems of thin film flow of a fourth grade fluid down a vertical cylinder [17–19]. The method has been used by many researchers for obtaining numerical approximations of linear and nonlinear differential equations [20–24]. The author has successfully applied the OHAM for deriving approximate solution of Equal Width Wave equation, Burger equations, and tenth order boundary value problems [25, 26]. The convergence criterion of proposed method is similar to that of homotopy analysis method (HAM) and homotopy perturbation method (HPM), but this method is more efficient and flexible. To improve the efficiency and accuracy of OHAM, Herişanu and Marinca introduced more generalized and new advances in OHAM which shows that the auxiliary function includes the functions of physical parameters in addition to the convergence control parameter [27, 28].

Here, we investigate the approximate solution of the DGRLW equation with a variable coefficient using OHAM. The whole paper is divided into 3 sections. Section 2 is devoted to the analysis of the proposed method. In Section 3, solution of homogeneous and non-homogenous (DGRLW) equations is presented by OHAM, and absolute errors are also compared with VIM. The 3D and 2D images of the approximate solution and exact solution are also drawn. In all cases, the proposed method yields very encouraging results.

#### 2. Fundamental Theory of OHAM

Consider the partial differential equation of the following form: where is a linear operator and is nonlinear operator. is boundary operator, is an unknown function, and denote spatial and time variables, respectively, is the problem domain, and is a known function.

Using the basic idea of OHAM, the optimal homotopy is constructed which satisfies the following condition: where is an embedding parameter and is a nonzero auxiliary function for , . In such a case, (3) is called optimal homotopy equation. Clearly, we have when and , then and hold. Thus, as varies from , the solution approaches from to , where is obtained from (3) for : Next, we choose auxiliary function of the following general form:

Here are constants to be determined later.

To get an approximate solution, we expand in Taylor’s series about in the following manner,

Substituting (8) into (3) and equating the coefficient of like powers of , we obtain Zeroth order problem, given by (6); the first and second order problems are given by (9) and (10), respectively, and the general governing equations for are given by (11) as follows: where is the coefficient of in the expansion of about the embedding parameter :

Here for are a set of linear equations with the linear boundary conditions, which can be easily solved.

The convergence of the series in (8) depends upon the auxiliary constants . If it is convergent at , then

Substituting (13) into (1), it results with the following expression for residual:

If , then will be the exact solution.

For computing the optimal values of auxiliary constants, , there are many methods available like Galerkin’s, Ritz, Least Squares, and Collocation method. One can apply the method of Least Squares as under the following: where is the residual,

The constants can also be determined by another method as under: at any time , where . The convergence depends upon constants , which can be optimally identified and minimized by (18).

*Example 1. *Consider (1) with , , and which in the simplest form is given as
The initial condition is and exact solution given by

*Zeroth Order Problem*. Consider the following:
Its solution is given as under

*First Order Problem*. Consider the following:
Its solution is as follows:

*Second Order Problem*. Consider the following:
Its solution is under

*Third Order Problem*. Consider the following:

Its solution is given as follows: The third order approximate solution is given by the following equation:

The constants , , and are calculated using the Least Squares, we have their optimal values as follows:

The 3rd order OHAM solution yields very encouraging results after comparing with 3rd order VIM solution [5]. Tables 1(a–c), and Figures 1, 2, 3, and 4 show the effectiveness of OHAM for , and .

*Example 2. *Consider (1) with , and which in the simplest form is given as

The initial condition is and exact solution given by

*Zeroth Order Problem*. Consider the following:
Its solution is given by

*First Order Problem*. Consider the following:
Its solution is

*Second Order Problem*

Its approximate solution is obtained in similar manner. The second order approximate solution is given by

Using method of Least Squares, the optimal values of constants are computed and are given as under

Table 2(a) shows the effectiveness of OHAM for , , and , while Table 2(b), and Figures 5, 6, 7, and 8 shows the effectiveness of OHAM for various values of and .

*Example 3. *Consider (1) with , , and − which in the simplest form is given as
The initial condition is and exact solution given by

*Zeroth Order Problem*. Consider the following:
Its solution is as follows:

*First Order Problem*. Consider the following:

Its solution is

The first order approximate solution is given by The constants is calculated using the Least Squares that we have its optimal values as follows: The first order optimum solution using OHAM is as follows:

The first order OHAM solution yields very encouraging results after comparing with 2nd order VIM solution [5]. Tables 3(a–d), and Figures 9, 10, 11, and 12 show the effectiveness of OHAM for , , and .

*Example 4. *Let us consider the inhomogeneous DGRLW equation:
where and .

The initial condition is
and exact solution given by

According to OHAM scheme presented in Section 2.

*Zeroth Order Problem*. Consider the following:
Its solution is

*First Order Problem*. Consider the following:
Its solution is as follows:

*Second Order Problem*. Consider the following:

Its approximate solution is under The second order approximate solution is given by the following equation, Using the method of Least Squares the optimum values of and are computed which are as follows:

The 2nd order OHAM solution yields very encouraging results after comparing with 2nd order VIM solution [5]. Tables 4(a–d), and Figures 13, 14, 15, and 16 show the effectiveness of OHAM for , , , and .

#### 3. Conclusion

In this paper, the OHAM has been successfully implemented for the approximate solution of solutions of the Nonlinear Damped Generalized Regularized Long-Wave (DGRLW) equations. The results obtained by OHAM are very consistent in comparison with VIM.

#### References

- T. Achouri and K. Omrani, “Numerical solutions for the damped generalized regularized long-wave equation with a variable coefficient by Adomian decomposition method,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 14, no. 5, pp. 2025–2033, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - K. Al-Khaled, S. Momani, and A. Alawneh, “Approximate wave solutions for generalized Benjamin-Bona-Mahony-Burgers equations,”
*Applied Mathematics and Computation*, vol. 171, no. 1, pp. 281–292, 2005. View at Publisher · View at Google Scholar · View at MathSciNet - D. Kaya, “A numerical simulation of solitary-wave solutions of the generalized regularized long-wave equation,”
*Applied Mathematics and Computation*, vol. 149, no. 3, pp. 833–841, 2004. View at Publisher · View at Google Scholar · View at MathSciNet - J. I. Ramos, “Solitary waves of the EW and RLW equations,”
*Chaos, Solitons and Fractals*, vol. 34, no. 5, pp. 1498–1518, 2007. View at Publisher · View at Google Scholar · View at MathSciNet - D. D. Demir, N. Bildik, A. Konuralp, and A. Demir, “The numerical solutions for the damped generalized regularized long-wave equation by variational method,”
*World Applied SciencesJournal*, vol. 13, pp. 1308–1317, 2011. View at Google Scholar - S. A. Yousefi, Z. Barikbin, and M. Behroozifar, “Bernstein Ritz-Galerkin method for solving the damped generalized regularized long-wave (DGRLW) equation,”
*International Journal of Nonlinear Science*, vol. 9, no. 2, pp. 151–158, 2010. View at Google Scholar · View at MathSciNet - T. Achouri, M. Ayadi, and K. Omrani, “A fully Galerkin method for the damped generalized regularized long-wave (DGRLW) equation,”
*Numerical Methods for Partial Differential Equations*, vol. 25, no. 3, pp. 668–684, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - T. B. Benjamin, J. L. Bona, and J. J. Mahony, “Model equations for long waves in nonlinear dispersive systems,”
*Philosophical Transactions of the Royal Society of London A*, vol. 272, no. 1220, pp. 47–78, 1972. View at Google Scholar · View at MathSciNet - D. H. Peregrine, “Calculations of the development of an undular bore,”
*Journal of Fluid Mechanics*, vol. 25, no. 2, pp. 321–326, 1966. View at Publisher · View at Google Scholar - M. Mei, “Large-time behavior of solution for generalized Benjamin-Bona-Mahony-Burgers equations,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 33, no. 7, pp. 699–714, 1998. View at Publisher · View at Google Scholar · View at MathSciNet - A. Durán and M. A. López-Marcos, “Numerical behaviour of stable and unstable solitary waves,”
*Applied Numerical Mathematics*, vol. 42, no. 1–3, pp. 95–116, 2002. View at Publisher · View at Google Scholar · View at MathSciNet - H. Zhang, G. M. Wei, and Y. T. Gao, “On the general form of the Benjamin-Bona-Mahony equation in fluid mechanics,”
*Czechoslovak Journal of Physics*, vol. 52, no. 3, pp. 373–377, 2002. View at Publisher · View at Google Scholar · View at MathSciNet - J. L. Bona and R. Smith, “The initial-value problem for the Korteweg-de Vries equation,”
*Philosophical Transactions of the Royal Society of London A*, vol. 278, no. 1287, pp. 555–601, 1975. View at Google Scholar · View at MathSciNet - L. A. Medeiros and M. M. Miranda, “Weak solutions for a nonlinear dispersive equation,”
*Journal of Mathematical Analysis and Applications*, vol. 59, no. 3, pp. 432–441, 1977. View at Google Scholar · View at MathSciNet - L. A. Medeiros and G. P. Menzala, “Existence and uniqueness for periodic solutions of the Benjamin-Bona-Mahony equation,”
*SIAM Journal on Mathematical Analysis*, vol. 8, no. 5, pp. 792–799, 1977. View at Google Scholar · View at MathSciNet - C. Yong, L. Biao, and Z. Hongqing, “Exact solutions for two nonlinear wave equations with nonlinear terms of any order,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 10, no. 2, pp. 133–138, 2005. View at Publisher · View at Google Scholar · View at MathSciNet - V. Marinca and N. Herişanu, “Application of optimal homotopy asymptotic method for solving nonlinear equations arising in heat transfer,”
*International Communications in Heat and Mass Transfer*, vol. 35, no. 6, pp. 710–715, 2008. View at Publisher · View at Google Scholar · View at Scopus - V. Marinca, N. Herişanu, C. Bota, and B. Marinca, “An optimal homotopy asymptotic method applied to the steady flow of a fourth-grade fluid past a porous plate,”
*Applied Mathematics Letters*, vol. 22, no. 2, pp. 245–251, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - V. Marinca, N. Herişanu, and I. Nemeş, “Optimal homotopy asymptotic method with application to thin film flow,”
*Central European Journal of Physics*, vol. 6, no. 3, pp. 648–653, 2008. View at Publisher · View at Google Scholar · View at Scopus - S. Iqbal, M. Idrees, A. M. Siddiqui, and A. R. Ansari, “Some solutions of the linear and nonlinear Klein-Gordon equations using the optimal homotopy asymptotic method,”
*Applied Mathematics and Computation*, vol. 216, no. 10, pp. 2898–2909, 2010. View at Publisher · View at Google Scholar · View at MathSciNet - S. Haq, M. Idrees, and S. Isalam, “Application of optimal homotopy asymptotic method to eighth order boundary value problems,”
*Journal of Applied Mathematics and Computing*, vol. 2, no. 4, pp. 38–47, 2008. View at Google Scholar - M. Idrees, S. Haq, and S. Islam, “Application of optimal homotopy asymptotic method to fourth order boundary value problems,”
*World Applied Sciences Journal*, vol. 9, no. 2, pp. 131–137, 2010. View at Google Scholar - M. Idrees, S. Islam, S. Haq, and S. Islam, “Application of the optimal homotopy asymptotic method to squeezing flow,”
*Computers & Mathematics with Applications*, vol. 59, no. 12, pp. 3858–3866, 2010. View at Publisher · View at Google Scholar · View at MathSciNet - M. Idrees, S. Haq, and S. Islam, “Application of optimal homotopy asymptotic method to special sixth order boundary value problems,”
*World Applied Sciences Journal*, vol. 9, no. 2, pp. 138–143, 2010. View at Google Scholar - S. Islam, R. Nawaz, M. Arif, and I. A. Shah, “Application of optimal homotopy asymptotic method to the equal width wave and Burger equations,”
*Life Science Journal*, vol. 9, no. 4, pp. 2380–2386, 2012. View at Google Scholar - R, M. N. Nawaz Khalid, S. Islam, and S. Yasin, “Solution of tenth order boundary value problems using optimal homotopy asymptotic method,”
*Canadian Journal Computing Mathematics Natural Sciences*, vol. 1, no. 2, 2010. View at Google Scholar - N. Herişanu and V. Marinca, “Explicit analytical approximation to large-amplitude non-linear oscillations of a uniform cantilever beam carrying an intermediate lumped mass and rotary inertia,”
*Meccanica*, vol. 45, no. 6, pp. 847–855, 2010. View at Publisher · View at Google Scholar · View at MathSciNet - N. Herişanu and V. Marinca, “Accurate analytical solutions to oscillators with discontinuities and fractional-power restoring force by means of the optimal homotopy asymptotic method,”
*Computers & Mathematics with Applications*, vol. 60, no. 6, pp. 1607–1615, 2010. View at Publisher · View at Google Scholar · View at MathSciNet