Mathematical Problems in Engineering

Volume 2012 (2012), Article ID 208375, 14 pages

http://dx.doi.org/10.1155/2012/208375

## A General Solution for Troesch's Problem

^{1}Electronic Instrumentation and Atmospheric Sciences School, University of Veracruz, Xalapa, VER 91000, Mexico^{2}Department of Mathematics, Zhejiang University, Hangzhou 310027, China^{3}Departamento de Física y Matemáticas, Universidad Iberoamericana, Prolongación Paseo de la Reforma 880, 01219 Mexico DF, Mexico^{4}Micro and Nanotechnology Research Center, University of Veracruz, Calzada Ruiz Cortines 455, 94292 Boca del Rio, VER, Mexico^{5}National Institute for Astrophysics, Optics and Electronics Luis Enrique Erro No.1, 72840 Santa María Tonantzintla, PUE, Mexico

Received 23 August 2012; Accepted 6 October 2012

Academic Editor: Farzad Khani

Copyright © 2012 Hector Vazquez-Leal 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

The homotopy perturbation method (HPM) is employed to obtain an approximate solution for the nonlinear differential equation which describes Troesch’s problem. In contrast to other reported solutions obtained by using variational iteration method, decomposition method approximation, homotopy analysis method, Laplace transform decomposition method, and HPM method, the proposed solution shows the highest degree of accuracy in the results for a remarkable wide range of values of Troesch’s parameter.

#### 1. Introduction

Troesch’s equation is a boundary value problem (BVP) expressed as where prime denotes differentiation with respect to and is known as Troesch’s parameter.

Equation (1.1) arises in the investigation of confinement of a plasma column by a radiation pressure [1] and also in the theory of gas porous electrodes [2, 3]. This BVP problem has a pole [4] approximately located at which makes the solution of (1.1) a difficult task for numerical methods.

In order to overcome such difficulties, there are several reported numerical solutions for Troesch’s problem [4–10]. Recently, after some decades using wrong numerical results, in [8, 9] were reported the most accurate solutions for and . Besides, there are approximated analytical solutions obtained by using different methods like homotopy perturbation method (HPM) [11, 12], decomposition method approximation (DMA) [12, 13], homotopy analysis method (HAM) [14], variational iteration method (VIM) [15], and Laplace transform decomposition method (LTDM) [16]. The main disadvantage of aforementioned approximated solutions is that they are obtained for specific values of , like , or . In contrast, we propose a general approximate solution for Troesch’s problem, useful for , by using HPM method [17–27].

This paper is organized as follows. In Section 2, we provide a brief review of HPM method. In Section 3, we obtain the solution of Troesch’s problem employing HPM. Section 4 shows numerical simulations and discuss our findings. Finally, a brief conclusion is given in Section 5.

#### 2. Basic Idea of HPM Method

Basically, the HPM method [17–26, 31–53] introduces a homotopy parameter , which takes values ranging from up to . When parameter , the equation usually reduces to a simple, or trivial, equation to solve. Then, is gradually increased to 1, producing a sequence of deformations. Eventually, at , the homotopy equation takes the original form of the equation to solve, and the final stage of deformation provides the desired solution. Usually, few iterations are required to obtain good results [17–19].

In the HPM method, it is considered that a nonlinear differential equation can be expressed as with the boundary condition where is a general differential operator, is a known analytic function, is a linear initial/boundary operator, is the boundary of domain , and denotes differentiation along the normal drawn outwards from . The operator, generally, can be divided into two operators, and , where is the linear operator, and is the nonlinear operator. Hence, (2.1) can be rewritten as

Now, a possible homotopy formulation is where is the initial approximation for the solution of (2.3) which satisfies the boundary conditions, and is known as the perturbation homotopy parameter.

We assume that the solution of (2.4) can be written as a power series of as

When , results that the approximate solution for (2.1) is

The series (2.6) is convergent on most cases, nevertheless, the convergence depends on the nonlinear operator [19, 21, 35, 36, 52].

#### 3. Solution of Troesch’s Problem by Using HPM Method

Straight forward application of HPM to solve (1.1) is not possible due to the hyperbolic term of dependent variable. However, in [11, 12] were reported HPM solutions, obtained by using a power series expansion of the term of (1.1), which are limited to , due to the truncate power series. Nevertheless, the polynomial type nonlinearities are easier to handle by the HPM method. Therefore, in order to apply HPM successfully for a wide range of Troesch’s parameter (), we convert the hyperbolic-type nonlinearity in Troesch’s problem into a polynomial type nonlinearity, using the variable transformation reported in [15]. First, we consider that from which we find where prime denotes differentiation with respect to .

