• Views 432
• Citations 1
• ePub 15
• PDF 548
`Journal of Applied MathematicsVolume 2013 (2013), Article ID 468909, 8 pageshttp://dx.doi.org/10.1155/2013/468909`
Research Article

## Numerical Scheme for Solving Singular Two-Point Boundary Value Problems

1School of Mathematical Sciences, Universiti Kebangsaan Malaysia, 43600 Bangi, Selangor, Malaysia
2Department of Mathematics, Faculty of Science, Hashemite University, 13115 Zarqa, Jordan

Received 25 December 2012; Revised 5 March 2013; Accepted 11 March 2013

Copyright © 2013 N. Ratib Anakira 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

Singular two-point boundary value problems (BVPs) are investigated using a new technique, namely, optimal homotopy asymptotic method (OHAM). OHAM provides a convenient way of controlling the convergence region and it does not need to identify an auxiliary parameter. The effectiveness of the method is investigated by comparing the results obtained with the exact solution, which proves the reliability of the method.

#### 1. Introduction

Consider singular two-point boundary value problems (BVPs) of the form subject to the boundary conditions where , and are continuous functions on and the parameters , and are real constants.

Problems of the forms (1) and (2) are encountered in the fields of fluid mechanics, reaction-diffusion processes, chemical kinetics, optimal control, and other branches of applied mathematics [1, 2]. Many different numerical methods have been proposed by various authors regarding singular two-point boundary value problems such as variational iteration method (VIM) [3], cubic splines [4], differential transformation method (DTM) [5], the Adomian decomposition method (ADM) [6], continuous genetic algorithm (CGA) [7], sinc-Galerkin method and homotopy perturbation method (HPM) [813], and homotopy analysis method (HAM) [1423]. The existence of a unique solution of (1) and (2) was discussed in [24, 25].

Recently, Marinca et al. introduced and developed optimal homotopy asymptotic method (OHAM) in a series of papers [2629] for the approximate solutions of nonlinear problems. OHAM does not depend on the presence of a small parameter. An advantage of OHAM is that it does not need to identify the -curve as in HAM for convergence region. In OHAM, the control and adjustment of the convergence region is provided in a convenient way. Furthermore, it has a built-in convergence criteria similar to HAM but with a greater degree of flexibility.

In this paper, OHAM is presented to create an approximate analytic solution of singular two-point boundary value problems. The method is directly applied without any linearization and discretizations and without splitting the nonhomogeneous term. The structure of this paper is organized as follows: Section 2 is devoted to the analysis of the proposed method, in Section 3, three examples are employed to illustrate the accuracy and computational efficiency of this approach, and lastly, conclusions are given in the last section.

#### 2. Analysis of the Method

To illustrate the basic idea of OHAM [30], we consider the following differential equation: where is the chosen linear operator, is the linear or nonlinear operator, is an unknown function, denotes an independent variable, is a known function, and is a boundary operator.

According to the basic idea OHAM we construct a homotopy which satisfies where and is an embedding parameter, is a nonzero auxiliary function for , , and is an unknown function. Obviously, when and it holds that = and , respectively. Thus, as varies from to , the solution approaches from to where is the initial guess that satisfies the linear operator and the boundary conditions Next, we choose the auxiliary function in the form where are constants which can be determined later. can be expressed in many forms as reported by Marinca et al. [2629].

To get an approximate solution, we expand in Taylor’s series about in the following manner: Substituting (7) into (4) and equating the coefficient of like powers of , we obtain the following linear equations. The zeroth-order problem is given by (5); the first- and second-order problems are given as The general governing equations for are where and is the coefficient of in the expansion of about the embedding parameter : It has been observed that the convergence of the series (7) depends upon the auxiliary constants . If it is convergent at , one has The result of the th-order approximation is given Substituting (12) into (3) yields the following residual: If , then will be the exact solution. Generally such a case will not arise for nonlinear problems, but we can minimize the functional where and are the endpoints of the given problem. The unknown constants can be identified from the conditions With these constants known, the approximate solution (of order ) is well determined.

#### 3. Numerical Examples

To illustrate the effectiveness of the OHAM we will consider three examples of singular two-point BVPs.

