Abstract
In this paper, optimal homotopy asymptotic method (OHAM) for the semianalytic solutions of nonlinear singular two-point boundary value problems has been applied to several problems. The solutions obtained by OHAM have been compared with the solutions of another method named as modified adomain decomposition (MADM). For testing the success of OHAM, both of the techniques have been analyzed against the exact solutions in all problems. It is proved by this paper that solutions of OHAM converge rapidly to the exact solution and show most effectiveness as compared to MADM.
1. Introduction
Nonlinear singular boundary values problem (SBVP) is of considerable importance in the modeling of many branches of mathematical physics and engineering i-e fluid mechanics, quantum mechanics, astrophysics and so forth. Nonlinear problems are difficult to solve as compared to linear problems but the singularity of such nonlinear problem becomes more difficult to solve by researcher. Such type of problems have been successfully solved by modified adomian decomposition method [1, 2], but the convergence rate of its series goes gradually to accurate solution. Recently a new growing method has been introduced by Marinca et al. [3–6], for approximate solution of the problems of thin film flow of a fourth grade fluid on a cylinder surface and for understanding the behavior of nonlinear mechanical vibration of an electrical machines. OHAM attracted great attention of many researchers to solve a large class of linear and nonlinear differential equations and so forth; Idrees et al. extended OHAM by applying on initial and boundary values problems of ODEs and PDEs [7–10]. Iqbal and Javed [11] applied OHAM for the analytical solution of singular Lane-Emden type equations. The goal of this paper is to apply the proposed method to compute best approximate solutions of nonlinear singular boundary values problems. This paper is arranged in four sections; the first section of the paper describes the introduction, the second section of the paper is about the basic idea of OHAM, and the third section presents numerical results to demonstrate the efficiency of OHAM, and the last section of the paper is the conclusion.
2. Basic Mathematical Theory of Optimal Homotopy Asymptotic Method
The framework of OHAM is given in few steps.
(i) Consider the undertaken differential equation where and are independent variables and and are linear and nonlinear part of (1), respectively. Simply linear operator is denoted by and nonlinear by , is taken as an unknown function, is considered as known function, and yield corresponding initial condition. Now
(ii) According to OHAM, an optimal homotopy, can be built, fulfilling homotopy equation: is auxiliary function with , an embedding parameter, such that varies from 0 to 1 and the solution varies from to the final solution , where is evaluated from (2) for : is reserved in the series as below:
(iii) Next can be spread out by Taylor’s series about the parameter : , cause the accurateness and convergence of (6) and if it converges at then we have the following (see [5]): Insert (6) into homotopy formula (2) and comparing like powers of , the original nonlinear problem transferred into a sequence of linear problems, we obtained zeroth order (4), first order (8), second order (9), and th order (10) as below:
(iv) The resultant linear problems can now be solved and their solutions are used to construct th order solution that comprises of the original problem through (7). Then by substituting (7) into (1), we get the following residual: when , for some values of then will agree with the exact solution. Though, this does not occur generally, especially, in nonlinear problems. Therefore, optimal values of the auxiliary constants are calculated for minimizing the following functional (see [5]): Therefore, the unknown constants () can be optimally identified from the following conditions (see [5]): With these known values of the auxiliary constants the approximate solution (7) is now well determined.
3. Application of OHAM
OHAM is implemented on three models of singular boundary value problems in order to show the accuracy of this method.
Model 3.1. Consider second order singular boundary value problem [1]: subject to boundary conditions with exact solution
Zeroth Order Problem. Consider the following: Its solution is
First Order Problem. Consider the following: Its solution is
Second Order Problem. Consider the following: Its solution is Solve (17) to (22) to obtain approximate solution in the form of To determine constant values, and , apply the method of least square (12), (13) where By substituting values of and in (22), we achieved the approximate solution of OHAM.
Model 3.2. Consider the following third order nonlinear boundary value problem [2]: with boundary conditions where
Zeroth Order Problem. Consider the following: Its solution is
First Order Problem. Consider the following: Its solution is
Second Order Problem. Consider the following: Its solution is Solving problem (28) to (33), in the form of We obtained the constant values of and as Taylor series of the exact solution with order 10, which is given below
Model 3.3. Consider the following third order singular boundary value problem [1]: subject to boundary conditions where Exact solution of (41) is
Zeroth Order Problem. Consider the following: Its solution is
First Order Problem. Consider the following: Its solution is
Second Order Problem. Consider the following: Its solution is Solving problem (41) to (47), we obtained We obtained the constant values of and as Substituting constant values of and in (48) to obtain approximate solution of OHAM.
4. Conclusion
Optimal homotopy asymptotic method has been applied to obtain the approximate solution of singular boundary value problems. Results have been compared with modified adomian decomposition method and with the exact solutions of proposed models. Numerical results present in Tables 1, 2, and 3 showed that this technique is fast convergent and has remarkable low error than MADM. It is observed through represented plots in Figures 1(a), 2(a), and 3(a) that the presented technique is effective and reliable within bounded domain. Absolute errors plots by Figures 1(b), 2(b), and 3(b) showed the precise closeness of OHAM to exact results.
(a)
(b)
(a)
(b)
(a)
(b)
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.