Now, by substituting (3.1) into the term of (1.1), we obtain

Then, equating (3.3) and (3.4), we achieve to the following transformed problem: where conditions are obtained by using variable transformation (see (3.1)) and substituting original boundary conditions and into above equation, results

From (2.4) and (3.5), we can formulate the following homotopy [17–19]: where is the homotopy parameter.

Substituting (2.5) into (3.8) and equating identical powers of terms, we obtain

We solve (3.9) by using Maple software, resulting where is and is

Next, calculating the limit when , we obtain the second-order approximated solution of (3.5)

Finally, from (3.1) and (3.13), the proposed solution of Troesch’s problem is

##### 3.1. Interval of Solution

The real branch of is restricted to the range . Therefore, (3.14) requires where is delimited by the boundary conditions as

In order to show that (3.14) fullfills the conditions (3.15) and (3.16), we plot (3.13) in Figure 1. From such figure, we can observe that is valid in the intervals of and . Now, from (3.13) and (3.14), we calculate the following limits: From (3.17), (3.18), (3.19), and Figure 1, we can conclude that the maximum value of (3.13) is 1 in the range of . Therefore, (3.13) fullfills (3.15) in the range given by (3.16). Additionally, limit (3.19) shows that for , the presented solution (3.14) becomes the exact/trivial solution for (1.1).

#### 4. Numerical Simulation and Discussion

In the case of (see Table 1), we can observe that the lowest average absolute relative error (A.A.R.E.) is for LDTM [16], followed closely by the proposed solution (3.14). A possible reason can lie in the fact that (3.14) is a second-order approximation, while LDTM is of third order. For (see Table 2), there is a change now the lowest A.A.R.E. is for the proposed solution (3.14), followed by LDTM solution. For both cases, the other approximations ADM [13], HPM [12], HPM [11], and HAM [14] have lower accuracy than (3.14). Equation (3.14) did not require an adjustment parameter, unlike LDTM solution, which required a specific adjustment parameter calculated for each value of Troesch’s parameter . Therefore, the proposed solution is easier to use than LDTM solution.

In Table 3, we can observe a comparison of (3.14) with numerical solution [15], and other solutions obtained by HAM [14], VIM [15], and HPM [12] for . Approximation (3.14) has the lowest A.A.R.E. from all above solutions, followed by VIM approximation. Besides, HAM [14] has an relatively poor value of A.A.R.E., despite the fact that it is a sixth-order approximation. Furthermore, HPM [12] shows divergence from the numerical solution. In addition, (3.14) do not require an adjustment parameter, nevertheless, VIM solution required a specific adjustment parameter calculated for each value of Troesch’s parameter . Therefore, the proposed solution is easier to use than VIM solution.

In Table 4, is presented a comparison of initial slope and the results reported in [4, 5, 15, 28, 29] for the range of ; resulting that the proposed solution is the only one reported in literature with high accuracy in the complete aforementioned range. Moreover, in order to compare the derivative of (3.14) for , we use the approximate reported in [29] resulting in a remarkable accuracy for the range .

In the same fashion, in Table 5, we compare for proposed solution and other reported numerical solutions [4, 5, 28, 30]. The results shows that (3.14) has a good accuracy for at least in the range .

Figure 2 results from Tables 1, 2, and 3 and Maple routines of numerical solution of differential equations. We can observe the overlap of numerical solution and (3.14), for different values of . Therefore, (3.14) is a high accurate approximate solution, valid for a wide range of values of Troesch’s parameter.

Further research should be performed in order to verify the accuracy of (3.14) for large values of . Nevertheless, as we know from above discussion and are valid and accurate in the ranges and , respectively, therefore, we can expect a high accuracy of (3.14) in the same ranges.

#### 5. Concluding Remarks

In this work, we obtained an approximate solution for Troesch’s problem. Besides, we presented a comparison between the numerical solution, the proposed solution, and other approximations reported in the literature. The numerical and graphical results show that the proposed solution is the most accurate one, for a wide range of values of Troesch’s parameter (). Moreover, Troesch’s problem approximation exhibit a remarkable accuracy for at least in the range . In the same fashion, exhibits a good accuracy at least in the range . Accordingly, we can expect a good accuracy of proposed solution for the range and maybe for even more larger values of . Additionally, the proposed approximated solution does not require any adjustment parameter as reported for solutions obtained by using LDTM and VIM methods, which makes our proposed solution easier to use than those approximations. Finally, further research is necessary in order to verify the accuracy of our proposed approximation for large values of .

