Skip to main content

Advertisement

SpringerLink
Log in

Energy-preserving trigonometrically fitted continuous stage Runge-Kutta-Nyström methods for oscillatory Hamiltonian systems

  • Original Paper
  • Published:
Numerical Algorithms Aims and scope Submit manuscript

Abstract

Recently, continuous-stage Runge-Kutta-Nyström (CSRKN) methods for solving numerically second-order initial value problem \(q^{\prime \prime }= f(q)\) have been proposed and developed by Tang and Zhang (Appl. Math. Comput. 323, 204–219, 2018). This problem is equivalent to a separable Hamiltonian system when f(q) = −∇U(q) with smooth function U(q). Symplecticity-preserving discretizations of this system were studied in that paper. However, as an important representation of the Hamiltonian system, energy preservation has not been studied. In addition, many Hamiltonian systems in practical applications often have oscillatory characteristics so we should design special algorithms adapted to this feature. In this paper, we propose and study energy-preserving trigonometrically fitted CSRKN methods for oscillatory Hamiltonian systems. We extend the theory of trigonometrical fitting to CSRKN methods and derive sufficient conditions for energy preservation. We also study the symmetry and stability of the methods. Two symmetric and energy-preserving trigonometrically fitted schemes of order two and four, respectively, are constructed. Some numerical experiments are provided to confirm the theoretical expectations.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Log in via an institution

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Similar content being viewed by others

