×

The series expansion and Chebyshev collocation method for nonlinear singular two-point boundary value problems. (English) Zbl 1476.65147

Summary: The solution of singular two-point boundary value problem is usually not sufficiently smooth at one or two endpoints of the interval, which leads to a great difficulty when the problem is solved numerically. In this paper, an algorithm is designed to recognize the singular behavior of the solution and then solve the equation efficiently. First, the singular problem is transformed to a Fredholm integral equation of the second kind via Green’s function. Second, the truncated fractional series of the solution about the singularity is formulated by using Picard iteration and implementing series expansion for the nonlinear function. Third, a suitable variable transformation is performed by using the known singular information of the solution such that the solution of the transformed equation is sufficiently smooth. Fourth, the Chebyshev collocation method is used to solve the deduced equation to obtain approximate solution with high precision. Fifth, the convergence analysis of the collocation method is conducted in weighted Sobolev spaces for linear singular equations. Sixth, numerical examples confirm the effectiveness of the algorithm. Finally, the Thomas-Fermi equation and the Emden-Fowler equation as some applications are accurately solved by the method.

MSC:

65L10 Numerical solution of boundary value problems involving ordinary differential equations
65L60 Finite element, Rayleigh-Ritz, Galerkin and collocation methods for ordinary differential equations
34B16 Singular nonlinear boundary value problems for ordinary differential equations
34B05 Linear boundary value problems for ordinary differential equations
65L20 Stability and convergence of numerical methods for ordinary differential equations
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] El-Gebeily, MA; Abu-Zaid, IT, On a finite difference method for singular two-point boundary value problems, IMA J Numer Anal, 18, 179-190 (1998) · Zbl 0913.65071
[2] Jamet, P., On the convergence of finite-difference approximations to one-dimensional singular boundary-value problems, Numer Math, 14, 355-378 (1970) · Zbl 0179.22103
[3] Kiguradze, IT; Shekhter, BL, Singular boundary-value problems for ordinary second-order differential equations, J Math Sci, 43, 2340-2417 (1988) · Zbl 0782.34026
[4] Lin, SH, Oxygen diffusion in a spherical cell with nonlinear oxygen uptake kinetics, J Theor Biol, 60, 449-457 (1976)
[5] Gray, BF, The distribution of heat sources in the human head-theoretical considerations, J Theor Biol, 82, 473-476 (1980)
[6] Wong, JSW, On the generalized Emden-Fowler equation, SIAM Rev, 17, 339-360 (1975) · Zbl 0295.34026
[7] Pikulin, SV, The Thomas-Fermi problem and solutions of the Emden-Fowler equation, Comput Math Math Phys, 59, 1292-1313 (2019) · Zbl 1430.34031
[8] Boyd, JP, Chebyshev spectral methods and the Lane-Emden problem, Numer Math Theor Meth Appl, 4, 142-157 (2011) · Zbl 1249.65158
[9] Hille, E., On the Thomas-Fermi equation, Proc Natl Acad Sci USA, 62, 7-10 (1969) · Zbl 0179.12903
[10] Hille, E., Some aspects of the Thomas-Fermi equation, J Anal Math, 23, 147-170 (1970) · Zbl 0208.11501
[11] Amore, P.; Boyd, JP; Fernández, FM, Accurate calculation of the solutions to the Thomas-Fermi equations, Appl Math Comput, 232, 929-943 (2014) · Zbl 1410.81015
[12] Flagg, RC; Luning, CD; Perry, WL, Implementation of new iterative techniques for solutions of Thomas-Fermi and Emden-Fowler equations, J Comput Phys, 38, 396-405 (1980) · Zbl 0456.65054
[13] Eriksson, K.; Thomée, V., Galerkin methods for singular boundary value problems in one space dimension, Math Comput, 42, 345-367 (1984) · Zbl 0546.65055
[14] Bender, CM; Orszag, SA, Advanced mathematical methods for scientists and engineers (1999), New York: Springer, New York · Zbl 0938.34001
[15] Ciarlet, PG; Natterer, F.; Varga, R., Numerical methods of high-order accuracy for singular nonlinear boundary value problems, Numer Math, 15, 87-99 (1970) · Zbl 0211.19103
[16] Chawla, MM; Katti, CP, Finite difference methods and their convergence for a class of singular two point boundary value problems, Numer Math, 39, 341-350 (1982) · Zbl 0489.65055
[17] Chawla, MM; Katti, CP, A uniform mesh finite difference method for a class of singular two-point boundary value problems, SIAM J Numer Anal, 22, 561-565 (1985) · Zbl 0578.65085
[18] Chawla, MM, A fourth-order finite difference method based on uniform mesh for singular two-point boundary-value problems, J Comput Appl Math, 17, 359-364 (1987) · Zbl 0614.65087
[19] Han, GQ; Wang, J.; Ken, H.; Xu, YS, Correction method and extrapolation method for singular two-point boundary value problems, J Comput Appl Math, 126, 145-157 (2000) · Zbl 0971.65066
[20] Abu-Zaid, IT; El-Gebeily, MA, A finite-difference method for the spectral approximation of a class of singular two-point boundary value problems, IMA J Numer Anal, 14, 545-562 (1994) · Zbl 0874.65059
[21] Kumar, M.; Aziz, T., A uniform mesh finite difference method for a class of singular two-point boundary value problems, Appl Math Comput, 180, 173-177 (2006) · Zbl 1103.65088
[22] Gustafsson, B., A numerical method for solving singular boundary value problems, Numer Math, 21, 328-344 (1973) · Zbl 0255.65032
[23] Kanth, ASVR; Reddy, YN, A numerical method for singular two point boundary value problems via Chebyshev economizition, Appl Math Comput, 146, 691-700 (2003) · Zbl 1035.65084
[24] Kadalbajoo, MK; Aggarwal, VK, Numerical solution of singular boundary value problems via Chebyshev polynomial and B-spline, Appl Math Comput, 160, 851-863 (2005) · Zbl 1062.65077
[25] Roul, P.; Goura, VMKP; Agarwal, R., A compact finite difference method for a general class of nonlinear singular boundary value problems with Neumann and Robin boundary conditions, Appl Math Comput, 350, 283-304 (2019) · Zbl 1429.65165
[26] Khuri, S.; Sayfy, A., A novel approach for the solution of a class of singular boundary value problems arising in physiology, Math Comput Model, 52, 626-636 (2010) · Zbl 1201.65135
[27] De Hoog, FR; Weiss, R., Collocation methods for singular boundary value problems, SIAM J Numer Anal, 15, 198-217 (1978) · Zbl 0398.65051
[28] Burkotová, J.; Rachunková, I.; Weinmüller, EB, On singular BVPs with nonsmooth data: convergence of the collocation schemes, BIT, 57, 1153-1184 (2017) · Zbl 1421.65021
[29] Auzinger, W.; Koch, O.; Weinmüller, EB, Analysis of a new error estimate for collocation methods applied to singular boundary value problems, SIAM J Numer Anal, 42, 2366-2386 (2005) · Zbl 1087.65082
[30] Huang, WZ; Ma, HP; Sun, WW, Convergence analysis of spectral collocation methods for a singular differential equation, SIAM J Numer Anal, 41, 2333-2349 (2003) · Zbl 1058.65076
[31] Russel, RD; Shampine, LF, Numerical methods for singular boundary value problems, SIAM J Numer Anal, 12, 13-36 (1975) · Zbl 0271.65051
[32] El-Gamel, M.; Sameeh, M., Numerical solution of singular two-point boundary value problems by the collocation method with the Chebyshev bases, SeMA, 74, 627-641 (2017) · Zbl 1380.65127
[33] Guo, BY; Wang, LL, Jacobi interpolation approximations and their applications to singular differential equations, Adv Comput Math, 14, 227-276 (2001) · Zbl 0984.41004
[34] Doedel, EJ; Reddien, GW, Finite difference methods for singular two-point boundary value problems, SIAM J Numer Anal, 21, 300-313 (1984) · Zbl 0565.65043
[35] Kumar, D., A collocation scheme for singular boundary value problems arising in physiology, Neural Parallel Sci Comput, 26, 95-118 (2018)
[36] Babolian, E.; Hosseini, MM, A modified spectral method for numerical solution of ordinary differential equations with non-analytic solution, Appl Math Comput, 132, 341-351 (2002) · Zbl 1024.65071
[37] Roul, P.; Thula, K., A new high-order numerical method for solving singular two-point boundary value problems, J Comput Appl Math, 343, 556-574 (2018) · Zbl 1453.65175
[38] Schreiber, R., Finite element methods of high-order accuracy for singular two-point boundary value problems with nonsmooth solutions, SIAM J Numer Anal, 17, 547-566 (1980) · Zbl 0474.65062
[39] Roul, P.; Warbhe, U., A novel numerical approach and its convergence for numerical solution of nonlinear doubly singular boundary value problems, J Comput Appl Math, 296, 661-676 (2016) · Zbl 1342.65164
[40] Roul, P.; Biswal, D., A new numerical approach for solving a class of singular two-point boundary value problems, Numer Algorithms, 75, 531-552 (2017) · Zbl 1385.65049
[41] Roul, P., A fast and accurate computational technique for efficient numerical solution of nonlinear singular boundary value problems, Int J Comput Math, 96, 51-72 (2019) · Zbl 1513.65237
[42] Roul, P.; Madduri, H., A new highly accurate domain decomposition optimal homotopy analysis method and its convergence for singular boundary value problems, Math Meth Appl Sci, 41, 6625-6644 (2018) · Zbl 1402.65087
[43] Singh, R., Analytic solution of singular Emden-Fowler-type equations by Green’s function and homotopy analysis method, Eur Phys J Plus, 134, 583, 1-17 (2019)
[44] Singh, R.; Kumar, J.; Nelakanti, G., Numerical solution of singular boundary value problems using Green’s function and improved decomposition method, J Appl Math Comput, 43, 409-425 (2013) · Zbl 1300.34036
[45] Singh, R.; Kumar, J., An efficient numerical technique for the solution of nonlinear singular boundary value problems, Comput Phys Commun, 185, 1282-1289 (2014) · Zbl 1344.65067
[46] Cen, ZD, Numerical study for a class of singular two-point boundary value problems using Green’s functions, Appl Math Comput, 18, 310-316 (2006) · Zbl 1105.65339
[47] Wang, TK; Liu, ZF; Zhang, ZY, The modified composite Gauss type rules for singular integrals using Puiseux expansions, Math Comput, 86, 345-373 (2017) · Zbl 1351.65018
[48] Wang, TK; Zhang, ZY; Liu, ZF, The practical Gauss type rules for Hadamard finite-part integrals using Puiseux expansions, Adv Comput Math, 43, 319-350 (2017) · Zbl 1373.65017
[49] Wang, TK; Gu, YS; Zhang, ZY, An algorithm for the inversion of Laplace transforms using Puiseux expansions, Numer Algorithms, 78, 107-132 (2018) · Zbl 1404.65316
[50] Hemmi, MA; Melkonian, S., Convergence of psi-series solutions of nonlinear ordinary differential equations, Can Appl Math Q, 3, 43-88 (1995) · Zbl 0835.34010
[51] Wang, TK; Li, N.; Gao, GH, The asymptotic expansion and extrapolation of trapezoidal rule for integrals with fractional order singularities, Int J Comput Math, 92, 579-590 (2015) · Zbl 1317.65066
[52] Liu, ZF; Wang, TK; Gao, GH, A local fractional Taylor expansion and its computation for insufficiently smooth functions, E Asian J Appl Math, 5, 176-191 (2015) · Zbl 1490.65002
[53] Zhao, TJ; Zhang, ZY; Wang, TK, A hybrid asymptotic and augmented compact finite volume method for nonlinear singular two point boundary value problems, Appl Math Comput, 392, 125745 (2021) · Zbl 1488.65206
[54] Shen, J.; Tang, T., Spectral and high-order methods with applications (2006), Beijing: Science Press, Beijing · Zbl 1234.65005
[55] Shen, J.; Tang, T.; Wang, LL, Spectral methods: algorithms, analysis and applications (2011), New York: Springer, New York · Zbl 1227.65117
[56] Canuto, C.; Hussaini, MY; Quarteroni, A.; Zang, TA, Spectral methods in fluid dynamics (1988), New York: Springer, New York · Zbl 0658.76001
[57] Quarteroni, A.; Valli, A., Numerical approximation of partial differential equations (2008), New York: Springer, New York · Zbl 1151.65339
[58] Rubio, G.; Fraysse, F.; Vicente, J.; Valero, E., The estimation of truncation error by \(\tau \)-estimation for Chebyshev spectral collocation method, J Sci Comput, 57, 146-173 (2013) · Zbl 1282.65125
[59] Baker, EB, The application of the Fermi-Thomas statistical model to the calculation of potential distribution in positive ions, Phys Rev, 36, 630-647 (1930) · JFM 56.1327.06
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.