Research Article  Open Access
Variational Iteration Method for Nonlinear Singular TwoPoint Boundary Value Problems Arising in Human Physiology
Abstract
The variational iteration method is applied to solve a class of nonlinear singular boundary value problems that arise in physiology. The process of the method, which produces solutions in terms of convergent series, is explained. The Lagrange multipliers needed to construct the correctional functional are found in terms of the exponential integral and Whittaker functions. The method easily overcomes the obstacle of singularities. Examples will be presented to test the method and compare it to other existing methods in order to confirm fast convergence and significant accuracy.
1. Introduction
In this paper, He’s variational iteration method (VIM) [1, 2] is applied to obtain an approximate solution for the following nonlinear singular twopoint boundary value problem (BVP): with the following two sets of boundary conditions: where Here , , and and are finite constants. It is assumed that is nonnegative, continuously differentiable on , and is analytic in the disk for some . It is also assumed that and are both continuous on and further .
Assuming that is analytic in for some , the existenceuniqueness has been established for problem (1) under the following restrictions (see [3] and the references therein):(1)BC (2) holds for , and such that , and(2)BC (3) holds for , , and such that .
Equation (1) with BC (3) arises in the study of tumor growth problems [4, 5] where , , and is of the form This problem also arises in the study of a steadystate oxygen diffusion in a cell with MichaleisMenten uptake kinetics when and [6]. Another application of this problem appeared in the study of heat sources distribution in the human head [7] given that , , and
Problems (1)(2) and (1)–(3) have been treated by several authors using different numerical schemes such as cubic Bsplines, finite difference, and pointwise solution bounds [3, 8, 9]. The VIM has been used by Wazwaz [10] to solve the nonlinear singular EmdenFowler boundary value problem, which is a special case of (1)(2), where . Also, Ravi Kanth and Aruna [11] used the variational iteration method to solve (1)(2) and (1)–(3), but that was also done under the assumption that .
The objective of this paper is to apply the VIM to obtain an approximate solution for the proposed nonlinear singular boundary value problems with more relaxed . The VIM has been intensively used to solve linear and nonlinear, ordinary and partial, delay and fractional order differential equations [12–18]. In a recent review article, He [19] shows how the VIM can be employed as an effective method in searching for wave solutions including solitons and compacton solutions without the need for linearization or weak nonlinearity assumptions. For nonlinear differential equations, the VIM gives an approximate solution in a series form that converges to the exact solution, if the latter exists. In case of linear differential equations, exact solution by the VIM can be readily obtained by a one iteration step because the exact Lagrange multiplier needed to construct the correctional functional can be identified.
2. Variational Iteration Method
Consider the general nonlinear differential equation given by where and are linear and nonlinear operators, respectively, and is a given analytical function. As per the variational iteration method, the correction functional is constructed as follows: where is a general Lagrange multiplier, which can be optimally identified by the variational theory, and is a restricted variation and hence . A proper choice of the initial approximation is crucial to obtain a successive approximations that converge as rapidly as possible to the exact solution.
For simplicity, (1) will be expressed in the form and its correction functional is given by
Taking the variation with respect to gives
Using integration by parts in (11) gives which in turn leads to the stationary relations
Depending on , the Lagrange multiplier can be determined by solving system (13). In the next section, we will use , and hence the choice of will lead to one of the following two Lagrange multipliers.
If , the solution of (13) is given in terms of the exponential integral [20]. For example, if , then where is the generalized complex exponential integral given by
If , the solution of (13) is given in terms of Whittaker function [20]. For example, if , then where and is given by in which is known as the Pochhammer symbol [20] and is defined by
Now that is constructed, an appropriate first guess will lead to the formation of the recurrence sequences . The solution for (1)(2) and (1)–(3) is, then, obtained from the relation provided that the limit exists. The proof that the sequences are convergent has been well established in the literature (see, e.g., [11, 21]).
3. Numerical Examples
Example 1. Consider the nonlinear singular BVP where
The analytic solution of (21) is for all . This problem is an application of oxygen diffusion corresponding to (1), (2), and (5) with , , , , , and . In Tables 1 and 3, the second and third iteration solutions are compared to the exact solution when and , respectively. In Tables 2 and 4, the absolute maximum errors obtained by the proposed VIM is compared with that of other methods.


