IJEM International Journal of Engineering Mathematics 2314-6109 2356-7007 Hindawi Publishing Corporation 10.1155/2015/949463 949463 Research Article Successive Complementary Expansion Method for Solving Troesch’s Problem as a Singular Perturbation Problem http://orcid.org/0000-0002-4345-1253 Cengizci Süleyman 1 http://orcid.org/0000-0002-2794-6715 Eryılmaz Aytekin 2 Dandapat Bhabani S. 1 Insitute of Applied Mathematics Middle East Technical University 06800 Ankara Turkey metu.edu.tr 2 Department of Mathematics Nevsehir Haci Bektas Veli University 50300 Nevsehir Turkey nevsehir.edu.tr 2015 26102015 2015 08 05 2015 04 10 2015 26102015 2015 Copyright © 2015 Süleyman Cengizci and Aytekin Eryılmaz. 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.

A simple and efficient method that is called Successive Complementary Expansion Method (SCEM) is applied for approximation to an unstable two-point boundary value problem which is known as Troesch’s problem. In this approach, Troesch’s problem is considered as a singular perturbation problem. We convert the hyperbolic-type nonlinearity into a polynomial-type nonlinearity using an appropriate transformation, and then we use a basic zoom transformation for the boundary layer and finally obtain a nonlinear ordinary differential equation that contains SCEM complementary approximation. We see that SCEM gives highly accurate approximations to the solution of Troesch’s problem for various parameter values. Moreover, the results are compared with Adomian Decomposition Method (ADM) and Homotopy Perturbation Method (HPM) by using tables.

1. Introduction

Troesch’s highly sensitive problem arises from a system of a nonlinear ordinary differential equations which occur in the investigation of the confinement of a plasma column by radiation pressure . It also arises in the theory of gas porous electrodes [2, 3]. The problem is defined by(1)y=λsinhλy,0x1,with the boundary conditions(2)y0=0,y1=1,where y=y(x) and λ>0. Roberts and Shipman  have shown that the closed form solution to problem (1) with the boundary conditions (2) in terms of the Jacobi elliptic function scnm is as follows:(3)yx=2λsinh-1y02scλx1-14y02,where y(0)=21-m. The constant m satisfies the transcendental equation(4)scλm=sinhλ/21-m,where scλm is the Jacobi function defined by scλm=sinϕ/cosϕ=tanϕ. Here ϕ,λ, and m are related by the integral(5)λ=0ϕdθ1-msin2θ.It has been shown in [4, 5] that y(x) has a singularity approximately located at(6)xs12λln161-m=1λln8y0.This singularity makes the problem very difficult to solve for large λ values.

The first explanation of Troesch’s problem was given and solved by Weibel . Because of its physical significance, Troesch’s problem has always been attractive to scientists and various methods have been implemented to solve it such as Monte Carlo method , Modified Newton Method [7, 8], Transformation Groups Method , Invariant Embedded Method , Shooting Method , Inverse Shooting Method , The quasilinearization Method [12, 13], and Matched Asymptotic Expansions Method . Recently, approximate solution techniques such as Adomian Decomposition Method [15, 16], Laplace transform and a modified decomposition technique , Variational Iteration Method [17, 18], Homotopy Perturbation Method , Differential Transform Method , B-spline collocation approach , The Sinc-Galerkin Method , Wavelet Analysis Method , Sinc-Collocation Method , and Jacobi Collocation Method  have been successfully applied to Troesch’s problem by various researchers. In  more detailed information and references on Troesch’s problem may also be seen.

In this paper, we consider Troesch’s problem as a singular perturbation problem. By using appropriate transformations, we apply Successive Complementary Expansion Method (SCEM) to it. The results obtained by SCEM are compared with HPM and ADM and the analytic solutions.

2. The Successive Complementary Expansion Method

In this section, we first give a short overview of asymptotic approximations. One can consult  for some definitions on the asymptotic approximations and singular perturbations in more detail.

