Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 362816, 8 pages

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

Research Article

Application of Optimal Homotopy Asymptotic Method to Doubly Wave Solutions of the Coupled Drinfel’d-Sokolov-Wilson Equations

H. Ullah, S. Islam, M. Idrees, and R. Nawaz

Department of Mathematics, Abdul Wali Khan University, Mardan, Pakistan

Received 7 August 2013; Accepted 5 November 2013

Academic Editor: Hai-Feng Huo

Copyright © 2013 H. Ullah 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.

The approximate solution of the doubly periodic wave solutions of the coupled Drinfel’d-Sokolov-Wilson equations has been considered by using the optimal homotopy asymptotic method (OHAM). We obtained the numerical solution of the problem and compared that with the OHAM solution. The obtained solutions show that OHAM is effective, simpler, easier, and explicit and gives a suitable way to control the convergence of the approximate solution.

1. Introduction

The coupled nonlinear partial differential equations (NPDEs) are widely used in applied mathematics, physics, and engineering sciences to offer the description of complex phenomena. Here we consider doubly periodic wave solutions of the coupled Drinfel’d-Sokolov-Wilson equation of the form [1]

with

The exact and explicit solution of the NPDEs in mathematical physics, engineering, and science plays an important role. The exact solution of NPDEs cannot be found easily as all NPDEs have infinitely many solutions. The analytical and exact solution of such problems is either not available in the literature or may be found by using transformation based on the invariance group analysis method [2], the Lie infinitesimal criterion [3], the symbolic computation [4], and the Backlund transformation [5]. All these methods reduced the complex equations into simple equations by using the transformation. In the literature most of the methods like the variational iterative method (VIM) [6], Adomian decomposition method (ADM) [7], differential transform method (DTM) [8], and homotopy perturbation method (HPM) [9] have been used for the solution of weakly NPDEs and few for strongly NPDEs. To tackle the strongly NPDEs the perturbation methods were introduced [10, 11]. These methods contain a small parameter which cannot be found easily. New analytic methods such as the artificial parameters method [12], homotopy analysis method (HAM) [13], and homotopy perturbation method (HPM) [9] were introduced. These methods combined the homotopy with the perturbation techniques. Recently, Vasile Marinca et al. introduced OHAM [14–18] for the solution of nonlinear problems which made the perturbation methods independent of the assumption of small parameters and huge computational work.

The motivation of this paper is to boost OHAM for the solution of coupled NPDEs. In [19–24] OHAM has been proved to be valuable for obtaining an approximate solution of the single partial differential equation (PDE). Before these coupled NPDEs were not solved by OHAM. We have proved that OHAM is useful and reliable for NPDEs, showing its validity and great potential for the solution of transient physical phenomena in science and engineering.

In the succeeding section, the basic idea of OHAM is formulated for the solution of NPDEs. The effectiveness and efficiency of OHAM are shown in Section 3.

2. Fundamental Mathematical Theory of OHAM

Let us see the partial differential equation of the following form:

where is a differential operator, is an unknown function, and denote spatial and temporal independent variables, respectively, is the boundary of , and is a known analytic function. can be divided into two parts and such that

is the simpler part of the partial differential equation which is easier to solve, and contains the remaining part of .

According to OHAM, one can construct an optimal homotopy which satisfies

where the auxiliary function is nonzero for and . Equation (5) is called an optimal homotopy equation. Clearly, we have

Obviously, when and we obtain

respectively. Thus, as varies from to , the solution approaches from to , where is obtained from (5) for :

Next, we choose the auxiliary function in the form

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

Substituting (10) into (5) and equating the coefficient of the like powers of , we obtain the zeroth order problem, given by (8), the first and second order problems are given by (11)-(12), respectively, and the general governing equations for are given by (13):

where are the coefficients of in the expansion of about the embedding parameter :

It should be underscored that the for are governed by the linear equations with linear boundary conditions that come from the original problem, which can be easily solved.

