• Views 734
• Citations 2
• ePub 24
• PDF 603
`International Journal of Differential EquationsVolume 2014, Article ID 106934, 11 pageshttp://dx.doi.org/10.1155/2014/106934`
Research Article

## An Extension of the Optimal Homotopy Asymptotic Method to Coupled Schrödinger-KdV Equation

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

Received 26 December 2013; Accepted 13 March 2014; Published 7 May 2014

Academic Editor: Patricia J. Y. Wong

Copyright © 2014 Hakeem 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

We consider the approximate solution of the coupled Schrödinger-KdV equation by using the extended optimal homotopy asymptotic method (OHAM). We obtained the extended OHAM solution of the problem and compared with the exact, variational iteration method (VIM) and homotopy perturbation method (HPM) solutions. The obtained solution shows that extended OHAM is effective, simpler, easier, and explicit and gives a suitable way to control the convergence of the approximate solution.

#### 1. Introduction

The nonlinear Schrödinger equations are of great interest due to their numerous applications in physical phenomena. The coupled Schrödinger-KdV equations are extensively used to model nonlinear dynamics of one-dimensional Langmuir and ion acoustic waves in the system of coordinates moving at the speed of ion acoustic. This problem remains under consideration from many years and has been investigated by many researchers. Many authors have investigated the nonlinear Schrödinger-KdV equation by various techniques such as the following: Wang [1] used finite difference method, Küçükarslan [2] used HPM, Bai and Zhang [3] used quadratic B-Spline finite element method, Fan and Hon [4] used extended tanh method, Kaya and El-Sayed [5] used adomian decomposition method (ADM), Doosthoseini and Shahmohamadi [6] used VIM, Alomari et al. [7] used homotopy analysis method (HAM), Qing et al. [8] used element free Galerkin method (EFG), and Golbabai and Safdari-Vaighani [9] used meshless method using RBF collocation scheme. The perturbation methods like HPM required a small parameter and are difficult to determine.

Recently, Marinca et al. introduced OHAM [1014] 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 extend the OHAM formulation for a system of three partial differential equations and to apply the extended OHAM formulation to coupled nonlinear Schrödinger-KdV equation. In [1517] OHAM has been proved to be valuable for obtaining an approximate solution of ordinary/partial differential equations (O/PDEs). Before, this system of nonlinear partial differential equations (NPDEs) was not solved by OHAM. We have proved that extended OHAM is useful and reliable for NPDEs, showing its validity and great potential for the solution of transient physical phenomenon in science and engineering.

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

#### 2. Extended Mathematical Formulation of OHAM

Consider a system of three partial differential equations: where , , are differential operators, , , are unknown functions, and denote spatial and temporal independent variables, respectively, is the boundary of , and , , are known analytic functions. , , can be divided into two parts: , , contain the linear parts while , , contain the nonlinear parts of the system of partial differential equations.

According to OHAM, we construct satisfying the following homotopies: where the auxiliary functions , , are nonzero for and , , . Equation (4) is called optimal homotopy equation. Clearly, we have

Obviously, when and we obtain respectively. When varies from to , the solution , , approaches from , , to , , , where , , are obtained from (4) for :

We choose auxiliary functions , , in the form

To get the approximate solutions, we expand , , by Taylor’s series about in the following manner: where . Now substituting (8)-(9) into (4) and equating the coefficient of like powers of , we obtain zeroth order system, given by (7) and the first and second order systems given by (10)-(11), respectively, and the general governing equations for are given by (12):

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

Substituting (14) into (1.1), the following expression for residuals results:

If , , then , , and will be the exact solutions of the problem. Generally it does not happen, especially in nonlinear problems.

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

The th order approximate solution can be obtained by these constants so-obtained. The more general auxiliary functions , , are useful for convergence, which depends upon constants , , , , , , can be optimally identified by (17), and is useful in error minimization.

#### 3. Application of Extended OHAM to Coupled Schrödinger-KdV Equation

To demonstrate the effectiveness of the extended OHAM for coupled Schrödinger- KdV equation taken from [6], we have with boundary conditions

The exact solution of (19) for and is given by

Applying the extended OHAM technique discussed in Section 2,

We consider Zeroth Order System. Consider with initial conditions

Its solution is First Order System. Consider with

Its solution is Second Order System. Consider with

The solution of second order system is

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

For the calculation of the constants , , , , and using (33) in (19) and applying the method of least square mentioned in (16)-(17) by taking, we get

The solution obtained by variational iteration method [6] is as follows:

The solution obtained by homotopy perturbation method [2] is as follows:

#### 4. Results and Discussions