Example 2. Consider the nonlinear singular BVP where , , and are as in (22).
This problem is an application of the nonlinear heat conduction model of the human head corresponding to (1), (3), and (6) with , , and . In Tables 5 and 7, the second and third iteration solutions are compared to the exact solution when and , respectively, and in Tables 6 and 8, the absolute maximum errors obtained by the proposed VIM is compared with that of other methods.


4. Conclusion
In this paper, the simplicity and reliability of He’s variational iteration method have been tested against a general form of nonlinear two point singular value problems that arise in physiology. The method is shown to be highly accurate and easily implemented even with the presence of singularities. For the two examples shown in Section 3, the second and third iteration solutions obtained by the VIM were substantially more accurate than the 16 and the 32 timestep approximate solutions obtained by cubic Bspline and the finite difference methods.
References
 J. H. He, “Variational iteration method—a kind of nonlinear analytical technique: some examples,” International Journal of NonLinear Mechanics, vol. 34, no. 4, pp. 699–708, 1999. View at: Publisher Site  Google Scholar
 J.H. He, “Variational iteration method—some recent results and new interpretations,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 3–17, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. K. Pandey and A. K. Singh, “On the convergence of a finite difference method for a class of singular boundary value problems arising in physiology,” Journal of Computational and Applied Mathematics, vol. 166, no. 2, pp. 553–564, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. A. Adam, “A simplified mathematical model of tumor growth,” Mathematical Biosciences, vol. 81, no. 2, pp. 229–244, 1986. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 H. P. Greenspan, “Models for the growth of a solid tumor by diffusion,” Studies in Applied Mathematics, vol. 4, pp. 317–340, 1972. View at: Google Scholar  Zentralblatt MATH
 S. H. Lin, “Oxygen diffusion in a spherical cell with nonlinear oxygen uptake kinetics,” Journal of Theoretical Biology, vol. 60, no. 2, pp. 449–457, 1976. View at: Publisher Site  Google Scholar
 B. F. Gray, “The distribution of heat sources in the human head—theoretical consideration,” Journal of Theoretical Biology, vol. 82, no. 3, pp. 473–476, 1980. View at: Publisher Site  Google Scholar
 N. S. Asaithambi and J. B. Garner, “Pointwise solution bounds for a class of singular diffusion problems in physiology,” Applied Mathematics and Computation, vol. 30, no. 3, pp. 215–222, 1989. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Abukhaled, S. A. Khuri, and A. Sayfy, “A numerical approach for solving a class of singular boundary value problems arising in physiology,” International Journal of Numerical Analysis and Modeling, vol. 8, no. 2, pp. 353–363, 2011. View at: Google Scholar  Zentralblatt MATH
 A. M. Wazwaz, “The variational iteration method for solving nonlinear singular boundary value problems arising in various physical models,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 10, pp. 3881–3886, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. S. V. Ravi Kanth and K. Aruna, “He's variational iteration method for treating nonlinear singular boundary value problems,” Computers & Mathematics with Applications, vol. 60, no. 3, pp. 821–829, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. He, “Variational iteration method for delay differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 2, no. 4, pp. 235–236, 1997. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J.H. He, “Variational iteration method for autonomous ordinary differential systems,” Applied Mathematics and Computation, vol. 114, no. 23, pp. 115–123, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J.H. He, “Variational principles for some nonlinear partial differential equations with variable coefficients,” Chaos, Solitons & Fractals, vol. 19, no. 4, pp. 847–851, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J.H. He, “Some asymptotic methods for strongly nonlinear equations,” International Journal of Modern Physics B, vol. 20, no. 10, pp. 1141–1199, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. M. Odibat and S. Momani, “Application of variational iteration method to nonlinear differential equations of fractional order,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 1, pp. 27–34, 2006. View at: Google Scholar
 S. Momani and Z. Odibat, “Numerical comparison of methods for solving linear differential equations of fractional order,” Chaos, Solitons & Fractals, vol. 31, no. 5, pp. 1248–1255, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Lu, “Variational iteration method for solving twopoint boundary value problems,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 92–95, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. H. He, “Asymptotic methods for solitary solutions and compactons,” Abstract and Applied Analysis, vol. 2012, Article ID 916793, 130 pages, 2012. View at: Google Scholar
 M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover, New York, NY, USA, 1972.
 D. Khojasteh Salkuyeh, “Convergence of the variational iteration method for solving linear systems of ODEs with constant coefficients,” Computers & Mathematics with Applications, vol. 56, no. 8, pp. 2027–2033, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
Copyright
Copyright © 2013 Marwan Abukhaled. 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.