SCEM was first introduced by Mauss and Cousteix (see ). It has always been applied to partial differential equations except for some introductory examples in . SCEM is based on the knowledge of generalized asymptotic expansions and the Method of Matched Asymptotic Expansions (MMAE). At first, a structure of the uniformly valid approximation (UVA) must be assumed and then the method to construct the UVA is deduced [32, 33]. In this method, any matching procedure is required in contrast to MMAE. Furthermore, the boundary conditions are satisfied exactly but not asymptotically.

Let us consider a sequence of real numbers f(n)n=1. We are familiar with the limit notation limnf(n). If this limit exists we can denote it as limnf(n)=fR. We can give the following definition: given any ε>0 there exists a number n0(ε) such that f(n)-f<ε for any nn0(ε). This definition contains information about the behaviour of f(n) as n but not about how f(n) approaches f. In order to describe the behaviour of sequences in a more precise way, we need some other definitions, Bachmann-Landau notations. Let us consider two continuous functions of real numbers that depend on a small parameter εf(ε)=O(g(ε)) for ε0 if there exists positive constants C and ε0 such that, in (0,ε0], f(ε)Cg(ε) for ε0. f(ε)=o(g(ε)) for ε0 if limε0f(ε)/g(ε)=0. f(ε)=OS(g(ε)) if f(ε)=O(g(ε)) and f(ε)o(g(ε)) as ε0. Let E be a set of real functions that depend on ε, strictly positive and continuous in (0,ε0], such that limε0δ(ε) exists and if δ1 and δ2E, δ1δ2E. A function that satisfies these conditions is called order function. Given two functions ϕ(x,ε) and ϕa(x,ε) defined in a domain Ω, they are asymptotically identical to order δε if their difference is asymptotically smaller than δε, where δε is an order function; that is,(7)ϕx,ε-ϕax,ε=oδε,where ε is small parameter arising from the physical problem under consideration. The function ϕa(x,ε) is named as asymptotic approximation of the function ϕ(x,ε). Asymptotic approximations in general form are defined by(8)ϕax,ε=i=1nδiεφix,ε,where the asymptotic sequence of order functions δi(ε) is an asymptotic sequence such that δi+1(ε)=o(δi(ε)), as ε0. In these conditions approximation (8) is named as generalized asymptotic expansion. If expansion (8) is written in the form of(9)ϕax,ε=E0ϕ=i=1nδi0εφi0x,then it is called regular asymptotic expansion, where the special operator E0 is outer expansion operator at a given order δ(ε). Thus ϕ-E0ϕ=o(δ(ε)). For more detailed information about the asymptotic approximations, [10, 11, 13, 19] can be studied. Interesting cases occur when the function is not regular in Ω so (8) or (9) is valid only in a restricted region Ω0Ω called the outer region. We are faced with a singular perturbation problem and we must introduce boundary layer domains. We introduce an inner domain which can be formally denoted as Ω1=Ω0-Ω and located near the point x=1 (for Troesch’s problem). The boundary layer variable is x¯=x-1/ξ(ε), with ξ(ε) being the order of thickness of this boundary layer. If a regular expansion can be constructed in Ω1, we can write down(10)ϕax,ε=E1ϕ=i=0nδi1εφi1x¯,where the inner expansion operator E1 is defined in Ω1 at the same order δ(ε) as the outer expansion operator E0; thus, ϕ-E1ϕ=o(δ(ε)) and(11)ϕa=E0ϕ+E1ϕ-E1E0ϕis clearly uniformly valid approximation (UVA) [35, 36]. In the Method of Matched Asymptotic Expansions (MMAE) , two distinct solutions of the singular perturbation problem are found for two distinct regions (inner and outer) and then to obtain uniformly valid approximation over the whole domain the solutions are matched using limit process(12)limx1y0x=limx¯-Yx¯,where Y(x¯) is boundary layer approximation. Finally uniformly valid approximation is obtained by simply adding the inner and outer solutions and subtracting the common limit. That is, using the procedure(13)yy0x+Y0x¯-Y0-,or equivalently(14)yy0x+Y0x¯-y01-,one can reach the composite solution which is uniformly valid over the whole domain. The uniformly valid SCEM approximation is in the regular form(15)ynscemx,x¯,ϵ=i=1nδiεyix+Ψix¯,where δi is an asymptotic sequence and Ψi(x¯) are the complementary functions that depend on x¯. Functions yi(x) are the outer approximation functions that have been found by MMAE and they only depend on x, not also ε. If the functions yi(x) and Ψi(x¯) depend also on ε, the uniformly valid SCEM approximation is called generalized SCEM approximation that is in the following form [30, 31]:(16)yngscemx,x¯,ϵ=i=1nδiεyix,ε+Ψix¯,ε.For the sake of simplicity, we adopt just one-term approximation in this study. That is, we look for an approximation in the form of(17)y0scemx,x¯,ε=y0x,ε+Ψ0x¯,ε.To improve the accuracy of approximation, (17) can be iterated using (16). It means that successive complementary terms will be added to the approximation.

