Abstract Comparison of Optimal Homotopy Asymptotic Method with Homotopy Perturbation Method of Twelfth Order Boundary Value Problems MukeshGrover grover.mukesh@yahoo.com Department of Mathematics G.Z.S.C.E.T Bathinda, Dr.ArunKumarTomer 2 tomer4@rediffmail.com Department of Mathematics, S.M.D.R.S.D College Pathankot, Punjab, India In this work, we consider special problem consisting of twelfth order two-point boundary value by using the Optimal Homotopy Asymptotic Method and Homotopy Perturbation Method. Now, we discuss the comparison in between Optimal Homotopy Asymptotic Method and Homotopy Perturbation Method. These proposed methods have been thoroughly tested on problems of all kinds and shows very accurate results. A numerical eample is present and approimate is compared with eact solution and the error is compared with Optimal Homotopy Asymptotic Method and Homotopy Perturbation Method to assess the efficiency of the Optimal Homotopy Asymptotic Method at 2 th order Boundary values problems. Keywords Twelfth order boundary value problems, Approimate analytical solution, homotopy perturbation Method, optimal homotopy Asymptotic method, Ordinary Differential Equations, Error Estimates.. Introduction In literature different techniques are available for the numerical solution of twelfth order boundary value problems. The motivation of this problem is to etend Optimal Homotopy Asymptotic Method to solve linear and non linear tenth order boundary value problems. We also compared the results obtained from these techniques with the available eact solution in different literatures. Some properties of solutions of a given differential equation may be determined without finding their eact form in especially in nonlinear behavior. If as self-contained formula for the solution is not available, the solution may be numerically approimated using computers. To overcome these difficulties, a modified form of the perturbation method called homotopy perturbation method.hpm has since then been effectively utilized in obtaining approimate analytical solutions to many linear and nonlinear problems arising in engineering and science, such as nonlinear oscillators with discontinuities [3], nonlinear Volterra Fredholm integral equations [], Twelfth order differential equations have several important applications in engineering. Solution of linear and nonlinear boundary value problems of twelfth-order was implemented by Wazwaz using Adomian decomposition method. Chandrasekhar [9] showed that when an infinite horizontal layer of fluid is put into rotation and simultaneously subjected to heat from below and a uniform magnetic field across the fluid in the same direction as gravity, instability will occur. Several researchers developed numerical techniques for solving twelfth order differential equations. The Adomian Decomposition Method [, 4], the Differential Transform Method [5], the Variational Iteration Method, the successive iteration, the splines [5, 6], the Homotopy Perturbation Method [7], the Homotopy Analysis Method etc Recently Vasile Marinca et al. [,2,4] introduced OHAM for approimate solution of nonlinear problems of thin film flow of a fourth grade fluid down a vertical cylinder. OHAM is straight forward, reliable and it does not need to look for h curves like HPM. Moreover, this method provides a convenient way to control the convergence of the series solution. Most recently, Javed Ali et al. used OHAM for the solutions of multi-point boundary value problems. The results of OHAM presented in this work are compared with those of eact solution HPM. ISSN : 975-3397 Vol. 3 No. 7 July 2 2739
2. Basic Idea of Homotopy-Perturbation Method To clarify the basic ideas of the homotopy perturbation method [3], we consider the following nonlinear differential equation with boundary conditions A (u) f (r) =, r Ω () B (u, ) =, r Γ (2) where A is a general differential operator, B is a boundary operator, u is a known analytical function, and Γ is the boundary of the domain Ω. The operator A can be divided into two parts L and N, where L is linear, while N is nonlinear. Therefore () can be rewritten as follows L (u) + N(u) f(r) =. (3) By the homotopy technique, we can construct a homotopy V (r, p): Ω [, ] R which satisfies H (v, p) = ( p) [L (v) L (u )] + p [A (v) f(r)] =, pϵ [, ] Ω (4) or H (v, p) = L (v) L (u ) + p L (u ) + p [N (v) f(r)] =, (5) where r Γ and p [, ] is an embedding parameter, u is an initial approimation of (), which satisfies the boundary conditions. Obviously, from Equations (4) and (5) we will have: H (v, ) = L (v) L (u ) =, (6) H (v, ) = A (v) f(r) =, (7) and the changing process of p from zero to unity is just that of H(v, p) from L(v) L(u ) to A(v) f(r). In topology, this is called deformation, L (v) L (u ) and A (v) f(r) is called homotopic. The embedding parameter p is introduced much more naturally, unaffected by artificial factors. Furthermore, it can be considered as a small parameter for < p. So it is very natural to assume that the solution of (4, 5) can be epressed as v = v + pv + p 2 v 2 +. Therefore, when p=, the approimate solution of above equation can be readily obtained as follows: u = lim v + v + v 2 +. p The combination of the perturbation method and the homotopy method is called the HPM, which eliminates the drawbacks of the traditional perturbation methods while keeping all its advantages. The above series is convergent for most cases. Using the transformation y = y, = y 2, = y 3,., = y 2n (8) Using the boundary conditions we can be written as a system of integral equations: ISSN : 975-3397 Vol. 3 No. 7 July 2 274
y = a + y, y 2 = b + y, y 3 = c + y 2n = d + f{t, y (t),y (t),y (t),.y (t)}dt 3. Optimal Homotopy Asymptotic Method Analysis y... The general differential equation is considered as, L (u()) +g() +N (u()), B (u, ) = (9) where L is a linear operator, denotes independent variable, u() is an unknown function, g() is a known function, N (u()) is a nonlinear operator and B is a boundary operator. A family of equations is constructed using (OHAM) which is given by ( p)l (f (, p) + g() = H (p) (L (f (, p) + g() + N (f (, p)) B (f (, p),f (, p) / ). () Where p [,] is an embedding parameter, H (p) is a nonzero auiliary function for p and H() =,f (, p) is an unknown function, respectively. Obviously, when p= and p =, it holds f (,) = u (), f (,) = u (), () respectively. Thus, as p increases from to, the solution f (, p) varies from u () to the solution u (), where u () is obtained from equation () for p=: L (u ()) +g() =, B (u, d u /d ) = (2) Auiliary function H (p) is chosen in the form H p = pc + p 2 c 2 +.. (3) where constants c, c 2,... can be determined latter by the method of least squares. Let us consider the solution of Eq. () in the form to get an approimate solution; we used Taylor s series to epand f(, p, C i ) in p which becomes f(, p, c ) = u () +u (, c,c,c, c )p (4) Now substituting Eq. (4) into Eq. () and equating the coefficients of like powers of p,, we obtain the following linear equations. Zeroth order problem is given by Eq. (2) and the first order problem is given by L (u ()) +g() = C N (u ()), B (u, du ) = (5) dn The mth -order problem is given by: L (u m ()) - L (u m- ()) = C m N (u () + m i c i [L u m ()+N m u (), u (), u 2 (),.,u m ()] (6) Here m varies 2, 3, 4..B (u m, du m dn ) = where Nm (u (), u (),, u m ()) is the coefficient of p m in the epansion of N(n(, p)) about the embedding parameter p. ISSN : 975-3397 Vol. 3 No. 7 July 2 274
Nf(, p, c i ) = N u ()+N m, c, c 2, c 3, c m p m It has been observed that the convergence of the series (4) depends upon the auiliary constants C, C 2,.. If it is convergent at p =, one has m f(, c i ) = u () +u m (c, c 2, c 3, c m ) m The result of the mth-order approimations are given uc, c 2, c 3, c m = u ()+u i, c, c 2, c 3, c i Substituting Eq. (7) into Eq. (9), it results the following residual: m i (7) R, c, c 2, c 3, c m = L(u, c, c 2, c 3, c m ) + g()+nu(, c, c 2, c 3, c m ) If R =, then u% will represent the eact solution. Generally it doesn t happen, especially in nonlinear problems. To have a minimum value of R for the values, C i here i =, 2, 3 various methods for instance the weighted residual method or variational methods etc. can be used. By using the Least squares method we obtain the following equation b R a, c, c 2, c 3, c m R c i d = (8) while the Galerkin s method gives b R a, c, c 2, c 3, c m R c i d = (9) where a and b are in the domain of the problem. For the given value of all constants we can easily determine the approimate solution of order m. 4. NUMERICAL EXAMPLES Eample : Consider the following linear twelfth order boundary value problem y 2 () = - y()- 3 e -23e -2e with following conditions: =, () =, (2) =, (3) = 3, (4) = 8, (5) = 5 =, () = e, (2) = 4e, (3) = 9e, (4) = 6e, (5) = 25e (2) Eact solution is y () = (- ) e Using the transformation (8) we can rewrite the twelve-order boundary value problem (2) as the system of integral equations: y = + y 2 = + y 3 = + y 2 y 3 y 4 ISSN : 975-3397 Vol. 3 No. 7 July 2 2742
y 4 = -3 + y 5 = -8 + y 6 = -5 + y 7 = a + y 8 = b + y 9 = c + y = d + y = e + y 5 y 6 y 7 y 8 y 9 y y y 2 y 2 = f + { y (t) 2e t 23 te t -t 3 e t }dt Comparing the coefficients of like powers of p, we have: Coefficients of p : y = y 2 = y 3 = y 4 = -3 y 5 = -8 y 6 = -5 y 7 = a y 8 = b y 9 = c y = d y = e y 2 = f Coefficients of p : y = y 2 = y 3 = -3 y 4 = -3-8 y 5 = -8-5 y 6 = -5+a Y 7 = a+ b Y 8 = b +c y 9 = c+ d y = d +e Y = e +f Y 2 = f-e ( 3-3 2 + 9+29) Coefficients of p 2 : Coefficients of p 3 : y 2 = y 22 = - 3 2 2 y 32 =-3-2 y 42 = -3-8- 5 2 2 y 52 =-8-5- a 2 2 y 62 = - 5+ a + b 2 2 y 72 = a+ b + c 2 2 y 82 = b+ c + d 2 2 y 92 = c+ d + e 2 2 Y 2 = d+ e + f 2 2 Y 2 = e + 5 +(9 +f) e (5+4-6 2 + 3 ) Y 22 = f +9-3 3 e (5+29-3 2 + 3 ) y 3 = - 3 2 y 23 = - 3 2 2 3 2 y 33 =-3-2 2-5 2 3 y 43 = -3-8- 5 2 2 + a 6 3 y 53 =-8-5- a 2 2 + b 6 3 y 63 = - 5+ a + b 2 2 + c 6 3 y 73 = a+ b + c 2 2 + d 6 3 y 83 = b+ c + d 2 2 + f 6 3 y 93 = c+ d + e 2 2 + f 6 3 Y 3 = -9 +d +(5 +e) +.5(9 +f) 2 e (-9+59-9 2 + 3 ) Y 3 = e + 5 +(9 +f) e (5+4-6 2 + 3 ) - 4 2 ISSN : 975-3397 Vol. 3 No. 7 July 2 2743
Y 23 = f +9-3 3 e (9+29-3 2 + 3 ) y 4 = - 3 2-4 3 Coefficients of p 4 : y = - 3 2 4 2 3 2-54 8 y 34 =-3-2 2-5 2 3 + a4 y 44 = -3-8- 5 2 2 + a 6 3 + b4 y 54 =-8-5- a 2 2 + b 6 3 + c4 y 64 = - 5+ a + b 2 2 + c 6 3 + d4 y 74 = a+ b + c 2 2 + d 6 3 + f 4 y 84 = b+ c + d 2 2 + e 6 3 + f 4 y 94 = -92 +c +(-9 +d) +.5(5 +f) 2 e (-92+83-2 2 + 3 )+(9 +f) 6 3 y 4 = -9 +d +(5 +e) +.5(9 +f) 2 e (-9+59-9 2 + 3 ) - 5 y 4 = e + 5 +(9 +f) e (5+4-6 2 + 3 ) - 4 y = f +9-3 3 e (9+29-3 2 + 3 ) + 5 6 2 6 Above procedure is proceeding up to coefficient of p 2.Now, adding above all coefficient of p to p 2, and then we obtain: y 2 () = - 2 3-3 4-8 5 + a 72 6 + b 54 7 + c 432 8 + d 36288 9 + e 36288 + f 39968-39968 2 435456 3-83 43589456 4 6762 5 + The coefficients a, b, c, d, e, f can be obtained using the boundary conditions at =, a = 23.9999985, b =35.57, c = 47.99896, d =63.83, e =79.935948, f= 99.737663. The series solution can, thus, be written as y 2 () = - 2 3-3 4-8 5 +.333333 6 -.694446 7 -.945 8 -.7364 9.22282 +. 845 435456 3 83-43589456 4 6762 5 + 39968 2 ISSN : 975-3397 Vol. 3 No. 7 July 2 2744
Analytical Solution Numerical Solution Errors......994653826.994653826 3. -.2.954443.954443..3.283473497.283473496 -. -.4.358379275.358379277 2. -.5.428378.428389. -9.6.4373852.43738564 4.4-9.7.422888685.42288882.35-8.8.356865485.356865853 3.68-8.9.2236428.2236437 9. -8...227 2.27-7 Table: [.] Eample 2: Consider the following linear twelfth order boundary value problem y 2 () = - y()- 3 e -23e -2e with following conditions: =, () =, () =, () = 3, () = 8, () = 5 =, () = e, () = 4e, () = 9e, () = 6e, () = 25e Eact solution is y () = (- ) e We construct the following zeroth and first-order problems. y () = - y()- 3 e -23e -2e with following conditions: () =, () () =, () () =, () () = 3, () () = 8, () () = 5 () =, () () = e, () () = 4e, () () = 9e, () () = 6e, () () = 25e First-Order Problem: y () = (+ C ) (2 + 23 + 3 ) e + C y () + ( + C ) y () With same above boundaries conditions Solutions to these problems are given by Equations. (2) and (2) respectively y () = (-288 +288 e -2276-5892 e -858 2 +432e 2-26 3-2 e 3-396 4-56 5 +977 6 +727376e 6-8368968 +3787725e 7 + 458966e 8 ISSN : 975-3397 Vol. 3 No. 7 July 2 2745
-28234 9 + 477329e 9 + 57487-2229665e -39868 +3825434e ) (2) y (, C ) = C (266-266e + (76352 3978 e ) + (745-2723 e ) 2 + (8795+ 84 e ) 3 + (3663 - e ) 4 + 558.833 5 + 69.47 6 +7.7738 7 +.6367 8 +.425263 9 + +.237265 +.89239-3.75782-7 3-4.23959-8 4-3.2863-9 5-2.6473-6 -.334-7 -5.85-2 8-2.758-4 9-7.898-6 2-2.4669-7 2-6.459-9 22 -.27589-2 23 -.5569-22 ) (2) Considering the OHAM first-order solution, Y app (, C ) = y () + y (, C ) (22) and using Eq.(8) with a =.5 and b =, we get C = -.2647. Using this value the first- order solution (22) is well-determined. Analytical Solution Numerical Solution Errors......99465383.99465383-7.565-4.2.95444.95444-2.7686-2.3.2834735.2834735 -.727 -.4.35837927.35837927-5.232 -.5.42838.42838-9.34 -.6.4373852.4373852 -.279 -.7.42288868.42288869 -.397 -.8.35686548.35686549 -.2278 -.9.2236428.2236428-7.4997 -.. -..9454-5. Conclusion Table: [. 2] In this paper, the Comparison of the results obtained by the Homotopy perturbation method and optimal homotopy asymptotic method of Twelfth order boundary value problems. The numerical results in the Tables [.-.2], show that the optimal homotopy asymptotic method provides highly accurate numerical results as compared to Homotopy perturbation method. It can be concluded that optimal homotopy asymptotic method is a highly efficient method for solving 2 th order boundary value problems arising in various fields of engineering and science. References [] Adomian, G., 99. A review of the decomposition method and some recent results for Nonlinear equations. Comput. Math. Appl., -27. [2] He, J.H. (). The Homotopy Perturbation Method for Non-Linear Oscillators with Discontinuities. Applied Maths. Computer. 287 292. [3] He, J.-H., 23. Variational approach to the sith order boundary value problems. Applied Math. Computer, 537-538. [4] Wazwaz, A.M., 2. Approimate solutions to boundary value problems of higher- order By the modified decomposition method. Computer.Math.Appl.679-69 [5] Siddique, S.S. and E.H. Twizell, 997. Spline solutions of linear twelfth-order boundary Value Problems Computer Math. Appl., 37-39. [6] Siddique, S.S. and Ghazala Akram, 28. Solutions of 2th order boundary value problems Using non-polynomial spline technique. Appl. Math.Comput.,559-57. [7] He, J.-H., 26. Homotopy perturbation methodfor solving boundary value problems.87-88. ISSN : 975-3397 Vol. 3 No. 7 July 2 2746
[8] Ghotbi, Abdoul R., Barari, A. and Ganji, D. D, (28). Solving Ratio-Dependent Predator-Prey System with Constant Effort Harvesting Using Homotopy Perturbation Method. Math. Prob.Eng., -8. [9] Chandrasekhar, S. (96). Hydrodynamic and Hydromagnetic Stability. Clarendon Press, Oford, Reprinted: Dover Books, New York, 98. [] Marinca, V., N. Herisanu and I. Nemes, 28.Optimal homotopy asymptotic method with Application to thin film flow. Cent. Eur. J. Phys. 648-653. [] Ghasemi, M., Tavassoli Kajani,M. and Babolian, E. (27). Numerical Solution of the Nonlinear Volterra Fredholm Integral Equations by using Homotopy Perturbation Method. Appl. Math. Computer 446-449. [2] Marinca, V. and N. Herisanu, 28. An optimal homotopy asymptotic method for solving Non linear equations arising in heat transfer. Int. Comm. Heat and Mass Tran. 7-75. [3] S. H. Mirmoradi, S. Ghanbarpour, I. Hosseinpour, A. Barari, (29). Application of Homotopy perturbation method and variational iteration method to a nonlinear fourth order Boundary value problem. Int. J.Math.Analysis., -9. [4] Marinca, V., N. Herisanu, C. Bota and B. Marinca,29. An optimal homotopy asymptotic Method applied to the steady flow of fourth-grade fluid pasta porous plate. App.Math.5-5 [5] Islam, S.-U., S. Haq and J. Ali, 29. Numerical solution of special 2th-order boundary Value Problems using differential transform method. Comm. Nonl. Sc. Nu. Sim.32-38. ISSN : 975-3397 Vol. 3 No. 7 July 2 2747