It has been observed that the convergence of the series equation (10) depends upon the auxiliary constants . If it is convergent at , one has

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

In actual computation .

If , then is the exact solution of the problem. Generally it does not happen, especially in nonlinear problems.

For the determinations of auxiliary constants, , , there are different methods like Galerkin’s method, the Ritz method, the least squares method, and, the collocation method. One can apply the method of least squares as under

The th order approximation can be obtained by these constants. The constants can also be determined by another method as under

at any time , where .

The more general auxiliary function is useful for convergence, which depends upon constants , can be optimally identified by (18), and is useful in error minimization.

3. Application of OHAM to Doubly Periodic Wave Solutions of the Coupled Drinfel’d-Sokolov-Wilson Equation

To demonstrate the effectiveness of the extended formulation of OHAM for coupled nonlinear partial differential equations (NPDEs), we consider the doubly periodic wave solutions of the coupled Drinfel’d-Sokolov-Wilson equations (1) with the boundary condition (2).

Applying the method formulated in Section 2 leads to the following:

We consider

3.1. Zeroth Order System

We have

with initial conditions

Its solution

3.2. First Order System

We have

Its solutions

3.3. Second Order System

We have

with

Its solutions

Adding (26), (29), and (28), we obtain

For the calculations of the constants , , , and using (30) in (1) and applying the procedure mentioned in (16)–(19), we get

4. Results and Discussions

The formulation presented in Section 2 provides highly accurate solutions for the problems demonstrated in Section 3. We have used Mathematica 7 for most of our computational work. In Tables 1 and 3, we have presented absolute errors for and at a spatial domain [0, 0.4] for , , , , , and , while in Tables 2 and 4 the convergence of the OHAM solution is given, through first and second order absolute errors at time and . Here we observe that the OHAM solution converges rapidly with increasing order of approximation. From Tables 1–4 and Figures 1, 2, 3, 4, 5, 6, 7, 8, 9, and 10 it is evident that the OHAM results are nearly identical to the numerical results. Here the results are very consistent with the increasing time.

Table 1: Absolute error of OHAM solution of corresponding to the numerical solution.

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

Table 2: Comparison of first order and second order errors of corresponding to the numerical solution at time , and .

First order absolute error
Second order absolute error

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

Table 3: Absolute error of OHAM solution of corresponding to the numerical solution.

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

Table 4: Comparison of first order and second order errors of corresponding to the numerical solution at time , and .

First order absolute error
Second order absolute error

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

Figure 1: 3D approximate solution of for .

Figure 2: 3D numerical solution of for .

Figure 3: 2D approximate solution of for .

Figure 4: 2D numerical solution for .

Figure 5: 2D numerical, zeroth, first, and second order solutions of for .

Figure 6: 3D approximate solution of for .

Figure 7: 3D exact solution of for .

Figure 8: 2D approximate solution of for .

Figure 9: 2D numerical solution of for .

Figure 10: 2D numerical, zeroth, first, and second order solution of for .

5. Conclusion

In this paper, we have seen the effectiveness of OHAM [16–20] in doubly periodic wave solutions of the coupled Drinfel’d-Sokolov-Wilson equation. By applying the basic idea of OHAM to doubly periodic wave solutions of the coupled Drinfel’d-Sokolov-Wilson equation, we found it simpler in applicability, more convenient to control convergence, and involving less computational overhead. Therefore, OHAM shows its validity and great potential for the solution of time dependant problems in science and engineering.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

M. Inc, “On numerical doubly periodic wave solutions of the coupled Drinfel'd-Sokolov-Wilson equation by the decomposition method,” Applied Mathematics and Computation, vol. 172, no. 1, pp. 421–430, 2006.

M. Torrisi, R. Tracinà, and A. Valenti, “A group analysis approach for a nonlinear differential system arising in diffusion phenomena,” Journal of Mathematical Physics, vol. 37, no. 9, pp. 4758–4767, 1996.

