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.

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.

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.

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.