3. Application of SCEM to Troesch’s Problem

In this section, we solve Troesch’s problems for different values of the parameter λ using the computer algebra system Matlab and make a comparison between our results and those ones reported in the literature to confirm the efficiency and accuracy of our method. In first step, let us turn Troesch’s problem into singular perturbation problem as follows:(18)y=λsinhλy.Letting λ=1/ε in (18), we obtain(19)y=1εsinhyε,εy-sinhyε=0.

Now, in order to avoid overflow or excessive error growth during numerical integration, removing the hyperbolic-type nonlinearity using the variable transformation, which was originally proposed in , y(x)=4εtanh-1(u(x)) or equivalently u(x)=tanhy(x)/4ε, we get(20)y=4ε1-u2u,y=4ε1-u2u+8εu1-u22u2or(21)1-u2u+2uu2=u1+u2ε2and finally Troesch’s problem has the form(22)ε2u+2ε2uu21-u2-u1+u21-u2=0with the boundary conditions(23)u0=0,u1=tanh14ε.We reach polynomial-type nonlinearity instead of hyperbolic-type nonlinearity. The problem has the outer solution (near the point x=0) y0(x)=0 (or equivalently u0(x)=0). So the complementary solution directly produces approximation to Troesch’s problem. If we substitute the inner variable x¯=x-1/ε into (22) and boundary conditions (23), using the chain rule, we reach(24)Ψ0+2Ψ0Ψ02-Ψ01+Ψ021-Ψ02=0and equivalently(25)Ψ0+Ψ02Ψ02-Ψ02-11-Ψ02=0with the boundary conditions(26)Ψ0-1ε=0,Ψ00=tanh14ε,where Ψ0 is the complementary function. Here, the complementary function of SCEM, Ψ0, is calculated numerically by Matlab bvp4c routine. Results which are obtained using SCEM are presented in the following section. Applying SCEM, we reach really highly accurate approximations to Troesch’s problem.

4. Conclusions

In this study, an efficient method called Successive Complementary Expansion Method (SCEM) is used to solve the nonlinear two-point boundary value problem with application to Troesch’s equation. SCEM is an easy-applicable and efficient method. It does not require any matching procedure in contrast to MMAE. Moreover, the boundary conditions are satisfied exactly, not asymptotically. In Tables 13, the absolute errors in solutions obtained by the presented method for ε=2, ε=1, and ε=0.1, respectively, are shown. In Figures 15, we compare exact and SCEM solutions. The results obtained here were compared with the exact solution, ADM , and HPM . It is shown that the method is computationally attractive and the presented method is much better than other reported ones in the literature in the sense of accuracy and efficiency.

Solution of Troesch’s problem for ε=2.

