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

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.


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]  (, )  + 3V (, ) V (, )  = 0, V (, )  + 2  3 V (, )  3   + 2 (, ) V (, )  +  (, )  V (, ) = 0, with  (, 0) = , V (, 0) = −. ( 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][15][16][17][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][20][21][22][23][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.

Fundamental Mathematical Theory of OHAM
Let us see the partial differential equation of the following form: where A 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.A can be divided into two parts L and N such that L is the simpler part of the partial differential equation which is easier to solve, and N contains the remaining part of A.
If (, ;   ) = 0, then ũ(, ;   ) is the exact solution of the problem.Generally it does not happen, especially in nonlinear problems.
For the determinations of auxiliary constants,   ,  = 1, 2, . . ., , 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  1 ,  2 , . .., can be optimally identified by (18), and is useful in error minimization.

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 Its solutions with For the calculations of the constants  11 ,  12 ,  21 , and  22 using (30) in ( 1) and applying the procedure mentioned in ( 16)-( 19), we get (31)

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 V(, ) at a spatial domain [0, 0.4] for  = 0.01,  = 0.015,  = 0.1,  = 0.2,  = 0.3, and  = 0.4, while in Tables 2 and 4 the convergence of the OHAM solution is given, through first and second order absolute errors at time  = 0.1 and 0 ≤  ≤ 1.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.

Conclusion
In this paper, we have seen the effectiveness of OHAM [16][17][18][19][20] in doubly periodic wave solutions of the coupled        shows its validity and great potential for the solution of time dependant problems in science and engineering.

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

Table 2 :
Comparison of first order and second order errors of  (, ) corresponding to the numerical solution at time  = 0.1, and 0 ≤  ≤ 1.

Table 3 :
Absolute error of OHAM solution of V (, ) corresponding to the numerical solution.

Table 4 :
Comparison of first order and second order errors of V (, ) corresponding to the numerical solution at time  = 0.1, and 0 ≤  ≤ 1.