Example 1. Consider the following singular two-point BVP [2, 31]: subject to the boundary conditions the exact solution of this problem in case of is given by According to the OHAM formulation described in the above section, we start with Now, apply (4) at to give the zeroth-order problem as follows: with conditions it gives us Now, apply (8) to give the first-order problem as follows: subject to the boundary conditions and having the solution The second-order problem can be defined by (9): subject to the boundary problem and has the solution
By applying (10), the third-order problem is defined as follows: subject to the boundary problem and has the solution Using (24), (27), (30), and (33), the third-order approximate solution by OHAM for is as follows: Following the procedure described in Section 2 on the domain between and , using the residual error, The less square method can be applied as Thus, the following optimal values of ’s are obtained: By considering these values our approximate solution becomes It is clear that the given solution is very close to the exact one since the other terms’ approach zero which leads to that the solution is converge. Moreover, Table 1 exhibits the approximate solution obtained by using the OHAM and CGA [7]. It is clear that the obtained results in a our method are in very good agreement with the exact solution, which proves the reliability of the method. In Figure 1 we plot the approximate solution and the exact solution.

Table 1: Comparison of the exact solution and the OHAM solution for Example 1.
Figure 1: Exact and approximate solution using OHAM for Example 1.

Example 2. Let us consider the singular two-point BVP [2, 3]: The exact solution of this problem in the case of is given by According to the OHAM formulation described in the previous section, we start with
Applying OHAM, we have the following zeroth-, first-, second-, and the third-order problem solutions:
Using (43) we obtain the following third-order approximate solution by OHAM: Following the procedure described in Section 2 on the domain between and , using the residual The following optimal values of ’s are obtained: By considering these values, our approximate solution becomes Therefore, we have the third-order approximate solution of Example 2.

Table 2 shows a comparison between the OHAM solution and the solutions of VIM [3] both with the exact solution. As it is an evident from the compared results, it was found that OHAM gives better results. In Figure 2, both the approximate solution by using OHAM and the exact solution have been plotted.

Table 2: Comparison of the exact solution and the OHAM solution for Example 2.
Figure 2: Exact and approximate solution using OHAM for Example 2.

Example 3. Consider the singular two-point BVP [2, 8] The exact solution of this problem in the case of is given by
According to the OHAM formulation described in the above section, we start with Applying OHAM, we have the following zeroth-, first-, and second-order solutions: Using (52) and also by adding the solutions of the third-order problems, we obtain the following third-order approximate solution by OHAM: Following the procedure described in Section 2 on the domain between and , using the residual The following optimal values of ’s are obtained: The approximate solution now becomes
Table 3 exhibits the approximate solution obtained by using the OHAM and He’s HPM [8]. It can be seen that the solution obtained by our procedure is nearly identical with that given by the exact solution, which proves the reliability of the method. In Figure 3 we compare the exact solution and the approximate solution obtained by OHAM.

Table 3: Comparison of the exact solution and the OHAM solution for Example 3.
Figure 3: Exact and approximate solution using OHAM for Example 3.

#### 4. Conclusions

In this work, OHAM has been applied successfully to solve singular two-point BVPs. The results which are obtained by using OHAM are in a good agreement with the exact solution as well as the results which are already presented in the literature like CGA, VIM, and HPM. This shows that the method is efficient and reliable for the solution of singular two-point boundary value problems.

#### Acknowledgment

The authors thank the anonymous referees for their comments which improved the paper.

#### References

