Abstract
The main theme of this paper is explicit Gautschi-type integrators for the nonlinear multi-frequency oscillatory second-order initial value problems of the form \(y^{\prime \prime }= -A(t,y)y+f(t,y),\ y(t_{0}) = y_{0},\ y^{\prime }(t_{0}) = y^{\prime }_{0}\). This work is important and interesting within the broader framework of the subject. In fact, the Gautschi-type methods for oscillatory problems with a constant matrix A have been investigated by many authors. The key question now is that the classical variation-of-constants approach is not applicable to the oscillatory nonlinear problems with a variable coefficient matrix A(t,y). We consider successive approximations or locally equivalent systems for the problems, and derive efficient explicit Gautschi-type integrators. The error analysis is presented for the local approximation accordingly. Accompanying numerical results demonstrate the remarkable efficiency of the new Gautschi-type integrators in comparison with some existing numerical methods in the scientific literature.
Similar content being viewed by others
References
Gautschi, W.: Numerical integration of ordinary differential equations based on trigonometric polynomials. Numer. Math. 3, 381–397 (1961)
Hersch, J.: Contribution à la méthode des équations aux différences. ZAMP 9a, 129–180 (1958)
Hochbruck, M., Lubich, C. h.: On Krylov subspace approximations to the matrix exponential operator. SIAM J. Numer. Anal. 34, 1911–1925 (1997)
Hochbruck, M., Lubich, C. h., Selhofer, H.: Exponential integrators for large systems of differential equations. SIAM J. Sci. Comput. 19, 1552–1574 (1998)
Hochbruck, M., Lubich, C. h.: A Gautschi-type method for oscillatory second-order differential equations. Numer. Math. 83, 403–426 (1999)
Franco, J.: New methods for oscillatory systems based on ARKN methods. Appl. Numer. Math. 56, 1040–1053 (2006)
Wang, B., Iserles, A., Wu, X.: Arbitrary–order trigonometric Fourier collocation methods for multi-frequency oscillatory systems. Found. Comput. Math. 16, 151–181 (2016)
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)
Yang, H., Wu, X., You, X., Fang, Y.: Extended RKN-type methods for numerical integration of perturbed oscillators. Comput. Phys. Commun. 180, 1777–1794 (2009)
Wu, X., You, X., Li, J.: Note on derivation of order conditions for ARKN methods for perturbed oscillators. Comput. Phys. Commun. 180, 1545–1549 (2009)
Mei, L., Liu, C., Wu, X.: An essential extension of the finite-energy condition for extended Runge–Kutta–Nyström integrators when applied to nonlinear wave equations. Commun. Comput. Phys. 22, 742–764 (2017)
Wu, X., You, X., Xia, J.: Order conditions for ARKN methods solving oscillatory systems. Comput. Phys. Commun. 180, 2250–2257 (2009)
Wu, X., Wang, B.: Comments on “embedded pair of extended Runge-Kutta-Nystrom type methods for perturbed oscillators”. Appl. Math. Modell. 34, 3708–3711 (2010)
Wu, X., You, X., Shi, W., Wang, B.: ERKN integrators for systems of oscillatory second-order differential equations. Comput. Phys. Commun. 181, 1873–1887 (2010)
Wu, X., Wang, B.: Multidimensional adapted Runge-Kutta-Nyström methods for oscillatory systems. Comput. Phys. Commun. 181, 1955–1962 (2010)
Wu, X., Liu, K., Shi, W.: Structure-preserving algorithms for oscillatory differential equations, vol. II. Springer, Heidelberg (2015)
Grimm, V.: On error bounds for the Gautschi-type exponential integrator applied to oscillatory second-order differential equations. Numer. Math. 100, 71–89 (2005)
Hale, J. K.: Ordinary differential equations. Roberte E. Krieger Publishing Company, Huntington (1980)
Wu, X., Wang, B., Shi, W.: Effective integrators for nonlinear second-order oscillatory systems with a time-dependent frequency matrix. Appl. Math. Modell. 37, 6505–6518 (2013)
Wu, X., You, X., Wang, B.: Structure-preserving algorithms for oscillatory differential equations. Springer, Berlin (2013)
Iserles, A.: On the global error of discretization methods for highly oscillatory ordinary differential equations. BIT 42, 561–599 (2002)
Van de Vyver, H.: A symplectic exponentially fitted modified Runge-Kutta-Nyström method for the numerical integration of orbital problems. New Astron. 10, 261–269 (2005)
Hairer, E., Nørsett, S.P., Wanner, G.: Solving ordinary differential equations I: Nonstiff problems. Springer, Berlin (1993)
Coleman, J. P., Duxbury, S. C.: Mixed collocation methods for y ″ = f(x,y). J. Comput. Appl. Math. 126, 47–75 (2000)
Van de Vyver, H.: An explicit Numerov-type method for second-order differential equations with oscillating solutions. Computers and Mathematics with Applications 53, 1339–1348 (2007)
Magnus, W.: On the exponential solution of differential equations for a linear operator. Commun. Pure Appl. Math. 7, 649–673 (1954)
Iserles, A., Munthe-Kaas, H. Z., Nørsett, S.P., Zanna, A: Lie group methods. Acta Numerica 9, 215–365 (2000)
Nyström, E. J.: Ueber die numerische Integration von Differentialgleichungen. Acta. Soc. Sci. Fenn. 50, 1–54 (1925)
Funding
The research was financially supported in part by the Natural Science Foundation of China under Grants 11501288, 11701271, and 11671200, by the Specialized Research Foundation for the Doctoral Program of Higher Education under Grant 20100091110033, by the 985 Project at Nanjing University under Grant 9112020301, by the Natural Science Foundation of Jiangsu Province under Grant BK20150934, and by the Natural Science Foundation of the Jiangsu Higher Education Institutions under Grants 16KJB110010 and 14KJB110009.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Shi, W., Wu, X. Explicit Gautschi-type integrators for nonlinear multi-frequency oscillatory second-order initial value problems. Numer Algor 81, 1275–1294 (2019). https://doi.org/10.1007/s11075-018-0635-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-018-0635-7
Keywords
- Gautschi-type integrators
- Oscillatory nonlinear second-order initial value problems
- Structure preservation
- Error analysis
- Kepler’s problem
- Oscillatory Hamiltonian systems