References

  1. Ixaru, L.G.: Exponential Fitting. Kluwer Academic Publishers, Dordrecht (2004)

    Book  MATH  Google Scholar 

  2. Gautschi, W.: Numerical integration of ordinary differential equation based on trigonometric polynomials. Numer. Math. 3, 381–397 (1961)

    Article  MathSciNet  MATH  Google Scholar 

  3. Lyche, T.: Chebyshevian multistep methods for ordinary differential equations. Numer. Math. 19, 65–75 (1972)

    Article  MathSciNet  MATH  Google Scholar 

  4. Simos, T.E.: A family of fifth algebraic order trigonometrically fitted Runge-Kutta methods for the numerical solution of the Schrödinger equation. Comput. Mater. Sci. 34, 342–354 (2005)

    Article  Google Scholar 

  5. Simos, T.E., Aguiar, J.V.: A modified phase-fitted Runge-Kutta method for the numerical solution of the Schrödinger equation. J. Math. Chem. 30, 121–131 (2001)

    Article  MathSciNet  MATH  Google Scholar 

  6. Franco, J.M.: Runge-Kutta methods adapted to the numerical integration of oscillatory problems. Numer. Math. Appl. 50, 427–443 (2004)

    Article  MathSciNet  MATH  Google Scholar 

  7. Van de Vyver, H.: An embedded exponentially fitted Runge-Kutta-Nyström method for the numerical solution of orbital problems. New Astron. 11, 577–587 (2006)

  8. Vanden Berghe, G., De Meyer, H., Van Daele, M., Van Hecke, T.: Exponentially fitted Runge-Kutta methods. Comput. Phys. Commun. 123, 7–15 (1999)

    Article  MathSciNet  MATH  Google Scholar 

  9. Vanden Berghe, G., De Meyer, H., Van Daele, M., Van Hecke, T.: Exponentially fitted explicit Runge-Kutta methods. J. Comput. Appl. Math. 125, 107–115 (2000)

    Article  MathSciNet  MATH  Google Scholar 

  10. Paternoster, B.: Runge-Kutta(-Nystrom) methods for ODEs with periodic solutions based on trigonometric polynomials. Appl. Numer. Math. 28, 401–412 (1998)

    Article  MathSciNet  MATH  Google Scholar 

  11. Van de Vyver, H.: Scheifele two-step methods for perturbed oscillators. J. Comput. Appl. Math. 224, 415–432 (2009)

  12. Li, J.Y., Wang, B., You, X., Wu, X.Y.: Two-step extended RKN methods for oscillatory systems. Comput. Phys. Commun. 182, 2486–2507 (2011)

    Article  MathSciNet  MATH  Google Scholar 

  13. Li, J.Y., Deng, S.: Trigonometrically fitted multi-step RKN methods for second-order oscillatory initial value problems. Appl. Math. Comput. 320, 740–753 (2018)

    MathSciNet  MATH  Google Scholar 

  14. Li, J.Y., Wang, X.F., Deng, S., Wang, B.: Symmetric trigonometrically-fitted two-step hybrid methods for oscillatory problems. J. Comput. Appl. Math. 344, 115–131 (2018)

    Article  MathSciNet  MATH  Google Scholar 

  15. Li, J.Y., Deng, S., Wang, X.F.: Extended explicit pseudo two-step RKN methods for oscillatory systems \(y^{\prime \prime } + My =f(y)\). Numer. Algor. 78, 673–700 (2018)

    Article  MathSciNet  MATH  Google Scholar 

  16. Li, J.Y., Wu, X.Y.: Adapted Falkner-type methods solving oscillatory second-order differential equations. Numer. Algor. 62, 355–381 (2013)

    Article  MathSciNet  MATH  Google Scholar 

  17. Wang, B., Iserles, A., Wu, X.: Arbitrary order trigonometric Fourier collocation methods for second-order ODEs. Found. Comput. Math. 16, 151–181 (2016)

    Article  MathSciNet  MATH  Google Scholar 

  18. Wang, B., Wu, X., Meng, F.: Trigonometric collocation methods based on Lagrange basis polynomials for multi-frequency oscillatory second-order differential equations. J. Comput. Appl. Math. 313, 185–201 (2017)

    Article  MathSciNet  MATH  Google Scholar 

  19. Zhong, G., Marsden, J.E.: Lie-Poisson Hamilton-Jacobi theory and Lie-Poisson integrators. Phys. Lett. A 133, 134–139 (1988)

    Article  MathSciNet  MATH  Google Scholar 

  20. Chartier, P., Faou, E., Murua, A.: An algebraic approach to invariant preserving integators: the case of quadratic and Hamiltonian invariants. Numer. Math. 103, 575–590 (2006)

    Article  MathSciNet  MATH  Google Scholar 

  21. Hairer, E., Lubich, C., Wanner, G.: Geometric Numerical Integration, Structure-Preserving Algorithms for Ordinary Differential Equations, 2nd edn. Springer, Berlin (2006)

    MATH  Google Scholar 

  22. Celledoni, E., McLachlan, R.I., Owren, B., Quispel, G.R.W.: Energy-preserving integrators and the structure of B-series. Tech. rep., NTNU preprint series: Numerics No.5/2009

  23. Hairer, E.: Energy-preserving variant of collocation methods. JNAIAM 5, 73–84 (2010)

    MathSciNet  MATH  Google Scholar 

  24. Miyatake, Y., Butcher, J.C.: Characterization of energy-preserving methods and the construction of parallel integrators for Hamiltonian systems. SIAM J. Numer. Anal. 54, 1993–2013 (2016)

    Article  MathSciNet  MATH  Google Scholar 

  25. Tang, W.S., Sun, Y.J., Cai, W.J.: Discontinuous Galerkin methods for Hamiltonian ODEs and PDEs. J. Comput. Phys. 330, 340–364 (2017)

    Article  MathSciNet  MATH  Google Scholar 

  26. Celledoni, E., Owren, B., Sun, Y.J.: The minimal stage, energy preserving Runge-Kutta method for polynomial Hamiltonian systems is the averaged vector field method. Math. Comput. 83, 1689–1700 (2014)

    Article  MathSciNet  MATH  Google Scholar 

  27. Tang, W.S., Zhang, J.J.: Symplecticity-preserving continuous-stage Runge-Kutta-Nyström methods. Appl. Math. Comput. 323, 204–219 (2018)

    MathSciNet  MATH  Google Scholar 

  28. Hairer, E., Nørsett, S.P., Wanner, G.: Solving Ordinary Differential Equations I: Nonstiff Problems, 2nd edn. Springer, Berlin (2002)

    MATH  Google Scholar 

  29. Butcher, J.C.: Numerical Methods for Ordinary Differential Equations, 2nd edn. Wiley, Chichester (2008)

    Book  Google Scholar 

  30. Butcher, J.C.: An algebraic theory of integration methods. Math. Comput. 26, 79–106 (1972)

    Article  MathSciNet  MATH  Google Scholar 

  31. Hairer, E., Wanner, G.: A theory for Nyström methods. Numer. Math. 25, 383–400 (1976)

    Article  MATH  Google Scholar 

  32. Franco, J.M.: Exponentially fitted symplectic integrators of RKN type for solving oscillatory problems. Comput. Phys. Commun. 177, 479–492 (2007)

    Article  MathSciNet  MATH  Google Scholar 

  33. Hoang, N. S., Sidje, R. B.: Functionally fitted Runge-Kutta-Nyström methods. BIT Numer. Math. 56, 129–150 (2016)

    Article  MATH  Google Scholar 

  34. Miyatake, Y.: An energy-preserving exponentially-fitted continuous stage Runge-Kutta method for Hamiltonian systems. BIT Numer. Math. 54, 1–23 (2014)

    Article  MathSciNet  MATH  Google Scholar 

  35. Ramos, H., Vigo-Aguiar, J.: On the frequency choice in trigonometrically fitted methods. Appl. Math. Lett. 23, 1378–1381 (2010)

    Article  MathSciNet  MATH  Google Scholar 

  36. Vigo-Aguiar, J., Ramos, H.: On the choice of the frequency in trigonometrically-fitted methods for periodic problems. J. Comput. Appl. Math. 277, 94–105 (2015)

    Article  MathSciNet  MATH  Google Scholar 

Download references

Acknowledgements

The authors are sincerely thankful to the anonymous referees for their constructive comments and valuable suggestions.

Funding

The research was supported in part by the Natural Science Foundation of China under Grant No: 11401164, by Hebei Natural Science Foundation of China under Grant No: A2014205136 and by Science Foundation of Hebei Normal University No:L2018J01.

Author information

Authors and Affiliations

  1. College of Mathematics and Information Science, Hebei Normal University, Shijiazhuang, 050024, People’s Republic of China

    Jiyong Li & Yachao Gao

  2. Hebei Key Laboratory of Computational Mathematics and Applications, Shijiazhuang, 050024, People’s Republic of China

    Jiyong Li & Yachao Gao

Authors
  1. Jiyong Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yachao Gao
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Jiyong Li.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

The research was supported in part by the Natural Science Foundation of China under Grant No. 11401164 and by Hebei Natural Science Foundation of China under Grant No. A2014205136.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, J., Gao, Y. Energy-preserving trigonometrically fitted continuous stage Runge-Kutta-Nyström methods for oscillatory Hamiltonian systems. Numer Algor 81, 1379–1401 (2019). https://doi.org/10.1007/s11075-019-00655-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s11075-019-00655-4

Keywords

Search

Navigation