1. S. Chandrasekhar, Hydrodynamic and Hydromagnetic Stability, Dover, New York, NY, USA, 1981.
2. A. S. Bataineh, M. S. M. Noorani, and I. Hashim, “Approximate solutions of singular two-point BVPs by modified homotopy analysis method,” Physics Letters A, vol. 372, no. 22, pp. 4062–4066, 2008.
3. J. Lu, “Variational iteration method for solving two-point boundary value problems,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 92–95, 2007.
4. A. S. V. R. Kanth and Y. N. Reddy, “Cubic spline for a class of singular two-point boundary value problems,” Applied Mathematics and Computation, vol. 170, no. 2, pp. 733–740, 2005.
5. A. S. V. R. Kanth and K. Aruna, “Solution of singular two-point boundary value problems using differential transformation method,” Physics Letters A, vol. 372, no. 26, pp. 4671–4673, 2008.
6. A. Ebaid and M. D. Aljoufi, “Exact solutions for a class of singular two-point boundary value problems using Adomian decomposition method,” Applied Mathematical Sciences, vol. 6, no. 121–124, pp. 6097–6108, 2012.
7. O. A. Arqub, Z. Abo-Hammour, S. Momani, and N. Shawagfeh, “Solving singular two-point boundary value problems using continuous genetic algorithm,” Abstract and Applied Analysis, vol. 2012, Article ID 205391, 25 pages, 2012.
8. K. Al-Khaled, “Theory and computation in singular boundary value problems,” Chaos, Solitons and Fractals, vol. 33, no. 2, pp. 678–684, 2007.
9. M. Jalaal, D. D. Ganji, and G. Ahmadi, “Analytical investigation on acceleration motion of a vertically falling spherical particle in incompressible Newtonian media,” Advanced Powder Technology, vol. 21, no. 3, pp. 298–304, 2010.
10. M. Jalaal and D. D. Ganji, “An analytical study on motion of a sphere rolling down an inclined plane submerged in a Newtonian fluid,” Powder Technology, vol. 198, no. 1, pp. 82–92, 2010.
11. M. Jalaal and D. D. Ganji, “On unsteady rolling motion of spheres in inclined tubes filled with incompressible Newtonian fluids,” Advanced Powder Technology, vol. 22, no. 1, pp. 58–67, 2011.
12. M. Jalaal, M. G. Nejad, P. Jalili et al., “Homotopy perturbation method for motion of a spherical solid particle in plane couette fluid flow,” Computers and Mathematics with Applications, vol. 61, no. 8, pp. 2267–2270, 2011.
13. M. Esmaeilpour and D. D. Ganji, “Application of He's homotopy perturbation method to boundary layer flow and convection heat transfer over a flat plate,” Physics Letters A, vol. 372, no. 1, pp. 33–38, 2007.
14. S. J. Liao, The proposed homotopy analysis techniques for the solution of nonlinear problems [Ph.D. thesis], Shanghai Jiao Tong University, 1992.
15. S. J. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Method, vol. 2 of CRC Series: Modern Mechanics and Mathematics, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2004.
16. S. J. Liao, “An explicit, totally analytic approximate solution for Blasius' viscous flow problems,” International Journal of Non-Linear Mechanics, vol. 34, no. 4, pp. 759–778, 1999.
17. S. J. Liao, “On the homotopy analysis method for nonlinear problems,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 499–513, 2004.
18. S. J. Liao, “An approximate solution technique not depending on small parameters: a special example,” International Journal of Non-Linear Mechanics, vol. 30, no. 3, pp. 371–380, 1995.
19. S. J. Liao, “A kind of approximate solution technique which does not depend upon small parameters—II: an application in fluid mechanics,” International Journal of Non-Linear Mechanics, vol. 32, no. 5, pp. 815–822, 1997.
20. S. J. Liao, “Comparison between the homotopy analysis method and homotopy perturbation method,” Applied Mathematics and Computation, vol. 169, no. 2, pp. 1186–1194, 2005.
21. S. J. Liao, “A new branch of solutions of boundary-layer flows over an impermeable stretched plate,” International Journal of Heat and Mass Transfer, vol. 48, no. 12, pp. 2529–2539, 2005.
22. T. Hayat and M. Khan, “Homotopy solutions for a generalized second-grade fluid past a porous plate,” Nonlinear Dynamics, vol. 42, no. 4, pp. 395–405, 2005.
23. M. Ayub, A. Rasheed, and T. Hayat, “Exact flow of a third grade fluid past a porous plate using homotopy analysis method,” International Journal of Engineering Science, vol. 41, no. 18, pp. 2091–2103, 2003.
24. I. T. Abu-Zaid and M. A. El-Gebeily, “A finite difference method for approximating the solution of a certain class of singular two-point boundary value problems,” Arab Journal of Mathematical Sciences, vol. 1, no. 1, pp. 25–39, 1995.
25. A. M. Siddiqui, M. Ahmed, and Q. K. Ghori, “Couette and poiseuille flows for non-newtonian fluids,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 1, pp. 15–26, 2006.
26. N. Herisanu, V. Marinca, T. Dordea, and G. Madescu, “A new analytical approach to nonlinear vibration of an electric machine,” Proceedings of Romanian Academy A, vol. 9, no. 3, 2008.
27. V. Marinca, N. Herisanu, 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.
28. V. Marinca, N. Herisanu, and I. Nemes, “Optimal homotopy asymptotic method with application to thin film flow,” Central European Journal of Physics, vol. 6, no. 3, pp. 648–653, 2008.
29. V. Marinca and N. Herisanu, “Application of optimal homotopy asymptotic method for solving nonlinear equations arising in heat transfer,” International Communications in Heat and Mass Transfer, vol. 35, no. 6, pp. 710–715, 2008.
30. J. Ali, S. Islam, S. Islam, and G. Zaman, “The solution of multipoint boundary value problems by the optimal homotopy asymptotic method,” Computers & Mathematics with Applications, vol. 59, no. 6, pp. 2000–2006, 2010.
31. M. Cui and F. Geng, “Solving singular two-point boundary value problem in reproducing kernel space,” Journal of Computational and Applied Mathematics, vol. 205, no. 1, pp. 6–15, 2007.