P. G. Drzain and R. S. Johnson, An Introduction Discussion the Theory of Solution and Its Diverse Applications, Cambridge University Press, Cambridge, Mass, USA, 1989.

B. Abdel-Hamid, “Exact solutions of some nonlinear evolution equations using symbolic computations,” Computers and Mathematics with Applications, vol. 40, no. 2, pp. 291–302, 2000.

G. Bluman and S. Kumei, “On the remarkable nonlinear diffusion equation (∂/∂x)[a(u + b)⁻²(∂u/∂x)]−(∂u/∂t) = 0,” Journal of Mathematical Physics, vol. 21, no. 5, pp. 1019–1023, 1980.

C. Chun, “Fourier-series-based variational iteration method for a reliable treatment of heat equations with variable coefficients,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 10, no. 11-12, pp. 1383–1388, 2009.

S. H. Chowdhury, “A comparison between the modified homotopy perturbation method and adomian decomposition method for solving nonlinear heat transfer equations,” Journal of Applied Sciences, vol. 11, no. 8, pp. 1416–1420, 2011.

H. Yaghoobi and M. Torabi, “The application of differential transformation method to nonlinear equations arising in heat transfer,” International Communications in Heat and Mass Transfer, vol. 38, no. 6, pp. 815–820, 2011.

D. D. Ganji, “The application of He's homotopy perturbation method to nonlinear equations arising in heat transfer,” Physics Letters A, vol. 355, no. 4-5, pp. 337–341, 2006.

R. Bellman, Perturbation Techniques in Mathematics, Physics, and Engineering, Holt, Rinehart and Winston, New York, NY, USA, 1964.

J. D. Cole, Perturbation Methods in Applied Mathematics, Blaisedell, Waltham, Mass, USA, 1965.

G. L. Liu, “New research direction in singular perturbation theory: artificial parameter approach and inverse perturbation technique,” in Proceedings of the 7th Conference on Modern Mathematics and Mechanics, 1997.

S. J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problems [Ph.D. thesis], Shanghai Jiao Tong University, Shanghai, China, 1992.

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.

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.

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.

N. Herisanu, V. Marinca, T. dordea, and Gh. Madescu, “A new analytic approach to nonlinear vibration of an electrical machine,” Proceeding of the Romanian Academy, vol. 9, pp. 229–236, 2008.

V. Marinca and N. Herişanu, “Determination of periodic solutions for the motion of a particle on a rotating parabola by means of the optimal homotopy asymptotic method,” Journal of Sound and Vibration, vol. 329, no. 9, pp. 1450–1459, 2010.

H. Ullah, S. Islam, M. Idrees, and M. Arif, “Solution of boundary layer problems with heat transfer by optimal homotopy asymptotic method,” Abstract and Applied Analysis, vol. 2013, Article ID 324869, 10 pages, 2013.

M. Idrees, S. Islam, S. I. A. Tirmizi, and S. Haq, “Application of the optimal homotopy asymptotic method for the solution of the Korteweg-de Vries equation,” Mathematical and Computer Modelling, vol. 55, no. 3-4, pp. 1324–1333, 2012.

R. Nawaz, S. Islam, I. A. Shah, M. Idrees, and H. Ullah, “Optimal homotopy asymptotic method to nonlinear damped generalized regularized long-wave equation,” Mathematical Problems in Engineering, vol. 2013, Article ID 503137, 13 pages, 2013.

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.

M. Idrees, S. Islam, S. Haq, and G. Zaman, “Application of the optimal homotopy asymptotic method to squeezing flow,” Computers and Mathematics with Applications, vol. 59, pp. 3858–3866, 2010.

R. Nawaz, H. Ullah, S. Islam, and M. Idrees, “Application of optimal homotopy asymptotic method to Burger equations,” Journal of Applied Mathematics, vol. 2013, Article ID 387478, 8 pages, 2013.