#### Acknowledgments

The authors gratefully acknowledge the financial support of the National Council for Science and Technology of Mexico (CONACYT) through Grant CB-2010-01 no. 157024. The authors would like to thank Rogelio-Alejandro Callejas-Molina and Roberto Ruiz-Gomez for their contribution to this project.

#### References

- E. S. Weibel, “On the confinement of a plasma by magnetostatic fields,”
*Physics of Fluids*, vol. 2, no. 1, pp. 52–56, 1959. View at Google Scholar · View at Scopus - D. Gidaspow and B. S. Baker, “A model for discharge of storage batteries,”
*Journal of the Electrochemical Society*, vol. 120, no. 8, pp. 1005–1010, 1973. View at Google Scholar · View at Scopus - V. S. Markin, A. A. Chernenko, Y. A. Chizmadehev, and Y. G. Chirkov, “Aspects of the theory of gas porous electrodes,” in
*Fuel Cells: Their Electrochemical Kinetics*, pp. 22–33, Consultants Bureau, New York, NY, USA, 1966. View at Google Scholar - S. M. Roberts and J. S. Shipman, “On the closed form solution of Troesch's problem,”
*Journal of Computational Physics*, vol. 21, no. 3, pp. 291–304, 1976. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - S.-H. Chang, “Numerical solution of Troesch's problem by simple shooting method,”
*Applied Mathematics and Computation*, vol. 216, no. 11, pp. 3303–3306, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - S. A. Khuri and A. Sayfy, “Troesch's problem: a b-spline collocation approach,”
*Mathematical and Computer Modelling*, vol. 54, no. 9-10, pp. 1907–1918, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. Zarebnia and M. Sajjadian, “The sinc-galerkin method for solving troesch's problem,”
*Mathematical and Computer Modelling*, vol. 56, no. 9-10, pp. 218–228, 2012. View at Publisher · View at Google Scholar · View at Scopus - U. Erdogan and T. Ozis, “A smart nonstandard finite difference scheme for second order nonlinear boundary value problems,”
*Journal of Computational Physics*, vol. 230, no. 17, pp. 6464–6474, 2011. View at Publisher · View at Google Scholar - Y. Lin, J. A. Enszer, and M. A. Stadtherr, “Enclosing all solutions of two-point boundary value problems for ODEs,”
*Computers and Chemical Engineering*, vol. 32, no. 8, pp. 1714–1725, 2008. View at Publisher · View at Google Scholar · View at Scopus - D. J. Jones, “Solution of Troesch's, and other, two point boundary value problems by shooting techniques,”
*Journal of Computational Physics*, vol. 12, no. 3, pp. 429–434, 1973. View at Google Scholar · View at Scopus - S. H. Mirmoradi, I. Hosseinpour, S. Ghanbarpour, and A. Barari, “Application of an approximate analytical method to nonlinear Troesch's problem,”
*Applied Mathematical Sciences*, vol. 3, no. 29-32, pp. 1579–1585, 2009. View at Google Scholar - X. Feng, L. Mei, and G. He, “An efficient algorithm for solving Troesch's problem,”
*Applied Mathematics and Computation*, vol. 189, no. 1, pp. 500–507, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - E. Deeba, S. A. Khuri, and S. Xie, “An algorithm for solving boundary value problems,”
*Journal of Computational Physics*, vol. 159, no. 2, pp. 125–138, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - H. N. Hassan and M. A. El-Tawil, “An efficient analytic approach for solving two-point nonlinear boundary value problems by homotopy analysis method,”
*Mathematical Methods in the Applied Sciences*, vol. 34, no. 8, pp. 977–989, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - S.-H. Chang, “A variational iteration method for solving Troesch's problem,”
*Journal of Computational and Applied Mathematics*, vol. 234, no. 10, pp. 3043–3047, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - S. A. Khuri, “A numerical algorithm for solving Troesch's problem,”
*International Journal of Computer Mathematics*, vol. 80, no. 4, pp. 493–498, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J.-H. He, “Comparison of homotopy perturbation method and homotopy analysis method,”
*Applied Mathematics and Computation*, vol. 156, no. 2, pp. 527–539, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J.-H. He, “An elementary introduction to the homotopy perturbation method,”
*Computers & Mathematics with Applications*, vol. 57, no. 3, pp. 410–412, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J.-H. He, “Homotopy perturbation technique,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 178, no. 3-4, pp. 257–262, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - H. Vázquez-Leal, R. Castaneda-Sheissa, U. Filobello-Nino, A. Sarmiento-Reyes, and J. Sanchez Orea, “High accurate simple approximation of normal distribution integral,”
*Mathematical Problems in Engineering*, vol. 2012, Article ID 124029, 22 pages, 2012. View at Publisher · View at Google Scholar · View at Scopus - H. Vázquez-Leal, U. Filobello-Niño, R. Castañeda-Sheissa, L. Hernández-Martínez, and A. Sarmiento-Reyes, “Modified HPMs inspired by homotopy continuation methods,”
*Mathematical Problems in Engineering*, vol. 2012, Article ID 309123, 19 pages, 2012. View at Publisher · View at Google Scholar · View at Scopus - U. Filobello-Nino, H. Vázquez-Leal, R. Castaneda-Sheissa et al., “An approximate solution of blasius equation by using hpm method,”
*Asian Journal of Mathematics and Statistics*, vol. 5, pp. 50–59, 2012. View at Google Scholar - Y. Khan, Q. Wu, N. Faraz, and A. Yildirim, “The effects of variable viscosity and thermal conductivity on a thin film flow over a shrinking/stretching sheet,”
*Computers & Mathematics with Applications*, vol. 61, no. 11, pp. 3391–3399, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - N. Faraz and Y. Khan, “Analytical solution of electrically conducted rotating flow of a second grade fluid over a shrinking surface,”
*Ain Shams Engineering Journal*, vol. 2, no. 3-4, pp. 221–226, 2011. View at Publisher · View at Google Scholar · View at Scopus - Y. Khan, Q. Wu, N. Faraz, A. Yildirim, and M. Madani, “A new fractional analytical approach via a modified Riemann-Liouville derivative,”
*Applied Mathematics Letters*, vol. 25, no. 10, pp. 1340–1346, 2012. View at Publisher · View at Google Scholar · View at Scopus - Y. Khan, H. Vázquez-Leal, and Q. Wu, “An efficient iterated method for mathematical biology model,”
*Neural Computing and Applications*. In press. View at Publisher · View at Google Scholar - U. Filobello-Nino, H. Vazquez-Leal, Y. Khan et al., “HPM applied to solve nonlinear circuits: a study case,”
*Applied Mathematical Sciences*, vol. 6, no. 85-88, pp. 4331–4344, 2012. View at Google Scholar - B. A. Troesch, “A simple approach to a sensitive two-point boundary value problem,”
*Journal of Computational Physics*, vol. 21, no. 3, pp. 279–290, 1976. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. Kubiček and V. Hlaváček, “Solution of Troesch's two-point boundary value problem by shooting technique,”
*Journal of Computational Physics*, vol. 17, no. 1, pp. 95–101, 1975. View at Google Scholar · View at Scopus - J. A. Snyman, “Continuous and discontinuous numerical solutions to the Troesch problem,”
*Journal of Computational and Applied Mathematics*, vol. 5, no. 3, pp. 171–175, 1979. View at Google Scholar · View at Scopus - H. Koçak, A. Yıldırım, D. H. Zhang, and S. T. Mohyud-Din, “The comparative Boubaker polynomials expansion scheme (BPES) and homotopy perturbation method (HPM) for solving a standard nonlinear second-order boundary value problem,”
*Mathematical and Computer Modelling*, vol. 54, no. 1-2, pp. 417–422, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. Yildirim and H. Koçak, “Homotopy perturbation method for solving the space-time fractional advection-dispersion equation,”
*Advances in Water Resources*, vol. 32, no. 12, pp. 1711–1716, 2009. View at Publisher · View at Google Scholar · View at Scopus - A. Yildirim and H. Koçak, “Series solution of the Smoluchowski's coagulation equation,”
*Journal of King Saud University. Science*, vol. 23, no. 23, pp. 183–189, 2010. View at Publisher · View at Google Scholar · View at Scopus - H. Koçak and A. Yildirim, “Numerical solution of 3D Green's function for the dynamic system of anisotropic elasticity,”
*Physics Letters, Section A*, vol. 373, no. 35, pp. 3145–3150, 2009. View at Publisher · View at Google Scholar · View at Scopus - J. Biazar and H. Aminikhah, “Study of convergence of homotopy perturbation method for systems of partial differential equations,”
*Computers & Mathematics with Applications*, vol. 58, no. 11-12, pp. 2221–2230, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J. Biazar and H. Ghazvini, “Convergence of the homotopy perturbation method for partial differential equations,”
*Nonlinear Analysis. Real World Applications*, vol. 10, no. 5, pp. 2633–2640, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J. Biazar, F. Badpeima, and F. Azimi, “Application of the homotopy perturbation method to Zakharov-Kuznetsov equations,”
*Computers & Mathematics with Applications*, vol. 58, no. 11-12, pp. 2391–2394, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J. Biazar and H. Ghazvini, “Exact solutions for non-linear Schrödinger equations by He's homotopy perturbation method,”
*Physics Letters. A*, vol. 366, no. 1-2, pp. 79–84, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J. Biazar, H. Ghazvini, and M. Eslami, “He's homotopy perturbation method for systems of integro-differential equations,”
*Chaos, Solitons and Fractals*, vol. 39, no. 3, pp. 1253–1258, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J. Biazar, M. Eslami, and H. Aminikhah, “Application of homotopy perturbation method for systems of Volterra integral equations of the first kind,”
*Chaos, Solitons and Fractals*, vol. 42, no. 5, pp. 3020–3026, 2009. View at Publisher · View at Google Scholar · View at Scopus - J. Biazar and H. Ghazvini, “He's homotopy perturbation method for solving systems of Volterra integral equations of the second kind,”
*Chaos, Solitons and Fractals*, vol. 39, no. 2, pp. 770–777, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J. Biazar and H. Ghazvini, “Numerical solution for special non-linear Fredholm integral equation by HPM,”
*Applied Mathematics and Computation*, vol. 195, no. 2, pp. 681–687, 2008. View at Publisher · View at Google Scholar - M. Madani, M. Fathizadeh, Y. Khan, and A. Yildirim, “On the coupling of the homotopy perturbation method and Laplace transformation,”
*Mathematical and Computer Modelling*, vol. 53, no. 9-10, pp. 1937–1945, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - N. Faraz, Y. Khan, and A. Yildirim, “Analytical approach to two-dimensional viscous flow with a shrinking sheet via variational iteration algorithm-II,”
*Journal of King Saud University. Science*, vol. 23, no. 1, pp. 77–81, 2011. View at Publisher · View at Google Scholar · View at Scopus - D. D. Ganji and A. Rajabi, “Assessment of homotopy-perturbation and perturbation methods in heat radiation equations,”
*International Communications in Heat and Mass Transfer*, vol. 33, no. 3, pp. 391–400, 2006. View at Publisher · View at Google Scholar · View at Scopus - M. Rafei, D. D. Ganji, and H. Daniali, “Solution of the epidemic model by homotopy perturbation method,”
*Applied Mathematics and Computation*, vol. 187, no. 2, pp. 1056–1062, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - D. D. Ganji, “The application of He's homotopy perturbation method to nonlinear equations arising in heat transfer,”
*Physics Letters. A*, vol. 355, no. 4-5, pp. 337–341, 2006. View at Publisher · View at Google Scholar - D. D. Ganji and A. Sadighi, “Application of homotopy-perturbation and variational iteration methods to nonlinear heat transfer and porous media equations,”
*Journal of Computational and Applied Mathematics*, vol. 207, no. 1, pp. 24–34, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. Sheikholeslami, H. R. Ashorynejad, D. D. Ganji, and A. Yildirim, “Homotopy perturbation method for three-dimensional problem of condensation film on inclined rotating disk,”
*Scientia Iranica*, vol. 19, no. 3, pp. 437–442, 2012. View at Google Scholar - D. D. Ganji, H. Tari, and M. B. Jooybari, “Variational iteration method and homotopy perturbation method for nonlinear evolution equations,”
*Computers & Mathematics with Applications*, vol. 54, no. 7-8, pp. 1018–1027, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - Y. Khan, H. Vazquez-Leal, and N. Faraz, “An auxiliary parameter method using adomian polynomials and laplace transformation for nonlinear differential equations,”
*Applied Mathematical Modelling*. In press. - H. Vazquez-Leal, “Rational homotopy perturbation method,”
*Journal of Applied Mathematics*, vol. 2012, Article ID 490342, 14 pages, 2012. View at Publisher · View at Google Scholar - Y. Khan, H. Vazquez-Leal, and L. Hernandez-Martinez, “Removal of noise oscillation term appearing in the nonlinear equation solution,”
*Journal of Applied Mathematics*, vol. 2012, Article ID 387365, 9 pages, 2012. View at Publisher · View at Google Scholar