Adv. Theor. Appl. Mech., Vol. 3, 21, no. 11, 513-52 An Analytic Study of the (2 + 1)-Dimensional Potential Kadomtsev-Petviashvili Equation B. Batiha and K. Batiha Department of Mathematics, Faculty of Science Qassim University, Kingdom of Saudi Arabia belalbatiha22@yahoo.com Department of Business Networking and Systems Management Philadelphia University, Jordan kh batiha@philadelphia.edu.jo Abstract In this paper, variational iteration method (VIM) is applied to obtain approximate analytical solution of the (2 + 1)-dimensional potential Kadomtsev-Petviashvili equation (PKP) without any discretization. Comparisons with the exact solutions reveal that VIM is very effective and convenient. Keywords: Variational iteration method (VIM); (2 + 1)-dimensional potential Kadomtsev-Petviashvili equation (PKP); Soliton-like solutions 1 Introduction We consider the (2+1)-dimensional potential Kadomtsev-Petviashvili equation (PKP) in the form u xt + 3 2 u xu xx + 1 4 u xxxx + 3 4 u yy =, (1) where the initial conditions u(x,,t) and u y (x,,t) are given. Nonlinear phenomena are of fundamental importance in various fields of science and engineering. The nonlinear models of real-life problems are still difficult to solve either numerically or theoretically. There has recently been much attention
514 B. Batiha and K. Batiha devoted to the search for better and more efficient solution methods for determining a solution, approximate or exact, analytical or numerical, to nonlinear models [1, 2, 3]. Senthilvelan [4] solved PKP using the homogenous balance method (HBM). Li and Zhang [5] made some improvement on the key steps of HBM, and applied it to the research of PKP equation. Various exact solutions, include soliton, multisoliton and rational-type solutions, of the equation are obtained. Li and Zhang [6] obtained new soliton-like solutions for (2 + 1)-dimensional (PKP) equation using the symbolic computation method. Kaya and El-Sayed [7] used Adomian decomposition method (ADM) to solve (2 + 1)-dimensional (PKP) equation. The variational iteration method (VIM) is a simple and yet powerful method for solving a wide class of nonlinear problems, first envisioned by He [8] (see also [9, 1, 11, 12, 13]). The VIM has successfully been applied to many situations. For example, He [9] solved the classical Blasius equation using VIM. He [1] used VIM to give approximate solutions for some well-known nonlinear problems. He [11] used VIM to solve autonomous ordinary differential systems. He [12] coupled the iteration method with the perturbation method to solve the well-known Blasius equation. He [13] solved strongly nonlinear equations using VIM. Soliman [14] applied the VIM to solve the KdV-Burger s and Lax s seventh-order KdV equations. The VIM has recently been applied for solving nonlinear coagulation problem with mass loss by Abulwafa et al. [15]. The VIM has been applied for solving nonlinear differential equations of fractional order by Odibat et al. [16]. Bildik et al. [17] used VIM for solving different types of nonlinear partial differential equations. Dehghan and Tatari [18] employed VIM to solve a Fokker-Planck equation. Wazwaz [19] presented a comparative study between the variational iteration method and Adomian decomposition method. Tamer et al. [2] introduced a modification of VIM. Batiha et al. [21] used VIM to solve the generalized Burgers-Huxley equation. Batiha et al. [22] applied VIM to the generalized Huxley equation. Abbasbandy [23] solved one example of the quadratic Riccati differential equation (with constant coefficient) by He s variational iteration method by using Adomian s polynomials. In this paper, we shall apply VIM to find the approximate analytical solution of PKP equation. Comparisons with the exact solution shall be performed.
VIM for PKP 515 2 Variational iteration method VIM is based on the general Lagrange s multiplier method [25]. The main feature of the method is that the solution of a mathematical problem with linearization assumption is used as initial approximation or trial function. Then a more highly precise approximation at some special point can be obtained. This approximation converges rapidly to an accurate solution [24]. To illustrate the basic concepts of VIM, we consider the following nonlinear differential equation: Lu + Nu = g(x), (2) where L is a linear operator, N is a nonlinear operator, and g(x) is an inhomogeneous term. According to VIM [1, 11, 13, 24], we can construct a correction functional as follows: u n+1 (x) =u n (x)+ x λ{lu n (τ)+nũ n (τ) g(τ)}dτ, n, (3) where λ is a general Lagrangian multiplier [25] which can be identified optimally via the variational theory, the subscript n denotes the nth-order approximation, ũ n is considered as a restricted variation [1, 11], i.e. δũ n =. 3 Analysis of the PKP equation In this section, we present the solution of Eq. (1) by means of VIM. To do so, we first construct a correction functional, u n+1 = u n + λ(s) [(u n ) ss + 4 3 (ũ n) xt +2(ũ n ) x (ũ n ) xx + 1 ] 3 (ũ n) xxxx ds, (4) where ũ n is considered as restricted variations, which means δũ n =. To find the optimal λ(s), we proceed as follows: δu n+1 = δu n + δ λ(s) [(u n ) ss + 4 3 (ũ n) xt +2(ũ n ) x (ũ n ) xx + 1 ] 3 (ũ n) xxxx ds, (5) and consequently δu n+1 = δu n + δ λ(s)[(u n ) ss ]ds, (6) which results in δu n+1 = δu n (1 λ (s)) + δ(u n ) s λ(s)+ δu n λ (s)ds =. (7)
516 B. Batiha and K. Batiha The stationary conditions can be obtained as follows: 1 λ (s) = s=y, λ(s) = s=y, λ (s) = s=y, (8) The Lagrange multipliers, therefore, can be identified as λ(s) =s y, (9) and the iteration formula is given as u n+1 = u n + (s y) [(u n ) ss + 4 3 (u n) xt +2(u n ) x (u n ) xx + 1 ] 3 (u n) xxxx ds, (1) 4 Applications We first consider the PKP equation (1) with initial conditions u(x,,t)=1+2kα tanh[k(αx ct)], u y (x,,t)=2αβk 2 sech 2 [k(αx ct)],(11) where c =(k 2 α 3 + 3β3 ). The exact solution for Eq. (1) was found to be [4]: [ ( ( ) )] u(x, y, t) =1+2kα tanh k αx + βy k 2 α 3 + 3β3 t (12) We can take an initial approximation as u = u(x,,t)+yu y (x,,t), (13) = 1+2kα tanh[k(αx ct)] + 2αβk 2 y sech 2 [k(αx ct)] (14) where c =(k 2 α 3 + 3β3 ). The first iterate is easily obtained from (1) and is given by: u[1] = 1 + 2kαm +2αβk 2 yn 2 2α 2 βk 3 n 2 m(8α 3 βk 4 n 2 m 2 (15) 3 βk 4 n 2 (1 m 2 ))y 4 +1/3(8yα 2 βk 3 n 2 m(8α 3 βk 4 n 2 m 2 3 βk 4 n 2 (1 m 2 ))+32/3α 2 βk 4 n 2 m 2 ( k 2 α 3 3/4β 3 /α) +4k 2 α 2 (1 m 2 )(8α 3 βk 4 n 2 m 2 3 βk 4 n 2 (1 m 2 )) 8/3α 5 βk 6 n 2 m 2 (1 m 2 ) 16/3α 2 βk 4 n 2 (1 m 2 )( k 2 α 3 3/4β 3 /α) +32/3α 5 βk 6 n 2 m 4 +32/3α 5 βk 6 n 2 (1 m 2 ) 2 )y 3 +1/2( y(32/3α 2 βk 4 n 2 m 2 ( k 2 α 3 3/4β 3 /α) +4k 2 α 2 (1 m 2 )(8α 3 βk 4 n 2 m 2 3 βk 4 n 2 (1 m 2 ))
VIM for PKP 517 u(x, y, t) 1.2 1.15 1.1 1.5 1.95.9.85.8-1 -5 x 5.2.15.1.5 -.5 y -.1 -.15 1 -.2 1.2 1.15 1.1 1.5 1.95.9.85.8 1.2 1.15 1.1 1.5 1.95.9.85.8 u(x, y, t) -1-5 x 5.2.15.1.5 -.5 y -.1 -.15 1 -.2 1.2 1.15 1.1 1.5 1.95.9.85.8 Figure 1: ccomparison between exact solution and the numerical results for u(x, y, t) by means of 1-iterate VIM solution when α =1,β =.1,k =.1 and c = k 2 α 3 + 3β3 8/3α 5 βk 6 n 2 m 2 (1 m 2 ) 16/3α 2 βk 4 n 2 (1 m 2 )( k 2 α 3 3/4β 3 /α) +32/3α 5 βk 6 n 2 m 4 +32/3α 5 βk 6 n 2 (1 m 2 ) 2 ) 16/3k 3 α 2 m(1 m 2 )( k 2 α 3 3/4β 3 /α) 16/3k 5 α 5 (1 m 2 ) 2 m 16/3k 5 α 5 m 3 (1 m 2 ))y 2 y 2 ( 16/3k 3 α 2 m(1 m 2 )( k 2 α 3 3/4β 3 /α) 16/3k 5 α 5 (1 m 2 ) 2 m 16/3k 5 α 5 m 3 (1 m 2 )) where m = tanh[k(αx ct)] n = sech [k(αx ct)] Fig. 1 shows the comparison between the 1-iterate of VIM and the exact solution (12). Now, we will find the approximate analytical solution of the PKP equation (1) with initial conditions u(x,,t)=1+2kα tan[k(αx ct)], u y (x,,t)=2αβk 2 sec 2 [k(αx ct)], (16) where c =(k 2 α 3 + 3β3 ). The exact solution for Eq. (1) is [4]: [ ( ( ) )] u(x, y, t) =1+2kα tan k αx + βy + k 2 α 3 3β3 t (17) We can take an initial approximation u = u(x,,t)+yu y (x,,t), (18) = 1+2kα tan[k(αx ct)] + 2αβk 2 y sec 2 [k(αx ct)] (19) The first iterate is easily obtained from (1) and is given by: u[1] = 1 + 2kαm +2αβk 2 yn 2 +2α 2 βk 3 n 2 m(8α 3 βk 4 n 2 m 2 (2) + 3 βk 4 n 2 (1 + m 2 ))y 4 +1/3( 8yα 2 βk 3 n 2 m(8α 3 βk 4 n 2 m 2
518 B. Batiha and K. Batiha Table 1: Absolute errors between the 1-iterate of VIM and the exact solutions (17) when α =1,β =.1,k =.1 and c = k 2 α 3 + 3β3 x i /y i.1.2.3.4.5.1 1.88111E 8 7.76727E 8 1.8227E 7 3.3118E 7 5.3987E 7.2 3.745E 8 1.5595E 7 3.44314E 7 6.21848E 7 9.86846E 7.3 5.53192E 8 2.23718E 7 5.8858E 7 9.1444E 7 1.4442E 6.4 7.36714E 8 2.97138E 7 6.748E 7 1.2818E 6 1.9311E 6.5 9.21218E 8 3.7955E 7 8.4199E 7 1.5356E 6 2.36474E 6 + 3 βk 4 n 2 (1 + m 2 ))+32/3α 2 βk 4 n 2 m 2 c +4k 2 α 2 (1 + m 2 )(8α 3 βk 4 n 2 m 2 + 3 βk 4 n 2 (1 + m 2 )) +272/3α 5 βk 6 n 2 m 2 (1 + m 2 )+16/3α 2 βk 4 n 2 (1 + m 2 )c +32/3α 5 βk 6 n 2 m 4 +32/3α 5 βk 6 n 2 (1 + m 2 ) 2 )y 3 +1/2( y(32/3α 2 βk 4 n 2 m 2 c +4k 2 α 2 (1 + m 2 )(8α 3 βk 4 n 2 m 2 α 3 βk 4 n 2 (1 + m 2 )) +272/3α 5 βk 6 n 2 m 2 (1 + m 2 )+16/3α 2 βk 4 n 2 (1 + m 2 )c +32/3α 5 βk 6 n 2 m 4 +32/3α 5 βk 6 n 2 (1 + m 2 ) 2 )+16/3k 3 α 2 m(1 + m 2 )c +8/3k 5 α 5 (1 + m 2 ) 2 m +16/3k 5 α 5 m 3 (1 + m 2 ))y 2 y 2 (16/3k 3 α 2 m(1 + m 2 )c +8/3k 5 α 5 (1 + m 2 ) 2 m +16/3k 5 α 5 m 3 (1 + m 2 )) where c =(k 2 α 3 + 3β3 ), m = tan[k(αx ct)], n = sec[k(αx ct)]. Table 1 shows the comparison between the 1-iterate of VIM and the exact solution (17). 5 Conclusions In this paper, the variation iteration method (VIM) has been successfully employed to obtain the approximate analytical solutions of the (2+1)-dimensional potential Kadomtsev-Petviashvili equation (PKP). The method has been applied directly without using linearization or any restrictive assumptions. Comparisons with the exact solution reveal that VIM is very effective and convenient. It may be concluded that VIM is very powerful and efficient in finding analytical as well as numerical solutions for a wide class of linear and nonlinear differential equations. VIM provides more realistic series solutions that converge very rapidly in real physical problems.
VIM for PKP 519 References [1] X.Y. Wang, Exact and explicit solitary wave solutions for the generalized Fisher equation Phys. Lett. A 131 (1988) 277-279. [2] A. Jeffrey, M.N.B. Mohamad, Exact solutions to the KdV-Burgers equation Wave Motion 14 (1991) 369 375. [3] M. Wadati, The exact solution of the modified Korteweg-de Vries equation J. Phys. Soc. Jpn. 32 (1972) 1681 1687. [4] M. Senthilvelan, On the extended applications of homogenous balance method, Appl. Math. Comput. 123 (21) 381 388. [5] D.S. Li, H. Q. Zhang, Symbolic computation and various exact solutions of potential Kadomstev-Petviashvili equation, Appl. Math. Comput. 145 (23) 351 359. [6] D.S. Li, H. Q. Zhang, New soliton-like solutions to the potential Kadomstev-Petviashvili (PKP) equation, Appl. Math. Comput. 146 (23) 381 384. [7] D. Kaya, S. M. El-Sayed, Numerical soliton-like solutions of the potential Kadomtsev-Petviashvili equation by the decomposition method, Phys Lett A 32 (23) 192-199. [8] J.H. He, A new approach to nonlinear partial differential equations. Commun. Nonlin. Sci. Numer. Simul. 2 (1997) 23-235. [9] J.H. He, Approximate analytical solution of Blasius equation. Commun. Nonlin. Sci. Numer. Simul. 3 (1998) 26 263. [1] J.H. He, Variational iteration method a kind of non-linear analytical technique: some examples. Int. J. Non-Linear Mech. 34 (1999) 699 78. [11] J.H. He, Variational iteration method for autonomous ordinary differential systems. Appl. Math. Comput. 114 (2) 115-123. [12] J.H. He, A simple perturbation approach to Blasius equation. Appl. Math. Comput. 14 (23) 217 222. [13] J.H. He, Some asymptotic methods for strongly nonlinear equations. Int. J. Modern Phys. B 2 (26) 1141 1199. [14] A.A. Soliman, A numerical simulation and explicit solutions of KdV- Burgers and Lax s seventh-order KdV equations. Chaos Solitons Fractals 29(2) (26) 294 32.
52 B. Batiha and K. Batiha [15] E.M. Abulwafa, M.A. Abdou and A.A. Mahmoud, The solution of nonlinear coagulation problem with mass loss. Chaos Solitons Fractals 29(2) (26) 313 33. [16] Z.M. Odibat, S. Momani, Application of variational iteration method to nonlinear differential equations of fractional order. Int. J. Nonlinear Sci. Numer. Simul. 7 (26) 27 34. [17] N. Bildik and A.Konuralp, The use of variational iteration method, differential transform method and Adomian decomposition method for solving different types of nonlinear partial differential equations. Int. J. Nonlinear. Sci. Numer. Simul. 7 (26) 65 7. [18] M. Dehghan and M. Tatari, The use of He s variational iteration method for solving a Fokker-Planck equation. Phys. Scr. 74 (26) 31 316. [19] A.M. Wazwaz, A comparison between the variational iteration method and Adomian decomposition method. J. Comput. Appl. Math. doi:1.116/j.cam.26.7.18. [2] H. El-Zoheiry, T.A. Abassy and M.A. El-Tawil, Toward a modified variational iteration method. J. Comput. Appl. Math. doi:1.116/j.cam.26.7.19. [21] B. Batiha, M.S.M. Noorani, I. Hashim, Application of variational iteration method to the generalized Burgers-Huxley equation Chaos Solitons Fractals doi: 1.116/j.chaos.26.6.8, in press [22] B. Batiha, M.S.M. Noorani, I. Hashim, Numerical simulation of the generalized Huxley equation by He s variational iteration method Appl. Math. Comput. 186 (27) 1322 1325. [23] S. Abbasbandy, A new application of He s variational iteration method for quadratic Riccati differential equation by using Adomian s polynomials. J. Comput. Appl. Math. doi:1.116/j.cam.26.7.12. [24] J.H. He, Non-perturbative methods for strongly nonlinear problems, Berlin: dissertation.de-verlag im Internet GmbH, 26. [25] M. Inokuti, H. Sekine and T. Mura, General use of the Lagrange multiplier in nonlinear mathematical physics. In: Nemat-Nassed S, editor. Variational method in the mechanics of solids. Pergamon Press, 1978. Received: August, 29