x Exact solution SCEM ADM  HPM  Abs. err. in SCEM
0.1 0.09594435 0.09594435 0.09593835 0.09593956 3.13441692 e ( - 9 )
0.2 0.19212874 0.19212875 0.19211805 0.19211932 7.32982755 e ( - 9 )
0.3 0.28879440 0.28879441 0.28878032 0.28878069 1.30811447 e ( - 8 )
0.4 0.38618485 0.38618487 0.38616870 0.38616754 1.98750830 e ( - 8 )
0.5 0.48454716 0.48454719 0.48453029 0.48452741 2.60076459 e ( - 8 )
0.6 0.58413325 0.58413328 0.58411697 0.58411278 2.88212714 e ( - 8 )
0.7 0.68520115 0.68520118 0.68518684 0.68518224 2.55610896 e ( - 8 )
0.8 0.78801652 0.78801654 0.78800556 0.78800183 1.51078886 e ( - 8 )
0.9 0.89285422 0.89285422 0.89284802 0.89284621 9.280242130 e ( - 10 )

Solution of Troesch’s problem for ε=1.

x Exact solution SCEM ADM  HPM  Abs. err. in SCEM
0.1 0.08466141 0.08466143 0.08424876 0.08438170 2.18232920 e ( - 8 )
0.2 0.17017167 0.17017172 0.16943070 0.16962076 5.16987398 e ( - 8 )
0.3 0.25739438 0.25739447 0.25641450 0.25659292 8.76146832 e ( - 8 )
0.4 0.34722348 0.34722359 0.34608572 0.34621073 1.09822865 e ( - 7 )
0.5 0.44060059 0.44060067 0.43940198 0.43944227 7.88214280 e ( - 8 )
0.6 0.53853523 0.53853518 0.53736570 0.53733006 5.56789394 e ( - 8 )
0.7 0.64212944 0.64212911 0.64108380 0.64101046 3.25396753 e ( - 7 )
0.8 0.75260880 0.75260812 0.75178800 0.75173354 6.76312419 e ( - 7 )
0.9 0.87136294 0.87136211 0.87090870 0.87088353 8.29947601 e ( - 7 )

Absolute error of SCEM approximation for ε=0.1.

x Exact solution SCEM Abs. err. in SCEM
0.10 0.000004211 0.000042111 3.79006047 e ( - 5 )
0.20 0.000129963 0.000129964 1.57300000 e ( - 10 )
0.30 0.000358977 0.000358978 4.20000000 e ( - 10 )
0.40 0.000977901 0.000977902 1.08799999 e ( - 9 )
0.50 0.002659017 0.002659019 2.60240000 e ( - 9 )
0.60 0.007228924 0.007228929 4.58729999 e ( - 9 )
0.70 0.019664060 0.019664057 2.41529999 e ( - 9 )
0.80 0.053730329 0.053730314 1.46944000 e ( - 8 )
0.90 0.152114078 0.152114033 4.55675999 e ( - 8 )
0.95 0.276267734 0.276267643 9.12888999 e ( - 8 )
0.99 0.574076498 0.574076080 4.17938999 e ( - 7 )

Exact solution and SCEM approximation of Troesch’s problem for ε=2.

Exact solution and SCEM approximation of Troesch’s problem for ε=1.

Exact solution and SCEM approximation of Troesch’s problem for ε=0.5.

Exact solution and SCEM approximation of Troesch’s problem for ε=0.1.

Exact solution and SCEM approximation of Troesch’s problem for ε=0.05.

Conflict of Interests

The authors declare that they have no conflict of interests regarding the publication of this paper.

