Hybrid Numerov-Type Methods with Coefficients Trained to Perform Better on Classical Orbits

Chenglian Liu1, Chieh-Wen Hsu2, Ch. Tsitouras3, T. E. Simos4,5,6,7,8
1Department of Computer Science and Technology, Neusoft Institute Guangdong, Foshan, China
2Economics and Management College, Zhaoqing University, Zhaoqing, China
3General Department, National and Kapodistrian University of Athens, Athens, Greece
4College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao, China
5Department of Mathematics, College of Sciences, King Saud University, Riyadh , Saudi Arabia
6Group Leader of Modern Computational Methods, Ural Federal University, Yekaterinburg, Russian Federation
7Section of Mathematics, Department of Civil Engineering, Democritus University of Thrace, Xanthi, Greece
8Athens, Greece

Tóm tắt

This manuscript constructed a class of explicit hybrid Numerov methods of sixth order for second-order ordinary differential equations. The methods have four stages at each step with coefficients expressed explicitly with respect to a couple of free parameters. The coefficients are “trained” so that the methods perform best on a single Kepler orbit. In the experiment, the new methods outperform the ones in the literature for various eccentricities, costs and time intervals.

Tài liệu tham khảo

Hairer, E.: Unconditionally stable methods for second order differential equations. Numer. Math. 32, 373–379 (1979) Cash, J.R.: High order P-stable formulae for the numerical integration of periodic initial value problems. Numer. Math. 37, 355–370 (1981) Chawla, M.M.: Two-step fourth order P-stable methods for second order differential equations. BIT 21, 190–193 (1981) Chawla, M.M.: Numerov made explicit has better stability. BIT 24, 117–118 (1984) Tsitouras, Ch.: Explicit Numerov type methods with reduced number of stages. Comput. Math. Appl. 45, 37–42 (2003) Chawla, M.M., Rao, P.S.: An explicit sixth-order method with phase-lag of order eight for \(y^{\prime \prime }=f(t, y)\). J. Comput. Appl. Math. 17, 365–368 (1987) Tsitouras, Ch.: Explicit eighth order two-step methods with nine stages for integrating oscillatory problems. Int. J. Mod. Phys. C 17, 861–876 (2006) Tsitouras, Ch.: Explicit two step methods for second order linear IVPs. Comput. Math. Appl. 43, 943–949 (2002) Simos, T.E., Tsitouras, Ch.: Evolutionary generation of high order, explicit, two step methods for second order linear IVPs. Math. Methods Appl. Sci. 40, 6276–6284 (2017) Simos, T.E., Tsitouras, Ch.: A new family of seven stages, eighth order explicit Numerov-type methods. Math. Methods Appl. Sci. 40, 7867–7878 (2017) Simos, T.E., Tsitouras, Ch., Famelis, ITh: Explicit Numerov type methods with constant coefficients: a review. Appl. Comput. Math. 16, 89–113 (2017) Berg, D.B., Simos, T.E., Tsitouras, Ch.: Trigonometric fitted, eighth-order explicit Numerov-type methods. Math. Methods Appl. Sci. 41, 1845–1854 (2018) Franco, J.M.: A class of explicit two-step hybrid methods for second-order IVPs. J. Comput. Appl. Math. 187, 41–57 (2006) Franco, J.M., Randez, L.: Explicit exponentially fitted two-step hybrid methods of high order for second-order oscillatory IVPs. Appl. Math. Comput. 273, 493–505 (2016) Franco, J.M., Randez, L.: Eighth-order explicit two-step hybrid methods with symmetric nodes and weights for solving orbital and oscillatory IVPs. Int. J. Mod. Phys. C 29(1), 1850002 (2018) Famelis, ITh, Tsitouras, Ch.: Symbolic derivation of order conditions for hybrid Numerov-type methods solving \(y^{\prime \prime }=f(x, y)\). J. Comput. Appl. Math. 218, 543–555 (2008) Famelis, ITh, Papakostas, S.N., Tsitouras, Ch.: Symbolic derivation of Runge–Kutta order conditions. J. Symb. Comput. 37, 311–327 (2004) Tsitouras, Ch., Famelis, ITh: Symbolic derivation of Runge–Kutta–Nyström order conditions. J. Math. Chem. 46, 896–912 (2009) Simos, T.E.: High order closed Newton–Cotes trigonometrically-fitted formulae for the numerical solution of the Schrödinger equation. Appl. Math. Comput. 209, 137–151 (2009) Tsitouras, Ch., Papageorgiou, G.: Runge–Kutta interpolants based on values from two successive integration steps. Computing 43, 255–266 (1990) Papageorgiou, G., Tsitouras, Ch., Papakostas, S.N.: Runge–Kutta pairs for periodic initial value problems. Computing 51, 151–163 (1993) Monovasilis, Th, Kalogiratou, Z., Simos, T.E.: A family of trigonometrically fitted partitioned Runge–Kutta symplectic methods. Appl. Math. Comput. 209, 91–96 (2009) Tsitouras, C., Famelis, ITh, Simos, T.E.: Phase-fitted Runge–Kutta pairs of orders 8(7). J. Comput. Appl. Math. 321, 226–231 (2017) Simos, T.E., Tsitouras, Ch.: Fitted modifications of classical Runge–Kutta pairs of orders 5(4). Math. Methods Appl. Sci. 41, 4549–4559 (2018) Papadopoulos, D.F., Simos, T.E.: The use of phase lag and amplification error derivatives for the construction of a modified Runge–Kutta–Nyström method. Abstr. Appl. Anal. 2013, Article ID 910624 (2013) Hui, Fei, Simos, T.E.: Four stages symmetric two-step P-stable method with vanished phase-lag and its first, second, third and fourth derivatives. Appl. Comput. Math. 15, 220–238 (2016) Simos, T.E.: Multistage symmetric two-step P-stable method with vanished phase-lag and its first, second and third derivatives. Appl. Comput. Math. 14, 296–315 (2015) Zhang, Wei, Simos, T.E.: A high-order two-step phase-fitted method for the numerical solution of the Schrödinger equation. Mediterr. J. Math. 13, 5177–5194 (2016) Dong, Ming, Simos, T.E.: A new high algebraic order efficient finite difference method for the solution of the Schrödinger equation. Filomat 31, 4999–5012 (2017) Bratsos, A.G., Tsitouras, Ch., Natsis, D.G.: Linearized numerical schemes for the Boussinesq equation. Appl. Numer. Anal. Comput. Math. 2, 34–53 (2005) Alolyan, I., Simos, T.E.: A family of high-order multistep methods with vanished phase-lag and its derivatives for the numerical solution of the Schrödinger equation. Comput. Math. Appl. 62, 3756–3774 (2011) Simos, T.E.: New stable closed Newton–Cotes trigonometrically fitted formulae for long-time integration. Abstr. Appl. Anal. 2012, Article ID 182536 (2012) Alolyan, I., Anastassi, Z.A., Simos, T.E.: A new family of symmetric linear four-step methods for the efficient integration of the Schrödinger equation and related oscillatory problems. Appl. Math. Comput. 218, 5370–5382 (2012) Simos, T.E.: On the explicit four-step methods with vanished phase-lag and its first derivative. Appl. Math. Inf. Sci. 8, 447–458 (2014) Medvedev, Maxim A., Simos, T.E., Tsitouras, Ch.: Explicit, two stage, sixth order, hybrid four-step methods for solving \(\text{ y }^{\prime \prime }=\text{ f }(\text{ x },\text{ y })\). Math. Methods Appl. Sci. 41, 6997–7006 (2018) Simos, T.E., Tsitouras, Ch.: High phase-lag order, four-step methods for solving \(y^{\prime \prime }=f(x, y)\). Appl. Comput. Math. 17, 307–316 (2018) Simos, T.E.: Optimizing a hybrid two-step method for the numerical solution of the Schrödinger equation and related problems with respect to phase-lag. J. Appl. Math. 2012, Article ID 420387 (2012) Simos, T.E., Tsitouras, Ch.: Fitted modifications of Runge–Kutta pairs of orders 6(5). Math. Methods Appl. Sci. 41, 6184–6194 (2018) Medvedev, M.A., Simos, T.E., Tsitouras, Ch.: Trigonometric fitted hybrid four-step methods of sixth order for solving y\(^{\prime \prime }=\) f(x, y). Math. Methods Appl. Sci. 42, 710–716 (2019) Medvedev, M.A., Simos, T.E., Tsitouras, Ch.: Hybrid, phase-fitted, four-step methods of seventh order for solving \(\text{ y }^{\prime \prime }=\text{ f }(\text{ x },\text{ y })\). Math. Methods Appl. Sci. 42, 2025–2032 (2019) Tsitouras, Ch., Simos, T.E.: Trigonometric fitted explicit Numerov type method with vanishing phase-lag and its first and second derivatives. Mediterr. J. Math. 15, 168 (2018) Butcher, J.C.: Implicit Runge–Kutta processes. Math. Comput. 18, 50–64 (1964) Butcher, J.C.: On Runge–Kutta processes of high order. J. Aust. Math. Soc. 4, 179–194 (1994) Tsitouras, Ch., Simos, T.E.: On ninth order, explicit Numerov type methods with constant coefficients. Mediterr. J. Math. 15, Article ID 46 (2018) Coleman, J.P.: Order conditions for a class of two-step methods for \(y^{\prime \prime } = f(x, y)\). IMA J. Numer. Anal. 23, 197–220 (2003) Papadopoulos, D.F., Simos, T.E.: A modified Runge–Kutta–Nyström method by using phase lag properties for the numerical solution of orbital problems. Appl. Math. Inf. Sci. 7, 433–437 (2013) Panopoulos, G.A., Simos, T.E.: A new optimized symmetric embedded predictor–corrector method (EPCM) for initial-value problems with oscillatory solutions. Appl. Math. Inf. Sci. 8, 703–713 (2014) Panopoulos, G.A., Simos, T.E.: An eight-step semi-embedded predictor-corrector method for orbital problems and related IVPs with oscillatory solutions for which the frequency is unknown. J. Comput. Appl. Math. 290, 1–15 (2015) Ramos, H., Kalogiratou, Z., Monovasilis, Th, Simos, T.E.: An optimized two-step hybrid block method for solving general second order initial-value problems. Numer. Algorithms 72, 1089–1102 (2016) Kalogiratou, Z., Monovasilis, Th, Ramos, Higinio, Simos, T.E.: A new approach on the construction of trigonometrically fitted two step hybrid methods. J. Comput. Appl. Math. 303, 146–155 (2016) Storn, R., Price, K.: Differential evolution—a simple and efficient heuristic for global optimization over continuous spaces. J. Glob. Optim. 11, 341–359 (1997) Lin, C., Chen, J.J., Simos, T.E., Tsitouras, Ch.: Evolutionary derivation of sixth order P-stable SDIRKN methods for the solution of PDEs with the method of lines. Mediter. J. Math. 16, Article ID 69 (2019) Famelis, ITh, Alexandridis, A., Tsitouras, Ch.: High accuracy hybrid DE-PSO algorithm for the construction of Runge–Kutta pairs. Eng. Optim. 50, 1364–1379 (2017) Monovasilis, T., Kalogiratou, Z., Simos, T.E.: Construction of exponentially fitted symplectic Runge–Kutta–Nyström methods from partitioned Runge–Kutta methods. Mediterr. J. Math. 13, 2271–2285 (2016) Hairer, E., Nørsett, S.P., Wanner, G.: Solving Ordinary Differential Equations I, Nonstiff Problems, 2nd edn. Springer, Berlin (1993) Kierzenka, J., Magherini, C., Mazzia, F.: Test set for IVP solvers, release 2.4, Feb (2008). http://pitagora.dm.uniba.it/~testset/ Papakostas, S.N., Tsitouras, Ch.: High phase-lag order Runge–Kutta and Nyström pairs. SIAM J. Sci. Comput. 21, 747–763 (1999)