The extended mathematical theory presented in Section 2 gives highly accurate solutions for the BVP presented in Section 3. We have used Mathematica 7 for most of our computational work. In Tables 1 and 2, the absolute errors of extended OHAM results corresponding to the exact solution of the problem given in (19) are given for , at , , , , . In Table 3, the extended OHAM results are compared with exact solution for . In Figures 1, 2, 3, 4, 5, and 6, 3D plots of the results obtained by extended OHAM and exact solutions are plotted for , , at fixed value of . While the 3D plots of the residuals at are given in Figures 7, 8, and 9, Figures 10, 11, 12, 13, 14, and 15 show the 2D plots of extended OHAM and exact solutions for , , at . The 2D residual plots are given in Figures 16, 17, and 18 at . In Figures 19, 20, and 21 the 2D plots of ordered solutions are given at for , , . The extended OHAM solutions are compared with VIM and HPM solution in 2D, Figures 22, 23, and 24. Here we observed from Figures 1921 that the extended OHAM solution converges rapidly with the increase in the order of approximation and produced better results. From Figures 2224 it is evident that the extended OHAM results are nearly identical to the exact, VIM, and HPM solutions. Here the results are very consistent with the decreasing time as evident from Tables 1-2.

Table 1: Absolute error of OHAM solution corresponding to the exact solution.
Table 2: Absolute error of OHAM solution corresponding to the exact solution.
Table 3: Comparison of solutions obtained by OHAM to the exact solution.
Figure 1: 3D, OHAM solution of at .
Figure 2: 3D, exact solution of at .
Figure 3: 3D, OHAM solution of at .
Figure 4: 3D, exact solution of at .
Figure 5: 3D, OHAM solution of at .
Figure 6: 3D, exact solution of at .
Figure 7: 3D, residual of at .
Figure 8: 3D, residual of at .
Figure 9: 3D, residual of at .
Figure 10: 2D, OHAM solution of at .
Figure 11: 2D, exact solution of at .
Figure 12: 2D, OHAM solution of at .
Figure 13: 2D, exact solution of at .
Figure 14: 2D, OHAM solution of at .
Figure 15: 2D, exact solution of at .
Figure 16: 2D, residual of at .
Figure 17: 2D, residual of at .
Figure 18: 2D, residual of at .
Figure 19: 2D, comparison of ordered solutions for at .
Figure 20: 2D, comparison of ordered solutions for at .
Figure 21: 2D, comparison of ordered solutions for at .
Figure 22: 2D, comparison of OHAM, VIM, HPM, and exact solutions for at .
Figure 23: 2D, comparison of OHAM, VIM, HPM, and exact solutions for at .
Figure 24: 2D, comparison of OHAM, VIM, HPM, and exact solutions for at .

#### 5. Conclusion

In this paper, we have seen the effectiveness of extended OHAM to coupled Schrödinger-KdV equation. By applying the basic idea of extended OHAM to coupled Schrödinger-KdV equation, we found it simpler in applicability and more convenient to control convergence and that it involved less computational overhead. Therefore, extended OHAM shows its validity and great potential for the solution of nonlinear system of PDEs 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

1. H. Wang, “Numerical studies on the split-step finite difference method for nonlinear Schrödinger equations,” Applied Mathematics and Computation, vol. 170, no. 1, pp. 17–35, 2005.
2. S. Küçükarslan, “Homotopy perturbation method for coupled Schrödinger-KdV equation,” Nonlinear Analysis: Real World Applications, vol. 10, no. 4, pp. 2264–2271, 2009.
3. D. Bai and L. Zhang, “Numerical studies on a novel split-step quadratic B-spline finite element method for the coupled Schrödinger-KdV equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 3, pp. 1263–1273, 2011.
4. E. Fan and Y. C. Hon, “Applications of extended tanh method to “special” types of nonlinear equations,” Applied Mathematics and Computation, vol. 141, no. 2-3, pp. 351–358, 2003.
5. D. Kaya and S. M. El-Sayed, “On the solution of the coupled Schrödinger-KdV equation by the decomposition method,” Physics Letters A, vol. 313, no. 1-2, pp. 82–88, 2003.
6. A. Doosthoseini and H. Shahmohamadi, “Variational iteration method for solving coupled Schrödinger-KdV equation,” Applied Mathematical Sciences, vol. 4, no. 17–20, pp. 823–837, 2010.
7. A. K. Alomari, M. S. M. Noorani, and R. Nazar, “Comparison between the homotopy analysis method and homotopy perturbation method to solve coupled Schrodinger-KdV equation,” Journal of Applied Mathematics and Computing, vol. 31, no. 1-2, pp. 1–12, 2009.
8. L. Y. Qing, C. R. Jun, and G. H. Xia, “An element free Galerkin (EFG) method for numerical solution of coupled Schrödinger-KdV equation,” Chinese Physics B, vol. 22, Article ID 100204, 2013.
9. A. Golbabai and A. Safdari-Vaighani, “A meshless method for numerical solution of the coupled Schrödinger-KdV equations,” Computing, vol. 92, no. 3, pp. 225–242, 2011.
10. 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.
11. 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.
12. 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.
13. V. Marinca, N. Herisanu, and I. Nemes, “A new analytic approach to nonlinear vibration of an electrical machine,” Proceedings of the Romanian Academy, vol. 9, pp. 229–236, 2008.
14. 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.
15. 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.
16. 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.
17. H. Ullah, S. Islam, M. Idrees, and R. Nawaz, “Application of optimal homotopy asymptotic method to doubly wave solutions of the coupled Drinfel d-Sokolov-Wilson equations,” Mathematical Problems in Engineering, vol. 2013, Article ID 362816, 8 pages, 2013.