Weibel E. S. On the confinement of a plasma by magnetostatic fields Physics of Fluids 1959 2 1 52 56 10.1063/1.1724391 2-s2.0-0039684959 Markin V. S. Chernenko A. A. Chizmadehev Y. A. Chirkov Y. G. Bagotskii V. S. Vasilev Y. B. Aspects of the theory of gas porous electrodes Fuel Cells: Their Electrochemical Kinetics 1966 New York, NY, USA Consultants Bureau 21 33 Gidaspow D. Baker B. S. A model for discharge of storage batteries Journal of the Electrochemical Society 1973 120 8 1005 1010 10.1149/1.2403617 2-s2.0-0015658471 Roberts S. M. Shipman J. S. On the closed form solution of Troesch's problem Journal of Computational Physics 1976 21 3 291 304 10.1016/0021-9991(76)90026-7 ZBL0334.65062 2-s2.0-0040384921 Troesch B. A. A simple approach to a sensitive two-point boundary value problem Journal of Computational Physics 1976 21 3 279 290 10.1016/0021-9991(76)90025-5 2-s2.0-0345899177 Tsuda T. Ichida K. Kiyono T. Monte Carlo path-integral calculations for two-point boundary-value problems Numerische Mathematik 1967 10 2 110 116 10.1007/BF02174142 2-s2.0-2942535457 Jones D. J. Solution of Troesch's, and other, two point boundary value problems by shooting techniques Journal of Computational Physics 1973 12 3 429 434 10.1016/0021-9991(73)90165-4 2-s2.0-33748111182 Kubiček M. Hlaváček V. Solution of Troesch's two-point boundary value problem by shooting technique Journal of Computational Physics 1975 17 1 95 101 10.1016/0021-9991(75)90066-2 Chiou J. P. Na T. Y. On the solution of Troesch's nonlinear two-point boundary value problem using an initial value method Journal of Computational Physics 1975 19 3 311 316 10.1016/0021-9991(75)90080-7 2-s2.0-49549147821 Scott M. R. Aziz A. K. On the conversion of boundary-value problems into stable initial-value problems via several invariant imbedding algorithms Numerical Solutions of Boundary-Value Problems for Ordinary Differential Equations 1975 New York, NY, USA Academic Press 89 146 Snyman J. A. Continuous and discontinuous numerical solutions to the Troesch problem Journal of Computational and Applied Mathematics 1979 5 3 171 175 10.1016/0377-0427(79)90002-5 ZBL0419.65047 2-s2.0-0018521558 Miele A. Aggarwal A. K. Tietze J. L. Solution of two-point boundary-value problems with Jacobian matrix characterized by large positive eigenvalues Journal of Computational Physics 1974 15 2 117 133 10.1016/0021-9991(74)90080-1 2-s2.0-45149092639 Vemuri V. Raefsky A. On a method of solving sensitive boundary value problems Journal of the Franklin Institute 1979 307 4 217 243 10.1016/0016-0032(79)90049-8 2-s2.0-0018452959 Chin R. C. Y. On the matched asymptotic solution of the Troesch problem Journal of Computational and Applied Mathematics 1981 7 3 181 185 10.1016/0771-050x(81)90016-4 2-s2.0-49049151579 Deeba E. Khuri S. A. Xie S. An algorithm for solving boundary value problems Journal of Computational Physics 2000 159 2 125 138 10.1006/jcph.2000.6452 2-s2.0-0005338184 Khuri S. A. A numerical algorithm for solving Troesch's problem International Journal of Computer Mathematics 2003 80 4 493 498 10.1080/0020716022000009228 2-s2.0-28244446638 Momani S. Abuasad S. Odibat Z. Variational iteration method for solving nonlinear boundary value problems Applied Mathematics and Computation 2006 183 2 1351 1358 10.1016/j.amc.2006.05.138 2-s2.0-33845753056 Chang S.-H. A variational iteration method for solving Troesch’s Problem Journal of Computational and Applied Mathematics 2010 234 10 3043 3047 10.1016/j.cam.2010.04.018 2-s2.0-77955279417 Feng X. Mei L. He G. An efficient algorithm for solving Troesch's problem Applied Mathematics and Computation 2007 189 1 500 507 10.1016/j.amc.2006.11.161 2-s2.0-34248395626 Chang S.-H. Chang I.-L. A new algorithm for calculating one-dimensional differential transform of nonlinear functions Applied Mathematics and Computation 2008 195 2 799 808 10.1016/j.amc.2007.05.026 2-s2.0-37349013500 Khuri S. A. Sayfy A. Troesch's problem: a B-spline collocation approach Mathematical and Computer Modelling 2011 54 9-10 1907 1918 10.1016/j.mcm.2011.04.030 2-s2.0-80051583591 Zarebnia M. Sajjadian M. The sinc-galerkin method for solving Troesch's problem Mathematical and Computer Modelling 2012 56 9-10 218 228 10.1016/j.mcm.2011.11.071 2-s2.0-84863599859 Nasab A. K. Atabakan Z. P. Kılıçman A. An efficient approach for solving nonlinear Troesch's and Bratu's problems by wavelet analysis method Mathematical Problems in Engineering 2013 2013 10 825817 10.1155/2013/825817 2-s2.0-84884243573 El-Gamel M. Numerical solution of Troesch's problem by sinc-collocation method Applied Mathematics 2013 4 4 707 712 10.4236/am.2013.44098 Doha E. H. Baleanu D. Bhrawi A. H. Hafez R. M. A Jacobi collocation method for Troesch's problem in plasma physics Proceedings of the Romanian Academy, Series A 2014 15 2 130 138 Bougoffa L. Al-khadhi M. A. New explicit solutions for Troeschs boundary value problem Applied Mathematics & Information Sciences 2009 3 1 8996 Zbl1191.34025 Mirmoradi S. H. Hosseinpour I. Ghanbarpour S. Barari A. Application of an approximate analytical method to nonlinear Troesch's problem Applied Mathematical Sciences 2009 3 29–32 1579 1585 2-s2.0-77951003631 Hassan H. N. El-Tawil M. A. An efficient analytic approach for solving two-point nonlinear boundary value problems by homotopy analysis method Mathematical Methods in the Applied Sciences 2011 34 8 977 989 10.1002/mma.1416 2-s2.0-79955804712 Mohyud-Din S. T. Solution of Troesch's problem using He's polynomials Revista de la Unión Matemática Argentina 2011 52 1 143 148 2-s2.0-80052465750 Cousteix J. Mauss J. Asymptotic Analysis and Boundary Layers 2007 Berlin, Germany Springer Scientific Computation 10.1007/978-3-540-46489-1 Cousteix J. Mauss J. Approximations of the Navier-Stokes equations for high Reynolds number flows past a solid wall Journal of Computational and Applied Mathematics 2004 166 1 101 122 10.1016/j.cam.2003.09.035 ZBL1107.76024 2-s2.0-1542634835 Mauss J. Cousteix J. Uniformly valid approximation for singular perturbation problems and matching principle Comptes Rendus Mécanique 2002 330 10 697 702 10.1016/s1631-0721(02)01522-x 2-s2.0-0036812613 Cousteix J. Mauss J. Interactive boundary layers in turbulent flow Comptes Rendus Mécanique 2007 335 9-10 590 605 10.1016/j.crme.2007.08.014 2-s2.0-35448968754 Mauss J. On matching principles Asymptotic Analysis 1979 711 Berlin, Germany Springer 1 8 Lecture Notes in Mathematics 10.1007/BFb0062944 Cathalifaud P. Mauss J. Cousteix J. Nonlinear aspects of high Reynolds number channel flows European Journal of Mechanics—B/Fluids 2010 29 4 295 304 10.1016/j.euromechflu.2010.02.002 2-s2.0-77955275315 Lagerstrom P. A. Matched Asymptotic Expansions: Ideas and Techniques 1988 76 New York, NY, USA Springer Applied Mathematical Sciences 10.1007/978-1-4757-1990-1