Abstract
We studied a new implicit method of order three based on two off-step points for the numerical solution of second-order nonlinear initial-value problems \(y^{\prime\prime} = f(t,y,y^{\prime}),\;y(t_{0} ) = \gamma_{0}\), \(y^{\prime}(t_{0} ) = \gamma_{1}\) on a graded mesh. In practice, only a monotonically decreasing mesh will be employed. We applied the proposed method to a test equation \(y^{\prime\prime} + 2\alpha y^{\prime} + \beta^{2} y = g(t), \quad \alpha > \beta \ge 0,\) and analyse the stability which is absolutely stable for a graded mesh and superstable for a uniform mesh. The proposed method is applicable to solve singular problems. Seven benchmark problems including three boundary layer problems are solved to demonstrate the stability and utility of the method discussed. Numerical results are given to confirm the usefulness of the developed method.
![](http://media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs00366-019-00857-3/MediaObjects/366_2019_857_Fig1_HTML.png)
![](http://media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs00366-019-00857-3/MediaObjects/366_2019_857_Fig2_HTML.png)
![](http://media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs00366-019-00857-3/MediaObjects/366_2019_857_Fig3_HTML.png)
![](http://media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs00366-019-00857-3/MediaObjects/366_2019_857_Fig4_HTML.png)
![](http://media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs00366-019-00857-3/MediaObjects/366_2019_857_Fig5_HTML.png)
![](http://media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs00366-019-00857-3/MediaObjects/366_2019_857_Fig6_HTML.png)
![](http://media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs00366-019-00857-3/MediaObjects/366_2019_857_Fig7_HTML.png)
Similar content being viewed by others
References
Lambert JD (1973) Computational methods in ordinary differential equations. Wiley, London
Lambert JD (1991) Numerical methods for ordinary differential systems: the initial value problem. Wiley, London
Jain MK (2014) Numerical solution of differential equations: finite difference and finite element methods, 3rd edn. New AGE International Publication, New Delhi
Lambert JD, Watson IA (1976) Symmetric multistep methods for periodic initial value problems. J Inst Math Appl 18:189–202
Dahlquist G (1978) On accuracy and unconditional stability of linear multistep methods for second order differential equations. BIT 18:133–136
Hairer E (1979) Unconditionally stable methods for second order differential equations. Numer Math 32:373–379
Chawla MM (1983) Unconditionally stable Numerov type methods for second order differential equations. BIT 23:541–542
Chawla MM (1985) Superstable two-step methods for the numerical integration of general second order initial value problem. J Comput Appl Math 12:217–220
Coleman JP (2003) Order conditions for a class of two-step methods for y″ = f(x, y). IMA J Numer Anal 23:197–220
Rai AS, Ananthakrishnaiah U (1996) Additive parameters methods for the numerical integration of y″ = f(x, y, y′). J Comput Appl Math 67:271–276
Rai AS, Ananthakrishnaiah U (1997) Obrechkoff methods having additional parameters for general second order differential equations. J Comput Appl Math 79:167–182
Simos TE, Famelis IT, Tsitouras Ch (2003) Zero dissipative explicit Numerov-type methods for second order IVPs with oscillating solutions. Numer Algorithms 34:27–40
Saldanha G, Saldanha DJ (2009) A class of explicit two-step superstable methods for second-order linear initial value problems. Int J Comput Math 86:1424–1432
Lin T-C, Schultz DH, Zhang W (2008) Numerical solution of linear and non-linear singular perturbation problems. Comput Math Appl 55:2574–2592
Kadalbajoo MK, Kapil KK (2008) A numerical method based on finite difference for boundary value problems for singularly perturbed delay differential equations. Appl Math Comput 197:692–707
Kadalbajoo MK, Kumar D (2009) Initial value technique for singularly perturbed two point boundary value problems using an exponentially fitted finite difference scheme. Comput Math Appl 57:1147–1156
Sirajul-ul-Islam, Aziz I, Šarler B (2010) The numerical solution of second order boundary-value problems collocation method with the Haar wavelets. Math Comput Model 52:1577–1590
Sirajul-ul-Islam, Šarler B, Aziz I, Fazal-i-Haq (2011) Haar wavelet collocation method for the numerical solution of boundary layer fluid flow problems. Int J Therm Sci 50:686–697
Mohanty RK, McKee S (2015) On the stability of two new two-step explicit methods for the numerical integration of second order initial value problem on a variable mesh. Appl Math Lett 45:31–36
Samarskii AA (2001) The theory of difference schemes. Marcel Dekker, New York
Jain MK, Iyengar SRK, Subramanyam GS (1984) Variable mesh methods for the numerical solution of two point singular perturbation problems. Comput Methods Appl Mech Eng 42:273–286
Mohanty RK (2005) A family of variable mesh methods for the estimates of (du/dr) and the solution of nonlinear two point boundary value problems with singularity. J Comput Appl Math 182:173–187
Mohanty RK (2006) A class of non-uniform mesh three point arithmetic average discretization for y″ = f(x, y, y′) and the estimates of y′. Appl Math Comput 183:477–485
Mohanty RK (2004) An unconditionally stable difference scheme for the one space dimensional linear hyperbolic equation. Appl Math Lett 17:101–105
Mohanty RK (2009) New unconditionally stable difference schemes for the solution of multi-dimensional Telegraphic equations. Int J Comput Math 86:2061–2071
Evans DJ, Mohanty RK (2005) On the application of the SMAGE parallel algorithms on a non-uniform mesh for the solution of non-linear two point boundary value problems with singularity. Intern J Comput Math 82:341–353
Kelly CT (1995) Iterative methods for linear and non-linear equations. SIAM Publications, Philadelphia
Hageman LA, Young DM (2004) Applied iterative methods. Dover Publication, New York
Acknowledgements
The authors thank the reviewers for their valuable suggestions, which substantially improved the standard of the paper. The second author is supported by ‘The South Asian University’, Grant no.: SAU/AM(P)/2017/04.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Mohanty, R.K., Ghosh, B.P. Absolute stability of an implicit method based on third-order off-step discretization for the initial-value problem on a graded mesh. Engineering with Computers 37, 809–822 (2021). https://doi.org/10.1007/s00366-019-00857-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00366-019-00857-3
Keywords
- Graded mesh
- Damped wave equation
- Region of absolute stability
- Interval of periodicity
- Interval of weak stability
- Superstability
- Singular coefficients
- Boundary layer problems