Abstract
The approximate solution to second-order Hamiltonian and stiff differential systems is obtained using an efficient hybrid Nyström method (HNM) in this manuscript. The development of the method considers three hybrid points that are selected by optimizing the local truncation errors of the main formulas. The properties of the proposed HNM are studied. An embedding-like procedure is explored and run in variable step-size mode to improve the accuracy of the HNM. The numerical integration of some second-order Hamiltonian and stiff model problems, such as the well-known Vander Pol, Fermi-Pasta-Ulam, and Duffing problems, demonstrate the improved impact of our devised error estimation and control strategy. Finally, it is essential to note that the proposed technique is efficient in terms of computational cost and maximum global errors.
Similar content being viewed by others
References
Anastassi ZA, Kosti A (2015) A 6 (4) optimized embedded Runge–Kutta–Nyström pair for the numerical solution of periodic problems. J Comput Appl Math 275:311–320
Ascher UM, Petzold LR (1998) Computer methods for ordinary differential equations and differential-algebraic equations. Society for Industrial and Applied Mathematics (SIAM), Philadelphia
Atkinson KE, Han W, Stewart D (2009) Numerical solution of ordinary differential equations. Wiley, Hoboken
Brugnano L, Trigiante D (1998) Solving differential problems by multistep initial and boundary value methods. Gordon and Breach Science Publishers, London
Brugnano L, Iavernaro F, Trigiante D (2010) Hamiltonian boundary value methods (energy preserving discrete line integral methods). J Numer Anal Ind Appl Math 5(1–2):17–37
Brugnano L, Montijano JI, Randez L (2019) On the effectiveness of spectral methods for the numerical solution of multi-frequency highly oscillatory Hamiltonian problems. Numer Algorithm 81:345–376
Butcher JC (2008) Numerical methods for ordinary differential equations. Wiley, New York
Butcher JC (2016) Numerical methods for ordinary differential equations. Wiley, New York
Channell PJ, Scovel C (1990) Symplectic integration of Hamiltonian systems. Nonlinearity 3:231
Chen ZX, You X, Wei S, Liu ZL (2012) Symmetric and symplectic ERKN methods for oscillatory Hamiltonian systems. Comput Phys Commun 183:86–98
Feng K, Qin M (2010) Symplectic geometric algorithms for Hamiltonian systems. Springer, Berlin
Franco JM, Gómez I (2009) Accuracy and linear stability of RKN methods for solving second-order stiff problems. Appl Numer Math 59(5):959–975
Hairer E, Lubich C (2001) Long-time energy conservation of numerical methods for oscillatory differential equations. SIAM J Numer Anal 38:414–441
Hairer E, Norsett SP, Wanner G (1993) Solving ordinary differential equations I, Nonstiff problems. Springer, Berlin
Hairer E, Lubich C, Wanner G (2006) Geometric numerical integration: structure-preserving algorithms for ordinary differential equations, 2nd edn. Springer, Berlin
He Y, Qin H, Sun Y (2015) Hamiltonian integration methods for Vlasov–Maxwell equations. Phys Plasmas 22:124503
Henrici P (1962) Discrete variable methods in ODEs. Wiley, New York
Jator SN, King KL (2018) Integrating oscillatory general second-order initial value problems using a block hybrid method of order 11. Hindawi Math Probl Eng 2018:3750274
Jator SN, Li J (2009) A self starting linear multistep method for the direct solution of the general second order initial value problems. Int J Comput Math 86(5):817–836
Jator SN, Oladejo HB (2017) Block Nyström method for singular differential equations of the Lane–Emden type and problems with highly oscillatory solutions. Int J Appl Comput Math 3:1385–1402
Keskin AU (2018) Ordinary differential equations for engineers: problems with MATLAB solutions
Khalsaraei MM, Shokri A (2020) An explicit six-step singularly P-stable Obrechkoff method for the numerical solution of second-order oscillatory initial value problems. Numer Algorithm 84:871–886
Kosti AA, Anastassi ZA (2015) Explicit almost P-stable Runge–Kutta–Nyström methods for the numerical solution of the two-body problem. Comput Appl Math 34:647–659
Lambert JD (1991) Numerical methods for ordinary differential systems. Wiley, New York
Leimkuhler B, Reich S (2005) Simulating Hamiltonian dynamics. Cambridge University Press, Cambridge
Li J, Wu X (2019) Energy-preserving continuous stage extended Runge–Kutta–Nyström methods for oscillatory Hamiltonian systems Appl. Numer Math 145:469–487
Mazzia F, Iavernaro F (2003) Test set for initial value problems solvers. University of Bari, Department of Mathematics
Mei L, Wu X (2016) The construction of arbitrary order ERKN methods based on group theory for solving oscillatory Hamiltonian systems with applications. J Comput Phys 323:171–190
Paternoster B (1998) Runge–Kutta–Nyström methods for ODEs with periodic solutions based on trigonometric polynomials. Appl Numer Math 28:401–412
Ramos H, Rufai MA (2018) Third derivative modification of k-step block falkner methods for the numerical solution of second order initial-value problems, Journal of. Appl Math Comput 333:231–245
Ramos H, Rufai MA (2022) An adaptive one-point second-derivative Lobatto-type method for solving efficiently differential systems. Int J Comput Math 99(8):1687–1705
Rufai MA, Mazzia F, Ramos H (2023) An adaptive optimized Nyström method for second-order IVPs. Math Meth Appl Sci 46(6):7543–7556
Rufai MA, Ramos H (2020) Numerical solution of second-order singular problems arising in astrophysics by combining a pair of one-step hybrid block Nyström methods. Astrophys Space Sci 365:96
Rufai MA, Ramos H (2022) A variable step-size fourth-derivative hybrid block strategy for integrating third-order IVPs, with applications. Int J Comput Math 99(2):292–308
Shafie S, Tran T (2020) Numerical studies: adaptive mixed finite element schemes for Benjamin-bona-mahony and burgers equations. Int J Eng Trends Technol 1:63–69
Shi W, Wu X (2012) On symplectic and symmetric ARKN methods. Comput Phys Commun 183:1250–1258
Tang W, Sun Y (2012) Time finite element methods: a unified framework for numerical discretizations of ODEs. Appl Math Comput 219:2158–2179
Tang W, Zhang J (2018) Symplecticity-preserving continuous stage Runge–Kutta–Nyström methods. Appl Math Comput 323:204–219
Wu XY, Wang B, Xia JL (2012) Explicit symplectic multidimensional exponential fitting modified Runge–Kutta–Nyström methods. BIT Numer Math 52:773–795
You X, Zhao J, Yang H, Fang Y, Wu X (2014) Order conditions for RKN methods solving general second-order oscillatory systems. Numer Algorithm 66:147–176
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
On behalf of all authors, the corresponding author declares that there is no conflict of interest.
Additional information
Communicated by Jose Alberto Cuminato.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Rufai, M.A., Tran, T. & Anastassi, Z.A. A variable step-size implementation of the hybrid Nyström method for integrating Hamiltonian and stiff differential systems. Comp. Appl. Math. 42, 156 (2023). https://doi.org/10.1007/s40314-023-02273-2
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s40314-023-02273-2
Keywords
- Hybrid Nyström method
- Hamiltonian and stiff differential systems
- Variable stepsize formulation
- Error estimation and control
- Collocation method
Keywords
- Hybrid Nyström method
- Hamiltonian and stiff differential systems
- Variable stepsize formulation
- Error estimation and control
- Collocation method